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Abstract 

Stochastic backgrounds or relic gravitons, if ever detected, will constitute a prima facie evi- 
dence of physical processes taking place during the earliest stages of the evolution of the plasma. 
The essentials of the stochastic backgrounds of relic gravitons are hereby introduced and reviewed. 
The pivotal observables customarily employed to infer the properties of the relic gravitons are 
discussed both in the framework of the ACDM paradigm as well as in neighboring contexts. The 
complementarity between experiments measuring the polarization of the Cosmic Microwave Back- 
ground (such as, for instance, WMAP, Capmap, Quad, Cbi, just to mention a few) and wide band 
interferometers (e.g. Virgo, Ligo, Geo, Tama) is emphasized. While the analysis of the microwave 
sky strongly constrains the low-frequency tail of the relic graviton spectrum, wide-band detectors 
are sensitive to much higher frequencies where the spectral energy density depends chiefly upon 
the (poorly known) rate of post-inflationary expansion. 
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1 The spectrum of the relic gravitons 

1.1 The frequencies and wavelengths of relic gravitons 

Terrestrial and satellite observations, scrutinizing the properties of the electromagnetic spectrum, are 
unable to test directly the evolution of the background geometry prior to photon decoupling. The 
redshift probed by Cosmic Microwave Background (CMB in what follows) observations is of the order 
of Zdec — 1087 and it roughly corresponds to the peak of the visibility function, i.e. when most of 
the CMB photons last scattered free electrons (and protons). After decoupling the ionization fraction 
drops; the photons follow null geodesies whose slight inhomogeneities can be directly connected with 
the fluctuations of the spatial curvature present before matter-radiation decoupling. 

The temperature of CMB photons is, today, of the order of 2.725 K. The same temperature at 
photon decoupling must have been of the order of about 2962 K, i.e. 0.25 eV. The CMB temperature 
increases linearly with the redshift: this fact may be tested empirically by observing at high redshifts 
clouds of chemical compounds (like CN) whose excited levels may be populated thanks to the higher 
value of the CMB temperature [UEj. 

The initial conditions for the processes leading to formation of CMB anisotropies are set well 
before matter radiation equality and right after neutrino decoupling (taking place for temperatures of 
the order of the MeV) whose associated redshift is around 10^°. The present knowledge of particle 
interactions up to energy scales of the order of 200 GeV certainly provides important (but still indirect) 
clues on the composition of the plasma. 

If ever detected, relic gravitons might provide direct informations on the evolution of the Hubble 
rate for much higher redshifts. In a rudimentary realization of the ACDM paradigm§|, the inflationary 
phase can be modeled in terms of the expanding branch of de-Sitter space. Assuming that, right after 
inflation, the Universe evolves adiabatically and is dominated by radiation, the redshift associated 
with the end of inflation can be approximately computed as 



where gp denotes the weighted □ number of relativistic degrees of freedom at the onset of the radiation 
dominated evolution and 106.75 corresponds to the value of the standard model of particle interactions. 
In Eq. (11.11) it has been also assumed, quite generally, that H ~ 10~^ Mp as implied by the CMB 
observations. For more accurate estimates the quasi-de Sitter nature of the inflationary expansion must 
be taken into account. In the ACDM paradigm, the basic mechanism responsible for the production 
of relic gravitons is the parametric excitation of the (tensor) modes of the background geometry and 
it is controlled by the rate of variation of space-time curvature. 

In the present article the ACDM paradigm will always be assumed as a starting point for any 
supplementary considerations. The reasons for this choice are also practical since the experimental 
results must always be stated and presented in terms of a given reference model. Having said this, 

^Natural units h — c — = I will be adopted. In this system, K = 8.617 x lO^'^eV. 

^In the acronym ACDM, A qualifies the dark-energy component while CDM qualifies the dark matter component. 

■^Fermions and bosons contribute with different factors to gp. By assuming that all the species of the standard model 
are in local thermodynamic equilibrium (for instance for temperature higher than the top quark mass), gp will be given 
by gp — 28 + (7/8)90 = 106.75 where 28 and 90 count, respectively, the bosonic and the fermionic contributions. 




(1.1) 
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most of the considerations presented here can also be translated (with the appropriate computational 
effort) to different models. 

Given a specific scenario for the evolution of the Universe (like the ACDM model), the relic graviton 
spectra can be computed. The amplitude of the relic graviton spectrum over different frequencies 
depends upon the specific evolution of the Hubble rate. The theoretical error on the amplitude increases 
with the frequency: it is more uncertain (even within a specified scenario) at high frequencies rather 
than at small frequencies. 

The experimental data, at the moment, do not allow either to rule in or to rule out the presence of 
a primordial spectrum of relic gravitons compatible with the ACDM scenario. The typical frequency 
probed by CMB experiments is of the order of § i/p = /Cp/(27r) 10^^^ Hz = 1 aHz where fcp is the 
pivot frequency at which the tensor power spectra are assigned. CMB experiments will presumably 
set stronger bounds on the putative presence of a tensor background for frequencies (9(aHz). This 
bound will be significant also for higher frequencies only if the whole post-inflationary thermal history 
is assumed to be known and specified. The typical frequency window of wide-band interferometers 

The relic graviton spectrum The eiectromagnetic spectrum 




log(v/Hz) log(v/Hz) 

Figure 1: A schematic view of the frequency range of the relic graviton spectrum (plot at the left) 
and of the electromagnetic spectrum (plot at the right). In both plots the common logarithms of the 
(comoving) frequency and of the (comoving) wavelengths are reported, respectively, on the horizontal 
and on the vertical axis. 

(such as Ligo and Virgo) is located between few Hz and 10 kHz, i.e. roughly speaking, 20 orders of 
magnitude larger than the frequency probed by CMB experiments. The frequency range of wide-band 
interferometers will be conventionally denoted by z/lv To compute the relic graviton spectrum over 
the latter range of frequencies, the evolution of our Universe should be known over a broad range of 
redshifts. We do have some plausible guesses on the evolution of the plasma from the epoch of neutrino 
decoupling down to the epoch of photon decoupling. The latter range of redshifts corresponds to an 
interval of comoving frequencies going from z/p ~ 10^^^ Hz up to z/bbn — 0.01 nHz ~ 10^^^ Hz (at 
most). 

In Fig. [1] (plot at the left) the frequency range of the relic graviton spectrum is illustrated, starting 
from the (comoving) frequency Up whose associated (comoving) wavelength is of the order of of 10^^ 

^We are here enforcing the usual terminology stemming from the prefixes of the International System of units: aHz 
(for atto Hz i.e. lO^^'* Hz), fHz (for femto Hz, i.e. 10^^^ Hz) and so on. 
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m, i.e. roughly comparable with the Hubble radius at the present time. The low-frequency branch 
of the spectrum can be conventionally defined between Up and I'hhn- The largest frequency of the 
relic graviton spectrum (i.e. z/max) is of the order of 0.1 GHz in the ACDM scenario. Thus the high- 
frequency branch of the graviton spectrum can be conventionally defined for i/bbn < z^ < t'max- In 
summary we can therefore say that 

• the range z/p < z/ < z/bbn "will be generically referred to as the low-frequency domain; in this range 
the spectrum of relic gravitons basically follows from the minimal ACDM paradigm; 

• the range z/bbn < z/ < z/max will be generically referred to as the high-frequency domain; in this 
range the spectrum of relic gravitons is more uncertain. 

The high-frequency branch of the relic graviton spectrum, overlapping with the frequency window of 
wide-band detectors (see shaded box in the left plot of Fig. [1]), is rather sensitive to the thermodynamic 
history of the plasma after inflation as well as, for instance, to the specific features of the underlying 
gravity theory at small scales. This is why we said that the theoretical error in the calculation of the 
relevant observables increases, so to speak, with the frequency. 

In Fig. [1] (plot at the right) the electromagnetic spectrum is reported in its salient features. It 
seems instructive to draw a simple minded parallel between the electromagnetic spectrum and the 
spectrum of relic gravitons. Consider first the spectrum of relic gravitons (see Fig. [1], plot at the 
left): between 10~^^ Hz (corresponding to z/p) and 10 kHz (corresponding to z/lv) there are, roughly, 
22 decades in frequency. A similar frequency gap (see Fig. [1], plot at the right), if applied to the 
well known electromagnetic spectrum, would drive us from low-frequency radio waves up to x-rays or 
7-rays. As the physics explored by radio waves is very different from the physics probed by 7 rays, 
it can be argued that the informations carried by low and high frequency gravitons originate from 
two very different physical regimes of the theory. Low frequency gravitons are sensitive to the large 
scale features of the given cosmological model and of the underlying theory of gravity. High frequency 
gravitons are sensitive to the small scale features of a given cosmological model and of the underlying 
theory of gravity. 

The interplay between long wavelength gravitons and CMB experiments will be specifically dis- 
cussed in the subsection 11.21 The main message will be that, according to current CMB experiments, 
long wavelength gravitons have not been observed yet. The latter occurrence imposes a very important 
constraint on the low frequency branch of the relic graviton spectrum of the ACDM scenario whose 
salient predictions will be introduced in subsection 11.31 According to the minimal ACDM paradigm 
a very peculiar conclusion seems to pop up: the CMB constraints on the low-frequency tail of the 
graviton spectrum jeopardize the possibility of any detectable signal for frequencies comparable with 
the window explored by wide band interferometers (see subsection ll.4p . The natural question arising 
at this point is rather simple: is it possible to have a quasi-flat low-frequency branch of the relic 
graviton spectrum and a sharply increasing spectral energy density at high-frequencies? This kind of 
signal is typical of a class of completions of the ACDM paradigm which have been recently dubbed 
TACDM (for tensor -ACDM). The main predictions of these models will be introduced in subsection 
II. 5[ We shall conclude this introductory section with a discussion of two relevant constraints which 
should be applied to relic graviton backgrounds in general, i.e. the millisecond pulsar and the big-bang 
nucleosynthesis constraint (see subsection II. 6p . 
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1.2 Long wavelength gravitons and CMB experiments 

The bounds on the backgrounds of rehc gravitons stemming from CMB experiments are phrased in 
terms of which is the ratio between the tensor and the scalar power spectra at the same conventional 
scale (often called pivot scale). While the use of tt is practical (see e. g. the 5- year WMAP data 
[21 m El El E]), it assumes the ACDM scenario insofar as the curvature perturbations are adiabatic. 
Within the ACDM model, the tensor and scalar power spectra can be parametrized as 

At = r^An, (1.2) 

^7^= (2.41 ±0.11) X 10"^ A;p = 0.002 Mpc"^ (1.3) 

where kp is the pivot wave-number and A-ji is the amplitude of the power spectrum of curvature 
perturbations!^ computed at kp] Ug and nx are, respectively, the scalar and the tensor spectral indices 
l|. The value of kp is conventional and it corresponds to an effective harmonic ics — 30. The figure 
for A-ji quoted in Eq. (11.31) corresponds to the value inferred from the WMAP 5-year data [31 IH 
El El IZ] in combination with the minimal ACDM models. In the ACDM model the origin of Atz 
stems from adiabatic curvature perturbations which are present after neutrino decoupling but before 
matter radiation equahty (taking place at a redshift = 3176li\ll according to the WMAP 5-yr data 
[SmiE])- The dominant component of curvature perturbations is adiabatic meaning that, over large 
scales, the fluctuations in the specific entropy are vanishing, at least in the minimal version of the 
model. The adiabatic nature of the fluctuations induces a simple relation between the first acoustic 
peak of the TT power spectra and the first anticorrelation peak of the TE power spectra f9j: this is, 
to date, the best evidence that curvature perturbations are, predominantly, adiabatic. It is useful to 
translate the comoving wave number kp into a comoving frequency 

k 

Vp = = 3.092 X 10"^^ Hz = 3.092 aHz, (1.4) 

so, as anticipated, z/p is of the order of the aHz. The amplitude at the pivot scal^ is controlled 
exactly by r^- The combined analysis of the CMB data, of the large-scale (LSS) structure data 
[T3l [H] and of the supernova (SN) data [151 HH] can lead to quantitative upper limits on rx which 
are illustrated in Tabs. [T] and [2] as they are emerge from the combined analyses of different data sets. 
The inferred values of the scalar spectral index (i.e. ng), of the dark energy and dark matter fractions 
(i.e., respectively, VL^ and f^Mo); and of the typical wavenumber of equality fceq are reported in the 
remaining columns. While different analyses can be performed, it is clear, by looking at Tabs. [1] and 
[21 that the typical upper bounds on r^ikp) range between, say, 0.2 and 0.4. More stringent limits 
can be obtained by adding supplementary assumptions. In Tab. [T] the quantity as determines the 
frequency dependence of the scalar spectral index. In the simplest case as = and the spectral index 

^The perturbations of the spatial curvature, conventionally denoted by TZ are customarily employed to characterize 
the scalar fluctuations of the geometry since TZ is approximately constant (in time) across the radiation- matter transition. 

^As it is clear from Eqs. (|1.2|) and (|1.3|) there is a difference in the way the scalar and the tensor spectral indices are 
assigned: while the scale-invariant limit corresponds to rig ^ 1 for the curvature perturbations, the scale invariant limit 
for the long wavelength gravitons corresponds to tit ~* 0. 

^See also [imUl mill] for earlier WMAP data releases. 

^In the first release of the WMAP data the scalar and tensor pivot scales were chosen to be different and, in 
particular, = 0.05 Mpc~^ for the scalar modes. In the subsequent releases of data the two pivot scales have been 
taken to coincide. 



PT(fc)=^T(^ 



Vn{k) = An 



k ^ "-^1 
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Data 


rx 


Us 


as 




^MO 


A;cqMpc 


WMAP5 alone 
WMAP5 + Acbar 
WMAP5+ LSS + SN 
WMAP5+ CMB data 


< 0.58 

< 0.54 

< 0.54 

< 0.64 


1.087°:[]^o^73 

±.UOO_Q 062 

-L-Uyo_o.069 
1 1 270-075 
-0.071 


-0.050 ±0.034 
-0.048 ±0.027 

—0 05^0.027 

-0 072+° °^^ 

U.U(Z_oo30 


n 799+0.054 
U. /ZZ_o.057 

0.7191°:°^? 

0.7141°:°!^ 

0.7040j055^, 


0.281°_;°:'ols 
0.2860_;0i6^7 
0.2960_;05o'55 


0.00993 
0.01004 
0.01008 
0.01013 



Table 1: The change in determination of the parameters of the tensor background for three different 
choices of cosmological data sets. 



Data 






as 




^MO 


/^cqMpC 


WMAP5 alone 


< 0.43 


0.986 ±0.22 





770+0-033 

U- ' ' U_o.032 


Q 2SnO-032 
U.ZOU_o 033 


0.00936 


WMAP5 ± Acbar 


< 0.40 


n qfic;0.019 

u.yoj_o 020 





0.767 ±0.032 


0.233 ±0.032 


0.00944 


WMAP5± LSS ± SN 


< 0.20 


0.968 ±0.015 





0.725 ±0.015 


0.275 ±0.015 


0.00999 


WMAP5± CMB data 


< 0.36 


0.979 ±0.020 





0.775 ±0.032 


0.225 ±0.032 


0.00922 



Table 2: Same as in Tab. [T]but assuming no running in the (scalar) spectral index (i.e. as = 0). 



is frequency- independent (i.e. ng does not run with the frequency). It can also happen, however, that 
as ^ which implies an effective frequency dependence of the spectral index. If the inflationary phase 
is driven by a single scalar degree of freedom (as contemplated in the minimal version of the ACDM 
scenario) and if the radiation dominance kicks in almost suddenly after inflation, the whole tensor 
contribution can be solely parametrized in terms of r^. The rationale for the latter statement is that 
tt not only determines the tensor amplitude but also, thanks to the algebra obeyed by the slow-roll 
parameters, the slope of the tensor power spectrum, customarily denoted by nx- To lowest order in 
the slow-roll expansion, therefore, the tensor spectral index is slightly red and it is related to (and 
to the slow-roll parameter) as 



riT 



' 8 



-2e, 



H 



(1.5) 



where e measures the rate of decrease of the Hubble parameter during the inflationary epoch 0. Within 
the established set of conventions the scalar spectral index ng is given by ns = (1 — 6e ± 27]) and it 
depends not only upon e but also upon the second slow-roll parameter rj = MpV^^^/V (where V is the 
inflaton potential, V^^ denotes the second derivative of the potential with respect to the inflaton field 
and Mp = 1/ VSirG) . It is sometimes assumed that also is not constant but it is rather a function 
of the wavenumber, i.e. 



nT{k) = -2e + —In (k/kp), ax = — 

Z o 



1) + ^ 



;i.6) 



where ax now measures the running of the tensor spectral index. 

As already mentioned, among the CMB experiments a central role is played by WMAP [31 HI O [6l [7] 
(see also [H [U [10] for first year data release and [HI [12] for the third year data release. In connection 
with [31 HI [SI [SI [Z] , the WMAP 5-year data have been also combined with observations of the Acbar 



^"^The overdot will denote throughout the paper a derivation with respect to the cosmic time coordinate t while the 
prime will denote a derivation with respect to the conformal time coordinate r. 



7 



experiment [TTJ [HI [T9l [20] . The TT, TE and, partially EE angular power spectra^ have been 
measured by the WMAP experiment. Other (i.e. non space-borne) experiments are now measuring 
polarization observables, in particular there are 

• the Dasi (degree angular scale interferometer) experiment [211 ESI EH] operating at south pole ; 

• the Capmap (cosmic anisotropy polarization mapper) experiment pll [25] ; 

• the Cbi (cosmic background imager) experiment [261 [2Z] ; 

• the Quad experiment [28] [291 [30] : 

as well as various other experiments at different stages of development In the near future the Planck 
explorer satellite [31] might be able to set more direct limits on tt by measuring (hopefully) the BB 
angular power spectra. 



1.3 The relic graviton spectrum in the ACDM model 

Having defined the frequency range of the spectrum of relic gravitons, it is now appropriate to illustrate 
the possible signal which is expected within the ACDM scenario. In Fig. [2] the spectrum of the relic 
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Figure 2: The spectral energy density of the relic gravitons is illustrated the context of the ACDM 
scenario and for different values of rx (see Eqs. fll.2p - fll.3p ). In the plot at the right the tensor 
spectral index ut does depend upon the frequency (q;t follows from Eq. (11.61) with the ACDM choice 
of parameters). The plot at the left corresponds to ax = 0. 



gravitons is reported in the case of the minimal ACDM scenario for different values of rx. On the 

^^The Arcminute Cosmology Bolometer Array Receiver (ACBAR) operates in three frequencies, i.e. 150, 219 and 274 
GHz. 

^^FoUowing the custom the TT correlations will simply denote the angular power spectra of the temperature autocor- 
relations. The TE and the EE power spectra denote, respectively, the cross power spectrum between temperature and 
polarization and the polarization autocorrelations. 

^•^Other planned experiments have, as specific target, the polarization of the CMB. In particular it is worth quoting 
here the recent projects Clover [35], Brain [33], Quiet [31] and Spider f35] just to mention a few. 
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horizontal axis the common logarithm of the comoving frequency is reported. The spectral energy 
density per logarithmic interval of frequency and in critical units is illustrated on the vertical axis. 
More quantitatively Qg^^{i',to) is defined as0 

^Gw('^,i"o) = J—^^^ (1.7) 
Pcrit dlnz/ 

2 

where pcrit = SifgMp is the critical energy density. Since pcrit depends upon Hq (i.e. the present 
value of the Hubble rate), it is practical to plot directly /iq ^gw('^5 To) at the present (conformal) time 
tq. The proper definition of f2Gw('^)To) in terms of the energy-momentum pseudo-tensor in curved 
space-time is postponed to section O The salient features of the relic graviton spectra arising in the 
context of the ACDM scenario can be appreciated by looking carefully at Fig. [2l 

The infra-red branch of the relic graviton spectrum (see also Fig. [1]) extends, approximately, from 
i/p up to a new frequency scale which can be numerically determined by integrating the evolution 
equations of the tensor modes and of the background geometry across the matter-radiation transition. 
A semi-analytic estimate of this frequency is given by 



= 1.281 X 10-17 (^^^V ^o^Ro y'/' 

VO.1326/ V4.15 X 10-V 



27r V0.1326A4.15 X 10-5. 

At intermediate frequencies Fig. [2] exhibits a further suppression which is due to the coupling of the 
tensor modes with the anisotropic stress provided by the collisionless species which are present prior to 
matter-radiation equality. This aspect has been recently emphasized in Ref. [35] (see also [2ZIES1ES])- 
Figure [2] assumes that the only colisionless species are provided by massless neutrinos, as the ACDM 
model stipulates and this corresponds, as indicated, to = 0.405. The quantity R^ measures the 
contribution of N^, families of massless neutrinos to the radiation plasma: 

r / N \ 

Ru = — r, = 0.681 , R^ + Ru = l. (1.9) 



The frequency range of the suppression due to neutrino free- streaming extends from z/gq up to fbbn 
which is given, approximately, by 

z^bbn = 2.252 X 10-" ( ( ^\ [ \ Hz ^ 0.01 nHz. (1. 10) 

V10.75/ V MeV/ V4.15 X 10-V ^ ' 

Both in Eqs. ( 11. 8p and ( ll.lOp VtyiQ and ^rq denote, respectively, the present critical fraction of matter 
and radiation with typical values drawn from the best fit to the WMAP 5-yr data alone and within 
the ACDM paradigm. In Eq. fll.lOp Qp denotes the effective number of relativistic degrees of freedom 
entering the total energy density of the plasma. While i^eq is still close to the aHz, z/bbn is rather in the 
nHz range. In Fig. [2] (plot at the left) the spectral index ut is frequency independent; in the plot at 
the right, always in Fig. [21 the spectral index does depend on the wavenumber. These two possibilities 
correspond, respectively, to ctx = and ctT 7^ in Eq. fll.6p . In the regime z/ < i/gq a numerical 
calculation of the transfer function is mandatory for a correct evaluation of the spectral slope. In 
the approximation of a sudden transition between the radiation and matter-dominated regimes the 
spectral energy density goes, approximately, as z/"^"*""'^. The spectra illustrated have been computed 
within the approach developed in [101 SI] and include also other two effects which can suppress the 
amplitude of the quasi-flat plateau. These effects are related to the contribution of the dark energy 



*In the present review the In wiU denote the natural logarithm while the log will always denote the common logarithm. 
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and to the evolution of the effective number of relativistic species. From a quantitative point of view 
both effects are, however, less relevant than neutrino free streaming. 

Apart from the modification induced by the neutrino free-streaming the slope of the spectral energy 
density for z/ > u^q is quasi flat and it is determined by the wavelengths which reentered the Hubble 
radius during the radiation-dominated stage of expansion. The suggestion that relic gravitons can 
be produced in isotropic Friedmann- Robert son- Walker models is due to Ref. [12] (see also [13]) and 
was formulated before the inflationary paradigm. After the formulation of the inflationary scenario 
the focus has been to compute reliably the low frequency branch of the relic graviton spectrum. In 
[HillSlllG] the low-frequency branch of the spectrum has been computed with slightly different analytic 
approaches but always assuming an exact de Sitter stage of expansion prior to the radiation-dominated 
phase. The analytical calculation (whose details will be described in section |6]) shows that in the range 
Up < u < z/gq, the spectral energy density of the rehc gravitons (see Eq. (11. 7p ) should approximately 
go as Qq-w^UjTq) ~ Within the same approximation, for z/ > z/gq the spectral energy density 

is exactly fiat (i.e. flGw{^,To) — i^^)- This result, obtainable by means of analytic calculations (see 
also [171 HHl HHl [50]), is a bit crude in the hght of more recent developments. To assess the accurately 
spectral energy density it is necessary to take into account that the infrared branch is gradually passing 
from a quasi-flat slope (for z/ > z/gq) to the slope z/~^ which is the one computed within the sudden 
approximation [Ul [13 [13 [50]. It is useful to quote some of the previous reviews which covered, in a 
more dedicated perspective, the subject of the stochastic backgrounds of relic gravitons. The review 
article by Thorne [5l] does not deal solely with relic graviton backgrounds while the reviews of Refs. 
[521 [S31 [SI] are more topical. 

The flat plateau of the spectral energy density extends, approximately, between z/gq and a certain 
J^max- Also the maximal amplified frequency can be computed once the model of smooth transition 
between inflation and radiation is known. The smoothness of the transition determines specifically the 
precise amount of exponential suppression for u > z/max- A simple estimate of z/max is given by 

= 0.346 (^) "7—^) ( ) GHz, (1.11) 

VO.Ol/ V2.41 X 10-V V4.15 x IQ-y ' ^ ^ 

where, as in Eqs. ( II. 2p and (11. 3|) . An denotes the amplitude of the power spectrum of curvature 
perturbations evaluated at the pivot wavenumber z/p. It is worth noticing that between z/bbn and z^max 
there are approximately 20 orders of magnitude in frequency. In the ACDM scenario the spectrum 
has, in this range, always the same slope (i.e. is frequency- independent in Eq. f ll.2l) ). 

Some details of the calculations leading to the spectral energy densities illustrated in Fig. [2] can be 
found in sections [5] and [B Without dwelling on the details it is however clear, as anticipated, that the 
constraints on the long wavelength gravitons make it difficult (if not impossible) to have a detectable 
spectral energy density at the scale of wide-band interferometers. The latter statement, valid in the 
minimal ACDM scenario, will be sharpened in the following subsection. 

1.4 Short wavelength gravitons and wide-band interferometers 

In the ACDM scenario the spectral energy density of the relic gravitons has its larger amplitude in 
the low-frequency branch. As the frequency increases the spectral energy density diminishes so that 
it is plausible to expect a rather small amplitude over the frequencies corresponding to wide-band 
interferometers (see, for instance. Fig. [2]for z/ ~ z/lv = 100 Hz). 
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Wide-band interferometers operate in a window ranging from few Hz up to 10 kHz (see also 
Fig. [Tj). The available interferometers are Ligo [55], Virgo [56], Tama [57| and Geo [58]. In loose 
terms these instruments are Michelson interferometers with two important differences: the mirrors are 
suspended and Fabry-Perot cavities are used to increase the optical path of the photons. It would 
be too pretentious to describe in detail, in the present script, also the experimental apparatus and 
we therefore suggest Ref. [59] where the basics of wide-band interferometers are introduced in a 
self-contined perspective. 

The sensitivity of a given pair of wide-band detectors to a stochastic background of relic gravitons 
depends upon the relative orientation of the instruments. The wideness of the band (important for 
the correlation among different instruments) is not as large as 10 kHz but typically narrower and, in 
an optimistic perspective, it could range up to 100 Hz. The putative frequency of wide-band detectors 
will therefore be indicated as u^y, i.e. in loose terms, the Ligo/Virgo frequency. There are daring 
projects of wide-band detectors in space like the Lisa [60], the Bbo [61] and the Decigo [62] projects. 
The common feature of these three projects is that they are all space-borne missions and that they are 
all sensitive to frequencies smaller than the mHz (1 mHz = 10~^ Hz). While wide-band interferometers 
are now operating and might even reach their advanced sensitivities during the incoming decade, the 
wished sensitivities of space-borne interferometers are still on the edge of the achievable technologies. 

Since z/bbn < ^lv < J^max, wide-band interferometers represent an ideal instrument to investigate 
the relic graviton spectrum at high-frequencies. The spectral energy density of the relic gravitons 
produced within the ACDM model is quite minute and it is undetectable by interferometers even in 
their advanced version where the sensitivity is expected to improve by 5 or even 6 orders of magnitude 
in comparison with the present performances [B31[Ml[SS| (see also [BB] and ^7\). In Fig. [3] the spectral 
energy density is reported for u = vym and always in the case of the prediction stemming from the 
minimal ACDM scenario. In Fig. [3l the common logarithm of the spectral energy density is illustrated 



h^n, =0.1326, R =0.405 

MO V 



n =0.1326, n = 0.963, R =0.405 

MO s V 





Figure 3: The spectral energy density of the relic gravitons in the context of the ACDM model 
evaluated at the Ligo/Virgo frequency as a function of the tensor-to-scalar ratio. To guide the intuition 
consider that (advanced) wide-band interferometers will achieve, optimistically, a sensitivity of 10~^° 
in VtQY4{v,TQ). Taking into account the current bounds on rx (see Tabs. [T]and[2]) it is clear that the 
ACDM scenario is definitively too small to provide a serious target for wide-band detectors. 



function of the common logarithm of Tt. 
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In Ref. [6lj (see also [631 ES]) the current limits on the presence of an isotropic background of 
relic gravitons have been assessed. According to the Ligo collaboration (see Eq. (19) of Ref. [M] ) 
the spectral energy density of a putative (isotropic) background of relic gravitons can be parametrized 



The parametrization of Eq. (11.121) fits very well with Fig. [3] where the pivot frequency z/lv = lOOHz 
coincides with the pivot frequency appearing in the parametrization (I1.12p . For the scale-invariant 
case (i.e. P = —3 in eq. (I1.12p ) the Ligo collaboration sets a 90% upper limit of 1.20 x 10~^ on the 
amplitude appearing in Eq. (11.121) . i.e. ficw.-a- Using different sets of data (see [S31E5]) the Ligo 
collaboration manages to improve the bound even by a factor 2 getting down to 6.5 x 10~^. Thus Fig. 
[3] together with the upper limit of Eq. (I1.12p shows that the current Ligo sensitivity is still too small 
to detect the relic graviton background arising within the ACDM paradigm. 

1.5 Beyond the ACDM paradigm and high-frequency gravitons 

In the case of an exactly scale invariant spectrum the correlation of the two (coaligned) LIGO detectors 
with central corner stations in Livingston (Lousiana) and in Hanford (Washington) might reach a 
sensitivity to a flat spectrum which is [681 EH EQ] 



where T denotes the observation time and SNR is the signal to noise ratio. Equation (11.131) is in 
close agreement with the sensitivity of the advanced Ligo apparatus [55] to an exactly scale-invariant 
spectral energy density [Til [721 [731 El] • Equation (11.131) together with the plots of Fig. [3] suggest 
that the relic graviton background predicted by the ACDM paradigm is not directly observable by 
wide-band interferometers in their advanced version. 

CMB observations probe the aHz region of the spectral energy density of Fig. [21 Wide-band 
interferometers probe a frequency range between few Hz and 10 kHz. In both ranges, the signal of 
the ACDM scenario might be too small to be directly detectable. In Fig. [H the spectral energy 
density is computed in an extension of the ACDM paradigm which has been dubbed tensor-ACDM 
(TACDM for short) [lOlllT]. In the TACDM scenario the transition from the inflationary phase to the 
radiation-dominated epoch is mediated by a rather long stiff phase. By stiff phase we mean a phase 
where the total sound speed of the plasma is larger than the sound speed of a radiation-dominated 
plasma (i.e. 1/v^ in natural units). In the simplest realization of the scenario the barotropic index 
Wt is constant during the stiff phase. For instance, in Fig. [H the cases Wt = I and Wt = 0.6 have 
been illustrated. Since Wt denotes the ratio between the total pressure and the total energy density, 
it is rather plausible to demand that Wt < 1. The latter requirement implies that the sound speed is 
always smaller than the speed of light. As suggested in [75] (see also [761 EZj) the presence of a stiff 
phase can have the effect of increasing the spectral energy density at high frequencies. The increase 

The variable /3 is used in Eq. ()1.12p just because this is the notation endorsed by the Ligo coUaboration and there 
is no reason to change it. At the same time, in the present review, f3 wiU be used also with different meanings. In section 
[SI (3 quantifies the theoretical error on the maximal frequency of the relic graviton spectrum(see e.g. Eq. (j6.48p and 
discussion therein). In section [7|/3 parametrizes a portion of the azimuthal structure of the Stokes parameters. Since 
none of these variables appear in the same context, potential clashes of conventions are avoided. 





(1.12) 




i/Lv = 0.1 kHz 



(1.13) 
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Figure 4: The spectral energy density of relic gravitons in the minimal TACDM scenario where the 
post-inflationary thermal history contemplates a phase where the sound speed of the plasma Cst is close 
to the speed of light c. In the plot at the right Cst = c (corresponding to a constant barotropic index 
Wt = 1); in the plot at the left Cgt = 0.77c (corresponding to a constant barotropic index Wt = 0.6). 
In both plots the other cosmological parameters have been chosen to coincide with the best fit to the 
WMAP 5-yr data alone. 

takes place for frequencies larger than the mHz and is typically maximal in the GHz region. The 
spectral energy densities illustrated in Fig. H] suggest that it is not impossible to imagine situations 
where the spectral energy density of the relic gravitons satisfies all the constraints demanded by CMB 
physics but, at the same time, it is sufficiently large to be observed by wide-band interferometers. The 
results reported in Fig. H] refer to the minimal TACDM model where only one post-inflationary phase 
stiffer than radiation is contemplated. The barotropic index could have, however, a more complicated 
dependence. Already in the examples of Fig. H] the numerical integration implies that the barotropic 
index does depend, effectively, upon the scale factor (see, e. g., discussions in section [6] on the transfer 
function for the spectral energy density). 

The comparison of Fig. [2] and H] suggests, in short, the following subjects of reflection: 

• the theoretical error in the estimate of the spectral energy density increases with the frequency; 

• departures from the standard post-inflationary thermal history can be directly imprinted in the 
primordial spectrum of the relic gravitons; 

• in the incoming decade the observations of wide-band interferometers could be analyzed in con- 
junction with more standard data sets (i.e. CMB data supplemented by large-scale structure 
data and by the observations of type la supernovae) to constrain the spectral energy density of 
the relic gravitons both at small and at high frequencies. 



The presence of post-inflationary phases stiffer than radiation is, after all, rather natural and this 
was the original spirit of [75]. We do not know which was the rate of the post-inflationary expansion 
and since guesses cannot substitute experiments it would be productive to use the TACDM paradigm 
as reference model for a unifled analysis of the low-frequency data stemming from CMB and of the 
high-frequency data provided by wide-band interferometers. Already in [75] (see also [761 [77]) ^ rather 
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special candidate for a post-inflationary phase stiffer than radiation was the case when the sound speed 
equals the speed of light, i.e. the case when the energy density of the sources driving the geometry 
is dominated by the kinetic term of a (minimally coupled) scalar field. This particular case was also 
prompted by various classes of quintessence models. A specific example of this dynamics was provided 
in [78]. 

A more detailed account of the techniques leading to Fig. H] will be swiftly presented in section 
[6] and can be found in [lOl HI]. Without going through the details it is however important to stress 
that the calculations should be accurate enough not only in the high-frequency region but also in the 
low-frequency part of the spectrum. Indeed, as stressed above, one of the purposes of the TACDM 
scenario is to convey the idea that low-frequency and high-frequency measurements of the relic graviton 
background can be analyzed in a single theoretical framework. 

1.6 The millisecond pulsar bound and the nucleosynthesis constraint 

The spectral energy density of the relic gravitons must be compatible not only with the CMB con- 
straints (bounding, from above, the value of tt) but also with the pulsar timing bound[79l [80] and the 
big-bang nucleosynthesis constraints [8T[ [821 [83] . The pulsar timing bound demands 

^^('^pulsar, T-o) < 1.9 X 10"^ Z/pu^ar ^ 10 uHz, (1.14) 

where z/puisar roughly corresponds to the inverse of the observation time during which the pulsars timing 
has been monitored. The spectral energy densities illustrated in Figs. [2] and [4] satisfy the pulsar timing 
bound. 

The most constraining bound for the high-frequency branch of the relic graviton spectrum is repre- 
sented by big-bang nucleosynthesis. Gravitons, being relativistic, can potentially increase the expan- 
sion rate at the BBN epoch. The increase in the expansion rate will affect, in particular, the synthesis 
of ^He. To avoid the overproduction of ^He the expansion rate the number of relativistic species must 
be bounded from above. 

The BBN bound is customarily expressed in terms of (equivalent) extra fermionic species. During 
the radiation-dominated era, the energy density of the plasma can be written as pt = fi'p(7i"^/30)T^ where 
T denotes here the common (thermodynamic) temperature of the various species. An (ultra) relativistic 
fermion species with two internal degrees of freedom and in thermal equilibrium contributes 2 ■ 7/8 = 
7/4 = 1.75 to Qp. Before neutrino decoupling the contributing relativistic particles are photons, 
electrons, positrons, and A^^i, = 3 species of neutrinos, giving 

= y + ^iV, = 10.75. (1.15) 

The neutrinos have decoupled before electron-positron annihilation so that they do not contribute 
to the entropy released in the annihilation. While they are relativistic, the neutrinos still retain 
an equilibrium energy distribution, but after the annihilation their (kinetic) temperature is lower, 
n = (4/11)V3t. Thus 

gp = 2 + ^N,(^^'^ = 2 + 0.454Ar, = 3.36, (1.16) 

after electron-positron annihilation. By now assuming that there are some additional relativistic 
degrees of freedom, which also have decoupled by the time of electron-positron annihilation, or just 
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some additional component px to the energy density with a radiation-hke equation of state (i.e. 
Px = Px/^)i the effect on the expansion rate will be the same as that of having some (perhaps a 
fractional number of) additional neutrino species. Thus its contribution can be represented by replacing 
Ny with Ny + ANy in the above. Before electron-positron annihilation we have px = {7/8)AN^p^ and 
after electron-positron annihilation we have px = {7 /8){4:/ll)^^^ AN^ p^ ^ 0.227 AA^j^p^. The critical 
fraction of CMB photons can be directly computed from the value of the CMB temperature and it is 
notoriously given by h^Q^ = p-y/ Pcrit = 2.47 x 10~^. If the extra energy density component has stayed 
radiation-like until today, its ratio to the critical density, Qx, is given by 

ftJn,.ftS^^5,61xlO-A7V„(^j^). (1,17) 

If the additional species are relic gravitons, then [HH ESI [83] : 

hi £jnGw{i^,To)dlnu = 5.61 X 10-6AiV,(^^^|^), (1.18) 

where z/bbn and z/^ax are given, respectively, by Eqs. (16.611) and (18. 4p . Thus the constraint of Eq. (11.181) 
arises from the simple consideration that new massless particles could eventually increase the expansion 
rate at the epoch of BBN. The extra-relativistic species do not have to be, however, fermionic [82] 
and therefore the bounds on AN^y can be translated into bounds on the energy density of the relic 
gravitons. 

A review of the constraints on AA^,^ can be found in [22] • Depending on the combined data sets 
(i.e. various light elements abundances and different combinations of CMB observations), the standard 
BBN scenario implies that the bounds on AA^j^ range from AN,y < 0.2 to AA^i, < 1. Similar figures, 
depending on the priors of the analysis, have been obtained in a more recent analysis [83|. All the 
relativistic species present inside the Hubble radius at the BBN contribute to the potential increase 
in the expansion rate and this explains why the integral in Eq. (I1.18P must be performed from i/bbn to 
^'max (see also [76j where this point was stressed in the framework of a specific model). The existence 
of the exponential suppression for z/ > z/max (see Figs. Hj) guarantees the convergence of the integral 
also in the case when the integration is performed up to z/ oo. The constraint of Eq. (I1.18P can be 
relaxed in some non-standard nucleosynthesis scenarios [S2], but, in what follows, the validity of Eq. 
(I1.18P will be enforced by adopting AA^,^ ^ 1 which implies, effectively 

hi £jnowii^,To)dlnu < 5.61 X 10'' 47° ^'iVO " ^^'^^^ 

The spectral energy densities illustrated in Figs. [2] and H] are both compatible with the big-bang 
nucleosynthesis bound. Thus the big-bang nucleosyntheis constraint does not forbid a potentially 
detectable signal in the high-frequency branch of the relic graviton spectrum. Potential deviations of 
the thermal history of the plasma must anyway occur before big-bang nucleosynthesis. 
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2 The polarization of relic gravitons and of relic photons 



2.1 Basic notations 

As discussed in the introduction, in the ACDM paradigm the background hne element can be written 

ds^ = y^j^dx^dx" = a^(r)[(ir^ —^^j(x)dx'^dx-'], (2-1) 

where, in the spatially flat case, will coincide with 6ij and the Friedmann-Lemaitre equations can 
be written as 

Sn^ = a^fppi, (2.2) 
2{n^-n')=aHl{pt+pt), (2.3) 
p[ + 3nip,+pt) = 0, (2.4) 

where Ti. = a' /a; the prime denotes a derivation with respect to the conformal time coordinate r. The 
Hubble rate is customarily defined in the synchronous frame where the time coordinate (conventionally 
denoted by t) obeys dt = a{T)dT. Denoting with a dot a derivation with respect to the cosmic time t, 
if = d/a, and, by definition, H = aH. In Eqs. fl2.2l) -( l2^ pt and pt are, respectively, the total energy 
density and the total pressure of the plasma, i.e. 

Pt = Ph + Pc + P-y + Pu + Pa, PM = Ph+Pc, (2.5) 

Pt = ^ + ^-PA- (2.6) 

The total matter fraction of the critical energy density, i.e. ^mo = PMo/Pcrit consists of baryons and 
(i.e. pb) and cold dark matter particles (i.e. pc). In Tabs. [1] and [2] the values of ^mo are given as they 
are inferred within the ACDM scenario. In similar terms = Pa/ Pent denotes the critical fraction of 
dark energy. In what follows, if not otherwise stated, the cosmological parameters will be fixed to the 
best fit of the WMAP-5yr data alone, i.e. 

(fibo, ^co, ^A, ho, n„ e) = (0.0441, 0.214, 0.742, 0.719, 0.963, 0.087). (2.7) 

where Q^o, ^co, denote, respectively, the (present) critical fractions of baryons, CDM particles and 
dark energy; ho fixes the present value of the Hubble rate; Hs, as already mentioned in section [T], is 
the spectral index of curvature perturbations and e is the reionization optical depth. 

At the beginning of the previous section we started by stressing analogies and differences between 
relic gravitons and relic photons. The most important one is that both gravitons and photons carry two 
polarizations. This observation is important for a quantitative understanding of the present endevours 
aimed at measuring the E-mode and the B-mode polarization of the CMB. In the present section 
the description of the polarization of the gravitons will be developed by stressing, when possible, the 
analogy with polarization observables of the electromagnetic field. 

2.2 Linear and circular tensor polarizations 

Recalling that i,j, k, ... are indices defined on the three-dimensional Euclidean sub-manifold, the tensor 
fiuctuations of the geometry are parametrized in terms of the rank-two tensor hij 

5l'^g,, = -a^h^,, /i^ = V./i}=0, (2.8) 
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where Vj is the covariant derivative with respect to ^^j] if = 6ij, Vi = di. In Eq. (12.81) the subscript 
refers to the tensor nature of the fluctuation while the superscript denotes the perturbative order. The 
tensor fluctuation hij{x,T) can be decomposed in terms of the two linear polarizations, i.e. 



h^Ax,T) = J2h^x){x,T)e\^\k), e^e^^ = 26''\ (2.9) 

A 



where A = ©, ® denote the two polarizations and where 

ef^ik) = {d,aj-kb,), (2.10) 

efjCk) = ia,b,+a,bi). (2.11) 

In Eqs. (12.101) and (12. lip , a, b and k = k/\k\ represent a triplet of mutually orthogonal unit vectors, 
i.e. 

a-b = 0, a-k = 0, b ■ k = 0. (2.12) 
If the direction of propagation coincides with the z, the unit vectors k, a and b can be chosen as: 

A; = (0,0,1), a = (1,0,0), 6 = (0, 1, 0). (2.13) 

Using Eq. (|27[3D . Eqs. (I2A0|) and (121T|) become 

e% = 1, S^. = 1. (2.15) 

If k coincides with the radial direction, the unit vectors k, a and b can be chosen, in spherical coordi- 
nates, as: 

k = (sin cos cp, sin ?9 sin cos t?) , (2-16) 

a= (cos-i^cosf/), cos-i^siny), — sin'(9), (2-17) 

b = {sirup, — cosip,0). (2-18) 



Since e\j {k)e\j (k) = 25^^\ it is straightforward to prove that 

hij{x, t)W^{x, t) = 2[h^{x, rf + rf]. (2.19) 

As in the case of electromagnetic waves, it is often desirable to pass from the linear to the circular 
polarizations: 

ef{k) = ^[4fi)+tef^{k)i (2.20) 
ef{k) = ^[e^fi)-t^^{k)]. (2.21) 

Equations (12.201) and (12.211) also imply that e^^^* = e^^^ and el^''* = efj'\ A rotation of a and b in the 
plane orthogonal to k = k/\k\ 



a' = cos a d + sin a 6, (2.22) 
6' = — sin a a + cos a 6, (2.23) 
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implies, using Eqs. f l2.10p and (12. lip . 



e| = efj cos 2a + e% sin 2a, (2.24) 
e% = -e% sin 2a + e% cos 2a, (2.25) 

where the tilde denotes the two transformed (linear) polarizations. Under the transformation given in 
Eqs. (12.221) and (I2.23P the two circular polarizations defined in Eqs. (I2.20p and (I2.2ip transform as 

~ef{k) = e-'^'^e^\k), ~e^\k) = e^^elf (fc). (2.26) 

The transformation properties of the circular polarization under a rotation in the plane orthogonal 
to the direction of propagation are closely analog to the transformation properties, under the same 
rotation, of the polarization of the electromagnetic field. This analogy will now be exploited to 
introduce the E-mode and B-mode polarization. 

Before proceeding with the discussion it is appropriate to recall a very basic aspect of rotations 
which can have, however, some confusing impact of the polarization analysis especially in the case of 
the tensor modes. Consider, for simplicity, a coordinate system characterized by two basis vectors, 
i.e. ix = fcos?? and iy = fsini). If we now perform a clockwise (i.e. right-handed) rotation of 
the axes ix and iy, the rotated basis (e^, e'y) will be given as in Eqs. (I2.22p and (I2.23P by replacing 
(d', b') (e!p, e'y) and (a, b) {cx, iy). Some authors, for different reasons, instead of rotating the 
coordinate system prefer to rotate the polarization vector. If angles are in the right-handed sense for 
the rotation of the axes, they are in the left-handed sense for the rotation of the vectors. 



2.3 Polarization of the CMB radiation field 

The radiation field can be described by the polarization tensor, i.e. 

P^J = E,E*, (2.27) 

where Ei and Ej are the electric components of the radiation field. Assuming, for sake of simplicity, 
that the radiation field propagate along the z axis, then the various entries of pij can be written in a 
matrix form 

"•'^Ie'e' II? )• (2.28) 

where x and y denote the components of the electric field orthogonal to the direction of propagation 
coinciding, in this set-up, with the third Cartesian axis. A full description of the radiation field can 
be achieved by studying the four Stokes parameters |84j conventionally named /, Q, U and V: 



I=\Ex\^ + \Ey\\ V = 2lm[E*Ey], (2.29) 



U = 2Re[E:Ey], Q = \Ex\^ - \Ey\\ (2.30) 



It is immediate, from the definitions (I2.29p and (I2.30p . to write the intensities of the radiation field 
along the different Cartesian axis as a function of the Stokes parameters, i.e. 

\E,f = \E,\'- = (2.31) 

e.e;^^, e;e,^^. (2.32) 
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Equations (I2.3ip and (12.321) can be inserted back into Eq. fl2.28p with the result that 

1/j + Q U-iV 
2[u + iV I-Q 

where 



1 / I + Q U-tV\ - 
PH = o rr , „-T/ T ^ = + (2-33) 



(2.34) 




From Eq. fl2.33p it also follows that 



p,, = ]-{Il + Ua^ + Va2 + Q a^) (2.35) 



2 



where 1 is the 2x2 unit matrix and 



-(!o> 

are the Pauli matrices. Consider now a rotation of an angle a on the plane orthogonal to the direction 
of propagation of the monochromatic wave. It is easy to show that / and V are left invariant while Q 
and U do transform by a rotation of 2a. By indicating with a tilde the transformed Stokes parameters 
the result can be expressed as 

Q = cos 2a Q + sin 2a [/, t7 = - sin 2a Q + cos 2a t/. (2.37) 

Equations fl2.24p - fl2.25p and (12.370 express the fact that the polarization of the graviton and of the 
radiation field do change for a rotation on the plane orthogonal to the direction of propagation of 
the radiation (either gravitational or electromagnetic). It is possible to construct polarization observ- 
ables which are invariant for rotations on the plane orthogonal to the direction of propagation of the 
radiation: because of their properties under parity transformations they are called E-and B-modes. 



2.4 E- and B-modes 

The fluctuations of the geometry induce fluctuations of the Stokes parameters whose spectral properties 
are, ultimately, the aim of CMB polarization experiments. In general terms the fluctuation of each 



Stokes parameter can be written as 

Ai(n,r) = AP(n,r) + Ap)(n,r) + Ar)(n,r), (2.38) 

AQ(n,r) = Ag'Hn,r) + Ag)(n,r) + A[J)(n,r), (2.39) 

Au(n,r) = A(7)(n,r) + A(j')(n,r) + Ai7)(ri,r). (2.40) 



In Eqs. f l2.38p - fl2.40p the superscript reminds that the various fluctuations of the Stokes parameters 
are induced, respectively, by the tensor, scalar and vector modes of the geometry. While some of the 
results of the present section will be generally valid, the focus, in what follows, will be on the tensor 
contribution. Defining the two linear combinations 

A±(n,r) = AQ(n,r)±2Au(n,r), (2.41) 

and denoting with a tilde the transformed quantities, Eq. (12.370 implies that A± (n, r) transform as 

A±(n, r) = e^2i« A±(n, r). (2.42) 
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In more general terms, consider a generic function of n (be it f{n)). Under a rotation of an angle 
a on the plane orthogonal to n, f{n) is said to transform as a function of (integral) spin-weight ±s 
provided 

Rn) = e^^^^ih). (2.43) 

In other words A+(n) and A_(n) transform, respectively, as functions of spin weight +2 and —2. 
The circular polarizations of the gravitons introduced in Eqs. fl2.24p - fl2.25p transform, respectively, as 
functions of spin weight +2 and —2. The brightness perturbations for the intensity of the radiation field 
(i.e. ^i{n)) transform, on the contrary, as quantities of spin weight 0. The fluctuations in the intensity 
of the radiation field, being a spin-0 quantity, can be expanded in ordinary spherical harmonics as 



Ai(n,r) = 



m 



n] 



Im 



a£ = JdhY;^ih)Mh,T) (2.44) 

The spin-s quantity will naturally be expanded in terms of a generalization of the ordinary spherical 
harmonics which are called spin-s spherical harmonics or also spin-weighted spherical harmonics. Ow- 
ing to this observation, A±{n,T) can be expanded in terms of spin-±2 spherical harmonics ±2^£m('^); 
i.e. 

A±(ri, r) = Y^ a±2,em ±2Yim{n). (2.45) 

em 

Given a quantity of spin-weight s it is possible to construct quantities of spin-weight by the 
repeated use of appropriate differential operators which can either raise or lower the spin-weight of a 
given function (see subsection 12.51 for a specific discussion). Consequently, from A+(n, r) and A_(n, r) 
it is possible to construct two fluctuations of spin which can be eventually expanded in ordinary 
spherical harmonics. By demanding that the two fluctuations of spin-weight are eigenstates of parity 
the E- and B- modes are defined as 



AE(n,r) = ^iV7\SF,^(n), iV, = ^ 

AB(n,r) = aP^Yemin), (2.46) 

im 

where ^ 

'^^m — ~2^'^2,^m + '2-2, £m) 5 = '^^(^'^^im — 0^-2, tm)- (2.47) 

From a^^, a£ and a£ the angular power spectra can be defined. In particular the EE, BB, TT 
and TE angular power spectra are given by: 



.(EE) _ ^ \^ ln^^>n^^h r^(BB) _ i \^ I n^^> n^^X ( A9i\ 

'i — 2£-\-l ^ ^™ ^ ~ 2£-\-l ^ "'im/j l^-^oJ 

m=—l m=—l 

,(TT) _ _J_ ^ m^ r^(TE) _ ^ , (T)* (E), 

'I ~ 2£ 4- I ^ ^ ~ 2^ + 1 ™ "'Imh yz.^it) 

m=—£ m=—i 



where (...) denotes the ensemble average. Two further power spectra can be defined and they are: 

^ = 2/ _|_ 1 ^ ^^^^"^ ^im.)^ ^ = 2£ -\- 1 ^ ^irl)- (2.50) 

m=—l m=—l 
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Overall, the existence of linear polarization allows for 6 different power spectra. 

In the minimal version of the ACDM paradigm the adiabatic fluctuations of the scalar curvature 
lead to a polarization which is characterized exactly by the condition a2jm = 0,-2, em, i-e. = 0. 
This observation implies that, in the ACDM scenario, the non- vanishing angular power spectra are 
given by the TT, EE and TE correlations. In the TACDM scenario the TT, EE and TE angular power 
spectra are supplemented by a specific prediction for the B-mode autocorrelation (see section [Tj). 



2.5 Spin-2 spherical harmonics 



Spherical harmonics of higher spin appear in matrix elements calculations in nuclear physics (see e.g. 
the classic treatise of Blatt and Weisskopf [85], and, in a similar perspective the book of Edmonds 
[86]). The comprehensive treatments of Biedenharn and Louk [87] and of Varshalovich et al. [88] can 
also be usefully consulted. 

The spin-s harmonics have been introduced, in their present form, by Newman and Penrose [89] 
and their group theoretical interpretation has been discussed in [90]. The spin-s spherical harmonics 
have been applied to the discussion of CMB polarization induced by relic gravitons in a number of 
papers [91], [921 [93] . They are rather crucial in the formulation of the so-called total angular momentum 
approach. Discussions of the spin-weighted spherical harmonics in a cosmological context can also be 
found in [9ll |95]. The spin weighted spherical harmonics will now be introduced by following the 
spirit of Ref . [9^ which has been also used, with different conventions, in [91] . In subsection 12.61 the 
(equivalent) approach of [92l |93] will be more specifically outlined. 

Th functions ±2Yim{n) appearing in Eq. (12.451) are the spin-2 spherical harmonics [90]. Consider 
the representations of the operator specifying three-dimensional rotations, i.e. R; this problem is 
usually approached within the matrix element, i.e. V^^^,{R) = (j, m'|_R|j, m) where j denotes the 



eigenvalue of and m denotes the eigenvalue of Jz- Now, if we replace m' 
definition of spin-s spherical harmonics in terms of the D^g^(a, /3, 0), i.e. 



-s, J 



we have the 



(2.51) 



where a, P and 7 (set to zero in the above definition) are the Euler angles defined as in [96]. If s = 0, 
V^^^^a, P,0) = ^ (2£ + l)/47rY£m(a, /3) where Yim{c(,P) are the ordinary spherical harmonics. The 
spin-s spherical harmonics can be obtained from the spin-0 spherical harmonics by using repeatedly 
the differential operators: 



1 

(sin^9)' 



•9,5 



sm 6 
sm v 



1 



(sin^9)^ 



(2.52) 
(2.53) 



The notation spelled out in Eqs. fl2.52p and f l2.53p (which is not usual) will be employed to emphasize 



the interpretation of K^' as ladder operators (see ^90j). The operator K^'^ raises the spin weight of a 
function by one unit. Consider, therefore, the ordinary spherical harmonics Yimin) defined as 



Yemin) 



m 



l 2i+ l 
An 



^ {i + m)\' 



fi = cos ^ 



(2.54) 
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2\m/2 



{l-m)\ 



(2.55) 



where P^{^) are the Legendre polynomials and P'^{^) the associated Legendre functions. It is ap- 
propriate to mention here that the factor (—1)™ (i.e. Condon-Shortley phase) can either be included 
in the normalization factor N]/^ or (as it has been done) in the definition of the associated Legendre 
functions appearing in Eq. (12.55^ . When using the recurrence relations of the associated Legendre 
functions the Condon-Shortley phase introduces a sign difference every time m is odd. The conven- 
tions expressed by Eqs. (12.541) and (12.551) will be followed throughout the present discussion and, in 
particular, in section [7] where the correlation functions of the 15-modes and of the B-modes will be 
specifically computed with different techniques. 

According to Eq. (12.431) . Yim transform with s = 0, i.e. they have spin weight 0. By applying once 
K^_^\n) to Yirn{n) we do get a function of spin weight 1, i.e. 



Kf\n)Y,^in) = - 



0^0 + 



smv 



Yemin). 



(2.56) 



We can then apply once mor K^\n) to K^\n)Yim{'n). The result of this simple manipulation will 
be 



dl — cot 19(9^ — 



^ + 

2 - ^ 



2i 



Yfmih) 



(2.57) 



sin^ ' sin -(9 ^ " 

The right hand side of Eq. (12.571) is -^2Yem{'&,f) (up to an overall normalization). Including the 
appropriate normalization factor, ±2Yem{'&,'^), i-e. the spherical harmonics of spin weight s = ±2 are 
given by: 



±2 



Yem{^,ip) 



(^-2)! 



dl - cot M<> - 



± 



2i 



sin sin 



■{d^d^ - cot dd^ 



Yfmin) 



(2.58) 



The spin weights s = ±2 are both needed since the transformation of the polarization involve both 
spin weights (see Eq. (I2.45P ). In fact, since ±2Yem{^,'^) form a complete and orthogonal basis on the 
sphere, i.e. 

(2.59) 
(2.60) 



J2 sY;^{^,<f )sYeU^',^') = - 6{cos^ - cosr) 



A-|-(n, r) can be expanded in terms of ±2^m(''^; V^) as in Eq. fl2.45p . The coefficients off the expansion 
will be given by 



0-2, em 



rfn2y;„(n)A+(n,r), (2.61) 
dn ^2Y;^{h)^.{n,T). (2.62) 
Integrating by parts in Eqs. (12.611) and (12.621) allows for a different form of the expansion coefficients 

0'±2,lm'- 

a2,,m = NiJdnY;^K^2\n)[K^'\n)A+{n,r)], (2.63) 
a-2,en. = Ne fdnY;^,KL-'\h)[Ki-^\n)A4h,T)], (2.64) 



This time, in K^^\h), s = 1 since K'^'' {n)Y£,in{n) is a quantity of spin weight 1. 



'(0) 
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where, as already mentioned, = — 2)!/(£ + 2)!. In Eqs. fl2.63p and fl2.64l) there appear only 
ordinary (i.e. spin-weight 0) spherical harmonics. This occurrence suggests a complementary approach 
to the problem: instead of expanding A±(n, r) in terms of spin-2 spherical harmonics, fluctuations of 
spin- weight can be directly constructed (in real space) from A±(n, r) by repeated application of the 
ladder operators defined in Eqs. ( 12.52^ and (12.53^ . 

The E-mode and B-mode polarization in real space are then, in explicit terms: 

AE(n,r) = -^{K^y{n)[K^^\h)A4n,r)] + K^~'\h)[Ki~'^^ (2.65) 

AB(ri,r) = ^{iri'^(n)[iri')(ri)A+(n,r)] -i^^')(n)[iri"'^(n)A_(fi,r)]}. (2.66) 

The quantities Ae(5^, t) and A^{n,T) can be expanded in terms of ordinary spherical harmonics, as 
already suggested in Eq. fl2.46p : 

Aj,{n,r) =Y.N;\f2yimin), ABin,T) =Y.N;\?J^Y,m{n), (2.67) 

£m Im 

The "electric" and "magnetic" components of polarization are eigenstates of parity and may be defined, 
from a±^ as already mentioned in Eq. fl2.47p : 

'^£m ~ ~-^{0'2,lm + Cl-2/m), O-im — lj{'^2,lm — 0,-2,em) ■ (2.68) 

Under parity the components appearing in Eqs. (12.681) transform 

a?2 - (-1)' «L, a?J (-1)^+1 aL- (2.69) 

Therefore, the E-modes have the same parity of the temperature correlations which have, in turn, 
the same parity of conventional spherical harmonics, i.e. (—1)^. On the contrary, the B-modes have 
(— 1)^+^ parity. The same analysis can be directly performed in real space, i.e. using Eqs. (I2.65P and 
f l2.66p . Denoting the radial direction with n and the tangential directions with ei and 62, the ladder 
operators defined in Eqs. f l2.52p and f l2.53p are consistent with the following choice of ei and 62'- 

ei = (cos-i^cosi/), — cos-i^sini/), — sin-i?), 62 = (— siny9, — cos(y9, 0). (2.70) 

A parity transformation (i.e. a space inversion) implies, in spherical coordinates, that 

r ^ r, ■}} ^ n — {}, ip ^ n + (p. (2-71) 

The transformation (12.710 implies that the two basis vectors defined in Eq. (I2.70p transform as ei Ci 
and 62 — > —62, i.e. while ei does not change 62 flips its sign under space inversion. It follows that 
space-inversion does not flip the sign of AE(n) but it does flip the sign of Ab(5t-), i.e. under the 
transformation (12. 71 p . Ae(5t.) — > AE{n) while AB(n) — — AB(n). 

Using Eqs. (I2.63P and (I2.64p inside Eq. (12.680 a more exphcit expression for and af^ can be 
obtained and it is: 

«S = ^/ dhYUnnK^-\n)[K^'\n)A^in,r)] - KL-'\n)[Ki-'\n)A_ih,T)]}, (2.72) 

«S = -^J t^^>^£™(^)*{i^i'^(^)[i^i'^(^)A+(n,r)] +i^i-^Hn)[Ki"')(n)A4n)]}. (2.73) 

^^As discussed at the end of subsection 12.11 the sign of (p can be flipped. This possibihty is not related to a parity 
transformation and it has to do with the way two-dimensional rotations are introduced. This aspect will also be relevant 
in section [7] for explicit derivations. 
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The contribution of long wavelength gravitons to Eqs. (12.721) and (I2.73P will be discussed in section 
[71 It is often useful to observe that the differential operators appearing in the definition of the spin- 
weighted spherical harmonics (see, e.g. Eq. (12.580 ) can be expressed in terms of the usual differential 
operators arising in the theory of the orbital angular momentum in non-relativistic quantum mechanics 
(see, e. g. [96]). Indeed, recalling that 

L± = e±'^[±9^ + icot§d^l = -id^, (2.74) 

L2 = 1(L+L_+L_L+) + L2, (2.75) 

it can be easily deduced that 

cot^99^ = -^[e-*'^L+ + e^^L_], (2.76) 

d^d^ = '-^e-'^ L+ - e*^ L_ L,] (2.77) 

- coi^d^ - = L^ + dl (2.78) 

sm V 

Equations (I2.76I) - (I2.78I) allow often to express combinations of spin-2 spherical harmonics in terms of 
ordinary (i.e. spin-weight 0) spherical harmonics using the properties of the ladder operators associated 
to the (orbital) angular momentum, i.e. L±: 

L±Ye^ = ,Ji{i + 1) - m{m ± l)F£,„±i, L.F,^ = mF,„, = £(£ + l)y,^, (2.79) 

where L± and obey the well known commutation relations [L±,Lz] = =F-^± and [L+,L_] = 2Lz. 

Looking at Eq. f l2.79p it is tempting draw a parallel between the (orbital) ladder operators and the 
ladder operators raising (or lowering) the spin weight of a given function (see Eqs. ( 12.521) and fl2.53p . 
This problem has been discussed and solved in [90] . It is possible to formulate the parallel in terms of 
a putative 0(4) group. Half of the generators will be connected with the orbital angular momentum 
operators, while the other half will allow to increase (or decrease) the spin weight of a given function. 
The two sets of generators commute. The operators are not directly, though, the ladder operators 
stemming from the second set of generators. This has to do with the fact that in Eq. fl2.5ip the third 
Euler angle (i.e. 7) has been fixed to zero. The are ladder operators defined within a putative 
0(4) group in the case 7 7^ 0. When 7 — the dependence upon 7 drops and we are left with Eqs. 
( 12321) and fl233D . 



2.6 Polarization on the 2-sphere 

In a more geometric perspective, the spin-2 harmonics are introduced by describing the polarization 
tensor on the 2-sphere which represents the microwave sky. In Eq. f l2.34p the tensor P^- describes 
the properties of the radiation field and it is symmetric and trace-free (i.e. Pij = Pji and P- = 0). 
Equation fl2.34p holds in fiat space-time. On the 2-sphere the line element can be written as 

ds^ = d^"^ + sin^ ^^(^^ g^^ = 1, g^^ = sin^ 1^. (2.80) 

The polarization matrix Pij will now be generalized as 

Qin) -Uin)smi} \ 
-Uin) sind -Q(ri) sin^t? J' ^ ' ' 
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satisfying Pab = Pba, and g°'^Pab = 0, where n is a unit vector in the direction (-(9, ip). The sign of the 
off-diagonal entries in Eq. (12.811) is opposite with respect to the one obtained in Eq. (12.341) . This is 
just because we want to match with the conventions adopted, for instance, in [93], [HU [95j. To avoid 
possible confusions, furthermore, the Latin indices a, 6, c, .... run over the two-dimensional space. 

As already mentioned, for scalar functions defined on the 2-sphere, such as the temperature 
anisotropies, the spherical harmonic functions Yim{n) are the complete orthonormal basis. For the 
2x2 tensors defined on the 2-sphere, such as Pab in Eq. (l2.8ip . the complete orthonormal set of tensor 
spherical harmonics can be written as [931 IMl ES] : 

1 



1 1, 



'm)a 



Nf ( VaVbY^tm) 



yC 

-' {em)ab 



Ne 
2 



where V, in this subsection, denotes a covariant derivation on the 2-sphere, e.g 



dadbY(^gmj 



ac 



^ g'"^{-ddgab + dbQdc + dc 



(2.82) 
(2.83) 

(2.84) 
(2.85) 

(2.86) 

:^-2)!/(£ + 2)!, while 

sine \ 
-l/sin^ j' 

is the Levi-Civita symbol on the 2-sphere. The differential operators acting in Eqs. fl2.82p and (12.831) 
are interpreted as a generalized gradient and curl operators, i.e. 



Using Eq. ( ]2.80p into Eq. ( 12.85^ . the Christoffel connections FJ^^ on the 2-sphere are 

1 



tant? ' 



sin cos 'd. 



In Eq. (I2.83P the normalization factor is given bj@ = \Fl 



e b 



(2.87) 



(2.88) 



The explicit form of the various components of Y^^-^^^^ and Y^ can be computed. For instance using 
Eqs. (ISJD, (12:81) and (JSJEI), the explicit components of Y|f^^^^^■. 





_ m 

2 


dlY(^ern) — 




_ 

2 








d^d^Yi^irn) 



d-dY(^lrn) 



siv? § tan 'd 
sin^ ddl Y(^ern) + sin {} cos {}d^Y(^ern)] 



(2.89) 
(2.90) 

The expressions obtained in Eqs. (12.891) . (12.901) and (I2.9ip can be simplified by recalling Eqs. (I2.74p . 
(EZZHD and fl2T9|) . Equations ^M), flCT]) and (ICTD can be simply rewritten as 

= ^[2a2r(,„) +£(£ + l)r(,„)], (2.92) 



^ (em 



(G 



Y, 



^ (em 



(^,y,) = -^[29^F(,„) 



£(£+l)y(,„)]. 



Yi 



(em) 



tant? 



(2.93) 
(2.94) 



i^Notice that Nt differs from defined in Eqs. ([^^ (see also Eqs. (PTM)) and by a factor V2. This difference 

will be ultimately relevant to relate {a'(^, ^Sn) ^tm)- 
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The same exercise can be conducted for the various components of Y(em)ab('^^ '^)- '^^^^ ^{ern)ab 



(C) 

^(£m)ab '^^^ Written in the form of 2 x 2 matrices: 



iV, / W^em)in) X(,„)(n)sin^ 
^ V X^im){n) sm-d -ly(^^)(n)sm i) 



and as 



where 



22771 

W(,„)(n) = [29| + ^(^+l)]y(,^)(n), = (9^-cot^)y(,„)(n). (2.97) 

Sin 7/ 

In terms of the spin-2 harmonics ±'2Y(irn) ip) 



U + 2)! 



(2.98) 



which is Eq. ( I2.58p . Using the orthonormahty of the spherical harmonics Yim{n) it is easy to prove 
the orthonormahty conditions, i.e. 

dnYl^Xbin)Y{Z''in) = 0, (2.99) 

where dn = sin-^d-ddip denotes, as usual, the integration over the solid angle. Since Y^^-^^f^ and Y^^^^^f^ 
form an appropriate orthonormal basis, the polarization can be expanded as 

oo e 

Pabih) = E E [&if2)abin) + afjY[2)Jn)] , (2.100) 

£=2 m=-£ 
(G) (O 

where expansion coefficients a\J and a)^ represent the electric and magnetic type components of the 
polarization, respectively. Note that the sum starts from i = 2, since relic gravitons generate only 
perturbations of multipoles from the quadrupoles up. The expansion coefficients are obviously 

= / dhPab{n)Y{^2;''*{h), afj = J dn P^,{n)Yl^;'\n). (2.101) 
In the notations of |9T] the a}^ and o!"^ can be related to the and already introduced in Eq. 



( ]2.68p . The relation between the two sets of expansion coefficients is simply: 

(G) _ ^ (E) _ t%n ^(C) _ ^ (B) _ ^9 109^ 

"£m — jVf ~ a/2' ™ ~ Ni ~ ^2 ' {^■J-'J^) 

The two approaches to the spin weighted spherical harmonics described in the present section are 
equivalent and can be used interchangeably depending upon the specific problem. 
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3 The action of the relic gravitons 



3.1 Second-order fluctuations of the Einstein-Hilbert action 

By perturbing the Einstein-Hilbert action, to second-order in the amphtude of the tensor fluctuations 
we have, formally, that: 



IGnG 



(2), 



-g)R + 



{3.1] 



where R denotes the background Ricci scalar; S^^^R and S^'^^R denote respectively, the first and second- 
order fluctuations of i? = g^'^R^y. In Eq. (13.11) the possible coupling to the anisotropic stress has been 
neglected. This is customary during the early evolution of the geometry since, in the context of the 
ACDM paradigm, during the early inflationary phase the sources of anisotropic stress can be safely 
ignored unless the number of effective e-folds is close to minimal. Later on the anisotropic stress of the 
fluid plays a role and cannot be neglected at least if we aim at a reasonable quantitative discussion of 
the relic graviton spectrum (see also section [1] and Fig. [3]). By introducing the first-order fluctuations 
of the background geometry 1)^^ we have that 



v-g = v-g 

Recalling now Eqs. and ([M]), Eqs. (Q-lM become 



(3.2) 
(3.3) 

(3.4) 



(3.5) 



The first- and second-order fluctuations of the Christoffel connections are: 



1 



^(2)pi 



2" 



Ml 



(3.6) 



where the prime denotes a derivation with respect to the conformal time coordinate. Using the result 
of Eq. (13. 6p the first- and second- order fluctuations of the Ricci tensor can be written in explicit 
terms: 



V 



00 



(3.7) 
(3.8) 



5r^R 



1~L 



-h'"^[dkdehij - dkdjhu - dkdih 
-Oj^-^ -..-^^ -.-^^..j 2 

jh% + - \h)'h[, - ^h]'K, 



n 

2 



1 
4 



'hi 



(3.9) 
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The Ricci scalar is zero to first order in the tensor fluctuations, i.e. S^^^'R = 0. This is due to the 
traceless nature of these fluctuations. To second-order, however, SI R ^ and its form is: 



a"^ 14 2 4 

+ ^[-ld^[h''\dehl-dkhl-d'hki)] 

- ^[dihidehl - d.hid'ha + dkh'%h1 - d'Kkd'h\ + d'Kkd^hu]Y (3.10) 

Using the resuhs of Eqs. fl3.7p - fl3.10p into Eq. (13. ip the second-order action for the tensor modes 
reads, up total derivatives, 

^gw = 5f ^5 = ^ / d^^^ d^Kjd,h^\ (3.11) 
where, as already mentioned in section [H 



Stt 



?2 = SttG = =j = — 2 , Mp = 1.221 X 10^^ GeV. (3.12) 



3.2 Lagrangian densities 

The action (13.111) can be written in various ways which differ by the addition (or subtraction) of a 
total conformal time derivative. Recalling the standard notations 



J dTLg^{T), Lg^{T) = J d^xCg^{x,T), (3.13) 
the Langrangian density C^^\x, r) can be recast in the form 

/:«(f , r) = ^[{drh)' - rV.hV.h], (3.14) 
where the canonical amplitude h has been introduced 

h = ^ = ^. (3.15) 
By now introducing the canonical amplitude as ah = fi, Eq. fl3.13p can be transformed as 



^) = ^iid^f^y + - - V./iV,/.]. (3.16) 



If a total r derivative is dropped, an equivalent form of the Lagrangian density can be obtained 

4w(^.^) = llidrf^y + (^' + ^')/i' - t^V./iV./x]. (3.17) 

All the three Lagrangian densities of Eqs. fl3.14p . fl3.16p and (13.171) lead to the same Euler-Lagrange 
equations. 
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3.3 Hamiltonian densities 



In Eq. fl3.14p the canonical field is h and the canonical momentum is 11 = ah' . Conversely, in Eq. 
(I3.16P the canonical field is /i and the associated canonical momentum is vf = /i' — Finally, 
according to Eq. (13.17^ the canonical momentum is vr = yu' while the canonical field is always /x. 
The three Lagrangian densities of Eqs. (13.141) . (13.161) and (13.171) will then lead to three corresponding 
Hamiltonians, i.e. 

Hf2{r) = ^ / ^'^[5 + «V^V,W,/i] , (3.18) 

hSI{t) = 1 1 (fx^k^ + 2n^l^ + r^v,/iV,/i], (3.19) 



n" 

2 " ,.2 



(3.20) 



a 

The Hamiltonians of Eqs. ( ]3.18p . ( 13.19^ and ( ]3.20p are related by successive canonical transformations. 
To prove this statement it is enough to show that Eq. (13.191) can be obtained from Eq. (I3.18P by means 
of an appropriate canonical tansformation and that, in turn, Eq. (13.201) can be obtained from Eq. (13.191) 
through another canonical transformation. To pass from the Hamiltonian of Eq. (13.180 to Eq. (13.190 
it is practical to consider a generating functional depending upon the new canonical fields (i.e. /i) and 
upon the old canonical momenta (i.e. H): 

6?(i)^(2)[/i,H] = - l^d'x. (3.21) 

By taking the functional derivative of ^(1)^(2) [/x, H] with respect to H and with respect to /i we get, 
up to a sign, the connection between the new and old pivot variables, namely: 

h = Jlm^ = ^^, ^ = -^%(H)=H (3.22) 
oH a ofx a 

Since the generating functional depends explicitly upon time, the new Hamiltonian will be related to 
the old one through a partial time derivative of the generating functional, i.e. 

dr 

as it can be explicitly verified by using Eqs. (I3.18p . (I3.19P and (I3.2ip into Eq. (I3.23p . A further 
canonical transformation allows to go from Eq. ( 13.190 to (]3.20p : the relevant generating functional is 

6^(2)^(3) [/u, ^] = - / d'x (^fiiT - ^) , (3.24) 

depending upon the old coordinates (i.e. fi) and upon the new momenta (i.e. it). The relations 
between the new and old variables are given by 

7f = %^ = 7r-H/., /. = %^ = /., (3.25) 

stipulating that, in this case, the canonical momentum gets shifted by Tin while the canonical field is 
left invariant. Since the generating functional depends explicitly upon the conformal time coordinate, 
we will simply have that 

HS!ir) = HS!ir) + ^^^ (3.26) 
as it can be explicitly verified by using Eqs. (13.190 . (I3.20p and (13.250 into Eq. (13.261) . 



<Hr) = /7«(r) + ^:^, (3.23) 
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3.4 Evolution equations in different regimes 



From Eq. fl3.1ip the evolution equations of will be given by 

h(' + 2nhf - V^hi = 0. (3.27) 

The canonical field h (see Eq. (13.151) ) will also obey Eq. fl3.27p . The Hamilton equations derived from 
Eq. (13:T8|) read: 

h' = ^, n' = a^V^h, (3.28) 

which has exactly the same content as Eq. (13.271) . In similar terms the Hamilton's equations can be 
derived from Eq. (13.19^ and the result is 

/i' = ^ + Ufi, = -Hn + VV- (3.29) 

Bearing in mind that ah = /i, Eqs. (13.281) and (I3.29P all reduce to Eq. (13.270 since the different 
Hamiltonians are related by canonical transformations. The same conclusion follows by deriving the 
Hamilton's equations using Eq. (13.200 . It is practical, for some applications, to change the time 
parametrization. For instance, in terms of the rescaled time coordinate a we will have that the 
evolution for the canonical amplitude h obeys the simple equation 

— - a^V^h = 0, a\T)dT = da. (3.30) 

Before concluding this section it should be pointed out that Eq. (I3.27P is accurate as long as the 
sources of anisotropic stress are totally absent. This approximation is, strictly speaking, unrealistic. 
Indeed we do know that there are sources of anisotropic stress. Typically, after neutrino decoupling, 
the neutrinos free stream and the effective energy-momentum tensor acquires, to first-order in the 
amplitude of the plasma fiuctuations, an anisotropic stress, i.e. 

6Ti = -6p6i + H| (3.31) 

where H^ is the contribution of the anisotropic stress, satisfying VjH* = and H* = 0. The presence of 
the anisotropic stress clearly affects the evolution of the tensor modes. To obtain the wanted equation 
we perturb the Einstein equations to first-order and we get: 

/if + 2nh{' - V^h{ = -16nGa''U{. (3.32) 

This form of the evolution equation for the tensor modes is the one required to compute the effects 
related to the finite value of the anisotropic stress. 
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4 Quantization of the tensor modes 



There are analogies between the quantum state of rehc gravitons and the quantum treatment of visible 
light. Quantum effects are not crucial to treat first-order interference of the radiation field (i.e. Young 
interferometry) [97]. First-order interference in quantum optics correspond to the calculation of the 
two-point function of the relic gravitons. Quantum effects arise, in optics, from second-order interfer- 
ence, i.e. when computing (and measuring) the interference between the intensities of the radiation 
field. Second-order interference effects are associated with the possibilities of counting photons and 
have been pioneered by Hanbury-Brown and Twiss in the early fifties [97l [98]. Hanbury-Brown-Twiss 
interferometry is based on photon counting statistics. 

Having said that we are not even close (experimentally) to study graviton counting statistics (as 
we do it with the photons), second order interference effects would allow, in principle, to assess the 
coherence properties of relic graviton backgrounds. The quantum state of the relic gravitons can be 
described in terms of a generalized coherent state usually called squeezed state. Squeezed states can 
be described in terms of quadrature operators where one of the modes of the radiation field is always 
broadened by the time evolution, while the other one is squeezed. 



4.1 Heisenberg description 



The quantization of the canonical Hamiltonian of Eq. (13.201) is performed by promoting the normal 
modes of the action to field operators in the Heisenberg description and by imposing (canonical) 
equal-time commutation relations: 



[il{x,T),7T{y,T)] =i6^^\x-y). 
The operator corresponding to the Hamiltonian f l3.20p becomes: 



TT^ A' + {difxY 



In Fourier space the quantum fields fi and n can be expanded as 



(4.1) 



(4.2) 



fi{x, t) 



2(27r)3/2 



2(27r)3/2 



(4.3) 



Demanding the validity of the canonical commutation relations of Eq. (14. ip . the Fourier components 
must obey: 



[Ag(r),4(r)] = i5^'\k-p), [f^l{r)Mr)] = i5^''\k-f>), 
[AfcM,%(r)] = i5^'\k + p), [4(r),4(r)] = i5^'\k+p). 



(4.4) 



Inserting now Eq. (14. 3 p into Eq. (13.201) the Fourier space representation of the quantum Hamiltonian 
can be obtained: 



Hir) 



d^k 



(4.5) 



19 



In order to derive the foUowing equation, the relations fi^ -. = jl^ and -. = should be used . 



it - 
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The evolution of /t and vr is therefore dictated, in the Heisenberg representation, by: 



ifi'=[fi,H], iTx'=[Ti,Hl (4.6) 

where, as usual, units h = 1 are assumed. Using now the mode expansion (14.31) and the Hamiltonian 
in the form (14.51) the evolution for the Fourier components of the operators is 

a" 



A'^ = %' ^'k = -[k--]h^ (4-7) 

implying 



a" 



a 



f^k = 0. (4.8) 



It is not a surprise that the evolution equations of the field operators, in the Heisenberg description, 
reproduces, for /t^ the classical evolution equation derived before. The general solution of the system 
is then 

AfcW = ki^o)fkir) + al.iro)mT), (4.9) 
Mr) = a^{ro)gkir) + alMgUr), (4.10) 



where the mode functions obey: 



f'k — 9k, 9k 



a" 



a 



fk. (4.11) 



If the form of the Hamiltonian is different by a time-dependent canonical transformation, also the 
canonical momenta will differ and, consequently, the relation of g to f may be different. For instance, 
in the case of the Hamiltonian of Eq. (I3.19P we will have, instead, 

f'k = 9k + nfk, g'k = -ngk-efk. (4.12) 

Consider now the canonical commutation relations expressed by Eq. (14. ip . Using Eqs. (14. 3 p together 
with Eqs. (14. 9 p and (I4.10p into Eq. (14.11) . the mode functions have to obey the condition: 

Mr)glir)-g,{r)n{r)=z. (4.13) 

Since, by construction, the Hamiltonians of Eqs. (I3.19P and (I3.20p are related by canonical transfor- 
mations, the mode functions of Eqs. (14. lip and (I4.12p will have both to obey Eq. (14.130 . In different 
terms, the commutation relations between field operators should be preserved by the time evolution 
and this is equivalent to the Wronskian normalization condition of Eq. (I4.13p . 

4.2 Generalized coherent states of relic gravitons 

Consider the Hamiltonian given in Eq. (13.191) in the spatially fiat case: 

^gw(^) = ^J d^xln^ + 'HifiTT + vr/i) + difid'fi], (4.14) 
dropping, for simplicity, the tilde from the momenta. Defining the creation and annihilation operators 



a 



I k / ^ i 
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and recalling that /t^^ = /i^ and tt^^ = vr^, Eq. (14. ISp imply 

+ at , . 

Since ag(r) and o.^-(r) obey [^^(t), at(T)] = 5*^^)(A; — p), inserting Eq. fl4.16p into Eq. fl4.14p . can 
be written as 

4w(r) = -J £k{k[a^al + ata^] + - a^k^'l^l)}. (4.17) 

The evolution of a^(r) and o|^g(r) obeys: 

—A = _^A;ar + Ha) ^, — ^ = ika^ ^ + ^a?- (4- 18) 



The solution of Eq. (14.181) is: 



OfcM = M^)^ki^o) +vl{r)d_^{To), (4.19) 
^U^) = vl{T)d^,{ro)+ulir)dl^^{r,), (4.20) 

where tq is the initial integration time. The unitarity of the time evolution demands that [^^(r)^ — 
|^fc(''")P = 1. A useful parametrization of ^^(t) and Vk{T) is given in terms of a real amplitude and 
two phases as: 

Mfc(r) = e-'^^^^^ cosh rfc(r), Vk^r) = e-*'^fcW-2^^feW sinh rfe(r). (4.21) 

Equation (14.180 determine the evolution equations for Uk{T) and Vk{T). Using then Eq. (14.210 the 
evolution equations for r^, (fk and -dk can be obtained: 

r'^ = ncos2^k, (4.22) 

, (4.23) 

tanh zvk 

^'^ = k + n tanh Tk sin 2ipk. (4.24) 

Note that dk does not appear neither in Eq. (14.221) nor in Eq. (14.231) . It is interesting, at this 
point, to compute the two-point functions connected with the two canonically conjugate operators, 
i.e. fi{x,T) and 7r(x, r). In terms of the creations and annihilations operators defined in Eqs. (14.151) 
and (I4.16p - (l4.17p the canonically conjugate operators can be written as 



r / d^k 

M-^, r) = ^ J -^H-^y-''-' + ^lir)e-''-% (4-25) 



^i^^r) = / rf3fc^^[a,^(r)e-^^^-^-at(r)e-^^i. (4.26) 

There is a slight difference in the normalizations adopted between Eqs. (I4.9p - (l4.10p and Eqs. (I4.25p - 
(I4.26p . This difference is due to the fact that, in Eqs. (l4.9p -( HrT0l) the mode functions fk are normalized, 
asymptotically, in such a way that fk — * l/v^2fc. In Eqs. (I4.15p - (l4.16p the factors \/2k and \Jk/2 have 
been included in the definition creation and annihilation operators. 
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After simple calculations the two-point functions of the field operators and of their related canonical 
momenta becomes : 

poc ^2 sill kv 

(0|/i(f,r)/i(y,r)|0) = — rflnA;-^[|nfc(r)|2 + |t;fc(r)p + M,(rK(r)+<(r)^;fe(r)](4.27) 

roo sin kv 

(0|7r(f,r)7r(f,r)|0) = —^^In/c -^[|w,(r)|2 + |t;,(r)|2 - w,(rK(r) - 4(rK(r)](4.28) 

where r = \x — y\. Again, as already remarked, the non-strandard pre-factors apperaing in the Fourier 
amplitudes of Eqs. (14.270 and (14.28^ are a consequence of the normalizations of Eq. (14.150 . In the 
limit \x — y\ —>■ (and making use of the definitions of Eq. (I4.2ip ). Eqs. (I4.27P and (14.281) lead to 

(0|[/i(f,r)]2|0) = / — din [cosh 2rfc(r) + sinh 2rfc(r) cos 2^fc(r)], (4.29) 

7r(f,r)]2|0) = / — rflnA;[cosh2rfc(r) -sinh2rfc(r) sin2(^fc(r)]. (4.30) 

Equations (14.291) and (14.301) show that the canonical field is broadened while the conjugate momentum 
gets squeezed by keeping constant the product of their respective root mean squares. The latter 
behaviour is evident as soon as the relevant wavelengths are larger than the Hubble radius. To 
demonstrate the two previous statements consider, indeed, Eqs. (14.221) and (I4.23p . Their solution for 
wavelengths larger than the Hubble radius (i.e. to leading order in kr < 1) is: 

rfc(r) ~ In a + C(A;r), (^^(r) ~ C(A;r). (4.31) 

Using Eq. (OB into Eqs. (091) and (lOOD we do get 

Aftk = -^[cosh2rfc(r) + sinh2rfc(r) cos2y;fc(r)] = ^e^''^'^ (4-32) 
ATTfe = i=[cosh2rfc(r) - sinh2r,,(r) cos2^fe(r)] = ^e-^'^'^M. (4.33) 

The rationale for the behaviour exhibited by Eq. (14.310 can be understood also from a slightly different 
perspective. The relation between Ukir) and Vkir) and the tensor mode functions /^(t) and 5'fc(r) is 

Mr) = '''^"^J^'^''^ ^ Mr) = ^[u,{r) - v,{r)]. (4.34) 



where the evolution of /^(r) and gki^r) is obtained by solving Eq. (14.12p . Indeed, by solving Eq. (14.120 
in the limit fc^ < (H'^ + H') we get 

-—- gk{T) = MT)-'Hh. (4.35) 
a"^(r ) 

From Eq. (14.35P the first derivative of fk with respect to r is nothing but 

fk{T)=nfk{T) + ^. (4.36) 
By computing 5'fc(r) from Eq. (14.350 it is clear that, in the limit kr 

fk{r)c^A^{k)a{T), ^,(r)^^, (4.37) 
34 



which has exactly the same physical content of Eqs. (14.321) and (14.331) . When the Universe expands, 
5ffc(r) decreases and that the solution associated with A2{k) becomes progressively subleading. How- 
ever, this observation does not imply that gkir) disappears since the evolution must be unitary. This 
feature of squeezed quantum state suggests the possibility of associating an effective entropy to the 
process of graviton production |104l 11051 1106j . 

In the Schrodinger description quantum state of the relic gravitons is closely related to the squeezed 
states of the radiation field [99] (see also |100l llOlj ). Defining, for practical reasons, 

Xk = fJ'k, Xk = fi_j:, Pk = TT^, Pk = 7r_g, (4.38) 
the wavefunction of the ground state will be, for a given mode k, 

vl/,(x,£) = (xx|S,(rfc,^,)|0) = &-*(^'+*'), a, = e-2^\ (4.39) 

V TT 

where, as previously discussed, ipk —>■ has been assumed. In Eq. (14.391) T,k{rk,Vk) is the so-called 
two-mode squeezing operator |102l I103j which will now be written in its generic form: 

E,(r,, (^,) = exp [rfc(r)(e-2*^'=a,^a_,^ - e^^^'^a^^a^)] (4.40) 

where the creation and destruction operators are the ones computed in tq, i.e. = a^(To) by definition 
of Schrodinger description. The state |0) is annihilated both by and by d_g. These two-modes 
appear simultaneously since gravitons are produced from the vacuum whose total momentum vanish. 
The X and x have been dubbed, in the literature, as superfiuctuant operators (see, e. g., |104[ll05(ll06j ). 

The statistical properties of squeezed states can be addressed by employing a useful analogy with 
quantum optics. The idea is to pretend to resolve single gravitons and to study the statistical properties 
of the second-order interference effects. In the case of the gravitons this problem is addressed by 
defining the (normalized) Glauber correlation function not for the photons (as customarily done) but 
for the gravitons. For simplicity let us consider a single mode of the field. In this case the normalized 
Glauber intensity correlation can be written as 

.<^.(M . AO ^ Mi±|m^, ,4.41) 

the colons denote normal ordering and / denotes the operator corresponding to the intensity of the 
radiation field. The normal ordering is related to the fact that, in the optical domain, most measure- 
ments of the electromagnetic field are based on the absorption of photons via the photoelectric effect 
0. In the case of a single mode of the field Eq. (14.411) can also be written as 

,(^)(t,t)=,(2)(0) = i^!^. (4.42) 

In Eq. (14.421) the normalized two point function is written for coincident spatial points (i.e. At = 0). 
The Hanbury-Brown and Twiss experiment, in some sense, probes directly the properties of g^'^\0). 
In the case of coherent states, g'^{0) = 1: this is the case when the photoelectric counts obey a 



^°Needless to say that there is no analog of photoelectric detection for (single) relic gravitons. In this sense the 
following considerations should be regarded as a conditional predictions based on the analogy between squeezed states 
of photons and squeezed states of gravitons. 
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Poisson statistics. In the case of chaotic hght, the joint detection probabihty greater than that for 
two independent events. This can be verified from Eq. (14.421) by assuming that the state of the 
photons/gravitons is given by a thermal mixture. Equation ( 14.42^ can also be recast in the form 

,(^.(0) - 1 = <(^'f)-<''> . (4.43) 

where ((An)^) = (n^) — (n)^ is the photon number variance and fi = a^a. From Eq. (I4.43p it is also 
customary to define the Mandel's Q-parameter, i.e. 

Q = {n)[g('\0) - 1], (4.44) 

which vanishes exactly in the case of a coherent state. The latter statement can be easily appreciated 
since, for a coherent state, a\a) = a\a). Thus, ((An)^) = |ap = (n). In the case of chaotic (thermal) 
light it turns out that g^^\0) = 2. This result can be easily drived by using the so-called Glauber- 
Sudarshan P-representation of the density matrix, i.e. 

p= f d^aP(a)\a)(a\, Pia) = J-e"!"!'/", (4.45) 
J Tin 

where n is the Bose-Einstein occupation number. In the P-representation we have that 

,(.,0)-l.»^lM!f!. (4.46) 



where = / P{a). By performing the required integrations it is easy to show that g^'^\0) — 

1 = 1, i.e. g^'^\0) = 2. So far it has been shown that while a purely coherent state implies g^'^^ = 1 
(i.e. Poissonian statistic) a thermal state implies that g^'^\0) = 2. In the case of the squeezed states 
it can be shown that ^ 

g^^\0) = 3 + j^, Q = 2{n) + 1, (4.47) 

where {n) = sinh^ is the multiplicity. The coherent state leads to a radiation field with Poissonian 
statistics. Thermal states (as well as squeezed states) have a statistics which is, according to the 
quantum optical terminology, superpoissonian. The latter statement is often dubbed by saying that if 
g^'^\0) > 1 photons are bunched while, in the opposite case (i.e. g^'^\0) < 1) the photons are said to be 
anti-bunched. The quantum optical language is much more effective for a mathematical description of 
the semi-classical limit than the usual considerations related to the limit ^ — >■ 0. Squeezed states are 
genuine quantum states with many particles. They are, in some sense, like coherent states with the 
crucial difference that their statistics is super-Poissonian. The possibility of scrutinizing the statistical 
properties of many-gravitons systems would rely on our ability of resolving single gravitons which is 
not even close to the present technological capabilities. 
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5 Relic graviton backgrounds: observables 



In the literature relic graviton backgrounds are characterized in terms of different quantities and, in 
particular, the most common ones ar 

• the power spectrum VT{k,T); 

• the spectral energy density of the relic gravitons QGw{k,T); 

• the spectral amplitude 5/i(z/, r). 

The three listed variables can be related in different regimes. For instance the power spectrum has 
a simple relation to the spectral energy density when the relevant wavelengths are inside the Hubble 
radius. In section [6] it will be argued that, for numerical applications, the transfer function for the 
spectral energy density is more practical to compute than the transfer function for the power spectrum 
or for the spectral amplitude itself. The power spectrum is actually a strongly oscillating function of 
the conformal time coordinate for wavelengths shorter than the Hubble radius (i.e. kr > 1); in the 
same limit the spectral energy density is asymptotically constant. 

5.1 The tensor power spectrum and the spectral ampUtude 

The two-point function of the tensor modes of the geometry is defined as: 

{hij{x,T)hij{y,T)) = {0\hij{x,T)hij{y,T)\0), (5.1) 

where the state |0) is annihilated by ak,x for A = 0, ©. Recalling Eqs. (12.91) and (14. 3p . fiij = ahij, the 
expansion of hij will then be: 



(5.2) 



where eij{k) has been defined in Eqs. (I2.10I) - (I2.11I) and where 

[H^raix'] = ^^'\k-P)hx'. (5.3) 
In Eq. (15. 2p ^ (and the associated Gk^x) are the (complex) tensor mode functions obeying 

Fkx = Gk,x, Gl^ = -2nGk,x - k^Fk,x. (5.4) 
It is immediate to realize that 

7^ fk,X ^ 1, 



Cfc.A = -{gk,x - yifk.x), (5.5) 



a 



where jk^x and g^^x have been introduced, for a single graviton polarization, in Eqs. (I4.10p - (l4.1ip . 
After computing the expectation value we get: 

(0|/.i(x,r)/.;.(y,r)|0) = (l^l/.Mre-^'-"- (5-6) 



^^It is understood that all the mentioned quantities can be expressed either in terms of the wave-number or in terms 
of the frequency since fc — I'kv. 
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which, thanks to spatial isotropy, can also be written as: 

m{S.r)h){y,rm = J dlnkV^ik^rf-^, (5.7) 
where Vrik, r) is, by definition, the tensor power spectrum: 

VT{k,T)=Ail^\F,{T)\\ F,{t) = (5.8) 



In Eq. (15.71) . V'T{k,T) is the tensor power spectrum which has been implicitly introduced in Eq. 
(11.21) when talking about the customary parametrization of tensor power spectra in the context of the 
ACDM model. It is often useful, for practical applications, to consider /ij,( classical field 

characterized by a Fourier amplitude obeying a specific stochastic average. Then we can write 



hi,{x, t) = ^^-p^ J Skhi,{k, r)e '''■^ h,j{k, t) = h*^{-k, r), (5.9) 
where the Fourier amplitudes obey: 

{h:^{k,T)h,,{p,T)) = —VT{k,T)S^'\k-p). (5.10) 

Following the suggestions of [91], it is useful to introduce, for some applications the stochastic fields 

e,ik) = ^^i^) -^^^{k) ^ ^^^^^ ^ ee(fc)+_ze^(fc)_ 

Equation (]5.10p implies that 

2{el{k)em = 2{e;{k)e2{p)) = '^-^V^{k)5\k - p) , (5.12) 
(e*(^)e2(p)) = (e;(^)ei(p))=0, (5.13) 

where V^lk) denotes the tensor power spectrum and where the factor 2 in front of the averages arises 
as a consequence of the y/2 appearing in Eq. (15. lip . In Eqs. (I5.1ip . (I5.12p and (I5.13P the conformal 
time coordinate is absent. In Eq. (15.101) the conformal time appears explicitly. Indeed, Eqs. (15. lip . 
( 15.12P and ( 15.13^ tacitly assume that /i©(fc, r) = e^{k)Te{k, r) and that h(^{k, r) = e^{k)Te{k, r). This 
factorization is related to the concept of transfer function for the amplitude which will be discussed 
in section [HI The decomposition of Eqs. (15.111) . (I5.12p and (15.131) is useful when all the polarization 
have to be treated simultaneously typically in problems involving long wavelength gravitons (see Eqs. 
(I7.85p - (l7.86p and discussion therein). Furthermore the decomposition of Eqs. (15.111) . (I5.12p and (15.130 
allows to factorize the dependence upon the initial spectrum which is useful for numerical applications. 

It is finally common to characterize the stochastic backgrounds of relic gravitons in terms of the so- 
called spectral amplitude. The definition of the spectral amplitude can be read-off from the definition 
of V-rlk, t). More precisely, taking the limit x ^ y in Eq. (15. 6p we will have that 

{0\h{{x,T)}/j{x,T)\0) = J d\nkVT{k) = 4 J u Sh{iy) dlnu, (5.14) 

where the second equivalence defines the spectral amplitude Shi^) by recalling, once more, that the 
comoving wavenumber is related to the comoving frequency as A; = 2tiv. 
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5.2 Energy- momentum tensors for the relic gravitons 

According to Eq. (13.111) . each polarization of the graviton obeys the action of a minimally coupled 
scalar field. This observation was discussed, in particular, in |109[ IllOj by Ford and Parker (see also, 
e.g. [IH])- Following then Refs. |109l IllOj and within our set of notations the energy- momentum 
tensor of the relic gravitons becomes 



rpu 



44 



(5.15) 



which can be obtained from the action of the gravitons by taking the functional derivative with respect 
to g^^,. By making explicit the sum over the polarizations, Eq. (15.151) becomes: 



rpu 



1 

2i' 



P A 



i.e., even more explicitly 
1 



rpu 



2£? 



24 



(5.16) 



(5.17) 



By now using the same rescalings defined in Eq. (13.151) for /i^ and /i® in terms of the canonical 
amplitude h we do get, from Eq. (15.171) . 



t; = 2 



(5.1^ 



From Eq. (15.181) it is clear that the total energy momentum pseudo-tensor summed over the two 
polarizations of the graviton is twice the energy-momentum tensor of a minimally coupled scalar 
degree of freedom provided the amplitudes of the two polarizations are defined in terms of h as in Eq. 
(I3.15p . The (00) and {ij) components of the energy-momentum pseudo-tensor of Eq. ( 15.18^ are: 



rp[) 



-\{drhf + {dmhf 



(5.19) 

(5.20) 



where dr denotes a derivation with respect to the conformal time coordinate r. Since, by definition, 

^ = pSv, V = -P^Si + l!l, (5.21) 
the energy density and pressure of the relic gravitons can then be written as 



Hgw 



(1) = _ 



gw 



{drhy-^{dmhf 



^7 ^ 

n = — 



d,hd^h--{dmhf6i 



(5.22) 
(5.23) 
(5.24) 



The superscript in the energy density and pressure (i.e. and Pgw) is convenient since different 

prescriptions for assigning the energy- momentum pseudo-tensor will be compared in a moment. 

The other possible definition of the energy-momentum pseudo-tensor stems from the generalization, 
to curved space-times, of the usual fiat space-time approach |111] (see also [1121 1113] ). The nonlinear 
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corrections to the Einstein tensor, will consist, to lowest order, of quadratic combinations of hij that 
can be formally expressed as 

4t; = -6i'^g;, (5.25) 

where the superscript at the right hand side denotes the second-order fluctuation of the corresponding 
quantity while the subscript refers to the tensor nature of the fluctuations. This procedure is essentially 
the one described in |112l 1113] and has been re-explored, in a cosmological context, in [1141 1115] (see 
also |116j ). Recalling the form of the Einstein tensor. 



v{2) 



we obtain 
where Pqo is a total derivative, i.e. 



2 Too = nh',,h''' + liK^h^'' + d,hud%^') + V, 
o 



00) 



^l,ki^ 



From Eqs. (13. 9p and (I3.10p it is also possible to write: 



n 

2 



n 



■^{h)'hik + h\'hkj) 



4 



3 
-( 

4 



where Vij and are further total derivative 



(5.26) 
(5.27) 

(5.28) 



(5.29) 



dk[h {dihij - djhif, - dih 



dj[h (dihik - dkhi - d.hke)] 



-deidkKh' + Kkd'h']. 



V 



R 



2 ' 



i ukt\ 



Therefore, up to total derivatives, the following result holds: 

4 Too = nK,h^' + \{K,h^'' + d.hud'h 

o 

and 



(5.30) 

(5.31) 
(5.32) 



To pass from doubly covariant indices to mixed ones, it is useful to recall that, to second order. 



By looking at the form of the specific terms arising in the previous equation it is clear that Tq = Tq 
and that = 'g^'^T^i. The expressions for and are 



00 



-'o 



r 



(5.34) 
(5.35) 
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where 



T 



-h' h''^' 



aHl 1 6 



(5.36) 



(5.37) 



with = 0. (see also Eqs. fl5.36p and (15.371) ). These expressions coincide with the ones obtained, 
for instance, in |114[ 1115] and are also consistent with the ones of |112[ 1113] . From Eqs. (15.341) and 
(I5.36P the components of the energy and pressure density can be easily obtained since, by definition, 
Pgw = = — T/3. As discussed in Eqs. (I5.16l) -( l5.23l) it is rather useful to derive the explicit 

form of and in terms of the normalized canonical amplitude defined in Eq. (I3.15p . The result 
of this calculation is simply: 



p(2) 

r gw 
p(2) 



3a' 



mhdrh + [dmhd^'h + (drhf 

■[5idrh)''-7{dmh){d"'h)]. 



(5.38) 
(5.39) 



By comparing Eqs. (15. 380 - 05. 390 with Eqs. (15.220 - 05.230 we can remark that the first term appearing 
in Eq. (15.380 is absent from Eq. (I5.22p . Moreover, also Pq-|y and seems to be superficially 
different. As it will be shown in a moment the equivalence of the two approaches is clear as soon as 
the relevant wavelengths are larger than the Hubble radius at a given time. 

Before proceeding with this step, it is relevant to remark that the components of the energy- 
momentum pseudo-tensor given in Eqs. (I5.34p and (I5.35P are not covariantly conserved. However, 
since the Bianchi identity V^^^ = should be valid to all orders, we will also have that: 



0, 



(5.40) 



whose explicit form is 



(5.41) 



Recalling now the components of the energy-momentum pseudo-tensor and the results for the fluctu- 
ations of the Christoffel symbols we have 



+ 3At(pgw + Pgw) ^ Ot 1 - U, 



that can also be written as 



where 



dpg^ 
dr 



37^(pgw + Pgw) = 



(n^ - H') ^, 



"Pgw Pgw + yhia? 



(5.42) 

(5.43) 
(5.44) 
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5.3 The energy density of the reUc gravitons 



By choosing the prescription of Eq. fl5.19p the energy density of the rehc gravitons will be given by: 

1 



(0|To"|0) 



Itt^IO) + 7i:2(o|/i2|0) - H(0|(/i7r + 7r/i)|0) + {0\difid'fi\0)]. 



(5.45) 



The expectation values appearing in Eq. (I5.45P can be computed in different ways. For instance, from 
Eqs. (gSI) and ra- fHTTUD . Eq. (KWf 



d\nk- 



\9k{r)\' + {k' + n')\Mr)\' - nmr)g,{r) + hW.ir)] , (5.46) 



where the functions fkij) and Qkij) are the mode functions obeying: 



f'k — 9k, 



k' 



a 
a 



(5.47) 



If the energy momentum pseudo-tensor has the form derived in Eq. (15.341) , then the energy density 
will have, apparently, a slightly different form: 



1 



dink 



27r2 



\9k{r)\' + {k' - m')\h{r)\' + m[n{r)gk{r) + h{r)gl{r)] . (5.48 



From Eqs. (15.461) and (I5.48P the corresponding critical fractions are 



1 dp^^ 



GW 



Pcrit dink 



1 c?Pgw 
Pcrit dink' 



(5.49) 



The two relevant physical limits are for modes inside the Hubble radius (i.e. kr > 1) and for modes 
outside the Hubble radius (i.e. kr < 1). If fc/Ti > 1, then fk{T) will be, in the first approximation, 
plane waves and gk{T) — ±ikfk{T). In this limit: 



^Gw(^' 



r 



1.5 p2 

" r|/.(r)P 



^Gw(^' 



k'il 



r|/.(r)P 



2P 



(5.50) 
(5.51) 



In the limit A;r > 1 we will have T-P ^ k"^. Thus Eqs. (15.501) and (15.511) coincide (up to corrections 
0{li? /k"^)). In this limit it is also possible to express f2Gw(^, t) solely in terms of the power spectrum. 
Indeed, since 

\fk{T)\' 



44A;3 



V^{k,T) 



we will also have 



n^^^ ik 



r 



^Gw(^' t) 



127^2 
12V? 



VT{k,T) 
VT{k,T) 



1 + 
1 - 



2P 
2F 



(5.52) 

(5.53) 
(5.54) 



In the opposite limit, i.e. when the given wavelengths are smaller than the Hubble radius, the same 
analysis implies ^lQy^{k,T) ^ r^Q^vl^) ''")• short the idea is that, for kr <^ 1, gk = Hfk', indeed, for 
kr <^ 1, Eq. (15.41) imphes 

dr' 



fk{r) = Afea(r) + 5fca(r) J 



9k ^ Hfk 



(5.55) 
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where the second term (going as Bk/a{T)) in Eq. (15.551) is actually negligible for large times. Therefore, 
in the limit fcr <^ 1 it can be easily shown that VL^^{k^ r) 7^ Vt^^{k^ r). It is curious to notice that, up 
to a factor of 2, the first energy density (i.e. Eq. (15.460 ) leads to the same value obtained for modes 
inside the Hubble radius. Therefore, for modes which are inside the Hubble radius: 



Pgw(r) = ^/ rfln^l^lM^ = ^ j d\nkk'VT{k,T). (5.56) 



.3 

Consequently, we can define the critical fraction of relic gravitons at a given time as 



1 fin i-2 

where, in the third equality it has been used that H = aH. Recalling Eq. (I5.14p we can also express 
the critical fraction of relic gravitons in terms of the spectral amplitude. Indeed, according to Eq. 

(EH 

VT{k,T)=AiyVT{i^,T), (5.58) 



and, therefore, 



An' 



n^ii^,T) = ^iy'Shii^,T). (5.59) 

As long as the relevant wavelengths are shorter than the Hubble radius at a given time, different 
prescriptions for assigning the energy-momentum pseudo-tensor lead to the same result (see also the 
discussion in section[H]). In the opposite limit different choices may exhibit quantitative differences. The 
limit of short wavelengths in comparison with the Hubble radius is the relevant one when discussing 
wide-band interferometers. Conversely, the initial conditions for the CMB anisotropics are set when 
the relevant wavelengths are larger than the Hubble radius before equality. 
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6 Relic gravitons from the ACDM scenario 



6.1 Inflationary power spectra 

During the inflationary phase, the tensor power spectrum can be easily computed by solving Eq. (15 ■4p : 



where H^^K^) = Mz) + iY^ (z) is the Hankel function of first kind |117l I118j and where e = —H/H^. 
To pass from Eq. (15. 4p to Eq. (16.11) it is useful to bear in mind the following pair of identities 

+ = a^H\2-e), aH = — -. (6.2) 

r(l -e) 

The second equality in Eq. (16.21) follows (after integration by parts) from the relation between cosmic 
and conformal times, i.e. a(r)(ir = dt. By substituting Eq. (16.11) into Eq. (15. 7p the standard 
expression of the tensor power spectrum can be obtained. In particular, when the relevant modes 
exited the Hubble radius during inflation the power spectrum becomes 

where the small argument limit of the Hankel functions has been taken. In the slow-roll approximation, 
~ V"/3; then Eq. ([63D implies tha1@ 

2 f V \ 128/ y 



The spectral index defined from Eq. (16. 4p is, by definition, 

where the second equality follows from the identities obeyed by the slow-roll parameters |119] . The 
tensor and scalar power spectra are customarily assigned at a reference scale (usually dubbed pivot 
scale) : 

VT^{k)=AT^Y) ' fcp = 0.002 Mpc-\ (6.6) 

where, by definition, is the amplitude of the tensor power spectrum evaluated at the pivot scale 
/Cp. In the case of single-field inflationary models, the power spectrum computed at the pivot scale 
/Cp (i.e. ^t) and the spectral index tt-t can be related. Bearing in mind that the power spectrum of 
curvature perturbations is given, in single field inflationary models, as |119] 

^'^w^K^)..„--^<|)""- '«-^) 

the ratio between the tensor and the scalar power spectra is given by 

VT{k) At _ 
V-jiik) An 

Equation (16.80 implies, recalling Eq. (16.50 . that rx = —Sut. Since there is a direct relation of the 
tensor spectral index to r^, the number of the parameters can be reduced from two to one. In Tabs. 
[1] and [2] the values of rx have been reported as they can be estimated in few different analyses of the 
cosmological data sets. 



Within the present notations, as already established in Eq. (|3.12p . £p = vSttG = 1/Mp = a/SttMp. 
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6.2 Transfer functions for inflationary power spectra 



Equation (16. 6p correctly parametrizes the spectrum only when the relevant wavelengths are larger than 
the Hubble radius before matter-radiation equality (i.e. kr ~ k/Ti. < 1 for r < Tcq). To transfer the 
spectrum inside the Hubble radius the procedure is to integrate numerically Eq. (15. 4p (as well as Eqs. 
(I2.2I) - (I2.4I) ) across the relevant transitions of the background geometry. While the geometry passes 
from inflation to radiation, Eq. (16. 6p implies that the tensor mode function is constant if the wavelength 
associated with the given Fourier mode is larger than the Hubble radius at the corresponding epoch: 



The term proportional to Bk in Eq. (16.91) leads to a decaying mode and Fk{T) is therefore determined, 
for \kT\ <^ 1, by the first term whose squared modulus coincides with the spectrum computed in Eq. 
(16.41) and parametrized as in Eq. (16. 6p . The latter statement is true if the inflationary phase is suddenly 
followed by the radiation dominated epoch since, during radiation, a(r) ~ r. The situation can be 
different if, after inflation, a stiff age takes over. The stiffer case still compatible with causality (i.e. 
Wi = 1) leads approximately to a scale factor a(r) ~ ^JT. In the latter situation the second term in 
Eq. (16. 9p grows logarithmically and cannot be neglected in comparison with the constant contribution. 
The evolution of the background (i.e. Eqs. (I2.2p - (l2.4p ) and of the tensor mode functions (i.e. Eq. 
(15.40 ) should therefore be solved across the radiation matter transition and the usual approach is to 
compute the transfer function for the amplitude |120j i.e. 



Th{k) 



\ (in(r)P)' 



In Eq. (I6.10p . Fkir) denotes the approximate form of the mode function (holding during the matter- 
dominated phase); -Ffc(r) denotes, instead, the solution obtained by fully numerical methods. The 
averages appearing in Eq. (I6.10p refer to the average over the oscillations: as the wavelengths are 
inside the Hubble radius, the solutions are all oscillating!^. The calculation of Th,{k) requires a careful 
matching over the phases between the numerical and the approximate (but analytic) solution. After 
matter-radiation equality, the scale factor is going, approximately, as a(r) ~ and, therefore, the 
(approximate) solution of Eq. (15.40 is given by 

/ N 3ji(A;r) sinkr coskr 

^'=(-) = ^:^^^' ii(^-)= (^-7^- (6-11) 

which is constant for kr < 1. The result of the numerical integration for the amplitude transfer 
function T'^ikjkeq) has been recently revisited in [lOl SI] and it turns out to be 



Th{k / k^q 

where 



l + cA^)+bi(^y. (6.12) 



\ '\k '\k y 

ci = 1.260, bi = 2.683. (6.13) 



^■^The numerical average over the phases introduces some arbitrariness which can be cured by computing directly the 
transfer function for the spectral energy density. 
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The latter result agrees with the findings of |120j who obtain ci = 1.34 and bi = 2.50. The value of 
fceq can be obtained directly from the experimental data (see, for instance, last column of Tabs. [T]and 
|2] implying keq — (9(0.009) Mpc~^). The WMAP 5-yr data combined with the supernova data and 
with the large-scale structure data would give k^q = 0.00999to;ooo27 Mpc~^. A rather good analytic 



estimate for keq is 



fceq = 0.0082879 fMlM^V ^0^^° V'^" Mpc'S (6.14) 
^ V0.1326/ V4.15 X 10-V ^ ' ^ ' 

where the typical value selected for /iq^ro is given by the sum of the photon component (i.e. hlfl^o = 
2.47 X 10~^) and of the neutrino component (i.e. hQ^l-yQ = 1.68 x 10~^): the neutrinos, consistently 
with the ACDM paradigm, are taken to be massless and their (present) kinetic temperature is just a 
factor (4/11)^/^ smaller than the (present) photon temperature. 

Equation (16.141) stems from the observation that the exact solution of Eqs. fl2.2l) - fl2.4p for the 
matter-radiation transition can be given as a(r) = aeq[y'^ + 2y] where y = t/ti. The time-scale 
Ti = Teq{\/2 + 1) is related to the equality time r^q which can be estimated as 

Tea = — --^TT^ = 120.658 ° ^° ° ^° , Mpc. 6.15 

Ho Quo VO. 13267 V4.15 x IQ-y ^ ^ ' 

In the case of the WMAP 5-yr data combined with the supernova and large-scale structure data 
hlVLuQ = 0.1368^g%^37. Consequently, Eqs. fl610l) . fl6lTl) and fl632l) imply that the spectrum of the 
tensor modes is given, at the present time, as 

VT{k, To) = ^^^^T',{k/keq)VT{k). (6.16) 

Within the standard approach, Eq. (I6.16P is customarily connected to the spectral energy density of 
the relic gravitons. In ^401 E] it has been observed that it is simpler and more accurate to compute 
directly the transfer function for the spectral energy density. In the following subsection this procedure 
will be illustrated in two different cases. 



6.3 Transfer function for the spectral energy 

Instead of computing the transfer function for the power spectrum it is more direct (and more accu- 
rate) to compute directly the spectral energy density and its related transfer function. As previously 
discussed there can be ambiguities in assigning the energy density because of possible different forms of 
the energy-momentum pseudo-tensor. In particular in Eqs. fl5.46p and fl5.48p two different expressions 
have been proposed on the basis of slightly different physical considerations. The result of the numer- 
ical calculation are reported in Fig. Oin terms of A^^\k,,x) (derived from Eq. (15.461) ) and in terms 
of the transfer function of the spectral energy density (denoted by Tp(K)). The quantities A^^\k,x) 
(and, analogously, A^p\k,,x)) are defined as 

^^Hk,T) = {k(r)P + (A:^ + 7^^)|^(r)r-7^[/:(r)^,(r) + ^(r)^*(r)]}, (6.17) 
^?Hk,r) = {k(r)r + (A:2-77^^)|^(r)r + 37^[/,*(r)^,(r) + ^(r)^*(r)]}. (6.18) 
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Figure 5: The functions given in Eqs. (16.171) and fl6.18p are numerically computed (plot at the left) 
for different values of n and in the case of the radiation-matter transition. In the plot at the right the 
transfer function for the spectral energy is illustrated. 

While Eq. f lGTTD follows from Eq. flCT]) . Eq. fl6J8|) follows from Eq. fICTD : Eqs. fl6Tri) - fl6A8|) are 
related to the spectral energy densities in critical units, i.e. 

^%{k,r) = J' Af\k,T), nSUk,r) = ^J' A?\k,r). (6.19) 

As a function oi x = kr and k = k/keq, both A^^\k, x) and A^^\k, x) go to a constant value when the 
relevant modes are evaluated deep inside the Hubble radius. This occurrence allows to introduce the 
energy transfer function which is defined as: 

hm A(i'2)(«:,x) = r2(«:)A(i'2)(«:,xO, x, « 1. (6.20) 
The specific form of the energy- momentum tensor is immaterial for the determination of T'^{k): 



K = k/k =10"^ K = k/k =10""^ 




log X log X 



Figure 6: The different definitions of energy-momentum pseudo-tensor (i.e. Eqs. (16.171) and (16.181) ) 
are compared in the determination of the asymptotic value of the energy transfer function. 

different forms of the energy-momentum tensor of the relic gravitons will lead to the same result. This 
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aspect can be appreciated by looking at Fig. [6] where A^^'^)(A;, r) has been reported for k = 10^^ (plot 
at the left) and for k = 10"'^ (plot at the right). The dashed and the dot-dashed curves (in both plots) 
correspond, respectively, to A^^\k,x) and to A^^\k,x). The full line, in both plots, corresponds to 
the combination 



(6.21) 



where c±{k) are the mixing coefficients which parametrize the solution for the tensor mode functions 
when the relevant wavelengths are, asymptotically, inside the Hubble radius, i.e. 



fkir) 



c+{k)e~'''^ - c_{k)e 



(6.22) 



In Eq. (16.221) /^(r) and (7fe(r) are the solutions to leading order in the limit /cr ^ 1. From Eq. (I6.22p . 
c±{k) are given by 



c+{k) 



[kfk{r) + iUr)i c4k) 



-ikr 



lkfki^)-Wki^)], 



(6.23) 



Using Eqs. flQ^jl - flOSj) . Eqs. fleTTjl - flerTSjl can be directly calculated in the limit x = A;r > 1 with 
the result that 

A«(«:,Xf) = Af\K,xi) = Ki\c4K)\' + \c4k)\') + o(-), (6.24) 

which proofs that the oscillating contributions are suppressed as x^'^ for Xf ^ 1. To get to the 
results illustrated in Figs. O and E] the evolution equations of the mode functions have been integrated 
by setting initial conditions deep outside the Hubble radius (i.e. x = fcr ^ 1), by following the 
corresponding quantities through the Hubble crossing (i.e. x ~ 1) and then, finally, deep inside the 
Hubble radius (i.e. x ^ 1). The integration of the mode functions is most easily performed in terms 
of f^ix) = kfk{T). Using f^ix), Ap(K, x) = A^^)(k, x) can be written as 



Ap(K,x) 



1 + 



din a 

dx 



+ 



dx 



d\na 

dx 



dx 



dx 



(6.25) 



In the case of the conventional ACDM scenario, the Universe is suddenly dominated by radiation as 
soon as infiation ends. Equation (15.41) implies that Fk is constant when kr <^1. The initial conditions 
are fixed by requiring that, at the initial time of the numerical integration, 



f^{xi) = kfkin) = ka{T)Fk 



Xi = kn <^1. 



(6.26) 



In Figs. [5] and [5|, Xi = 10 ^ even if, for practical reasons, the scale on the horizontal axis has been 
narrowed. Bearing in mind Eq. (I6.15p . we can also write, for x\ <^ 1, 



fniXi) 



2a, 



eq 



Fk Xii 



dx 



2a, 



(6.27) 



To avoid unnecessary complications, the initial condition of the integrations illustrated in Figs. [5] and 
in] have been set as fn{xi) = Xi, i.e. the initial spectrum has been rescaled. The transfer function, by 
definition, must always depend only on the dynamics of the transition and not upon the features (e.g 
slope, amplitude) of the initial power spectrum. 
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In the plot at the right of Fig. [5], the fit to the energy transfer function is reported with the full 
(thin) line on top of the diamonds defining the numerical points. The analytical form of the fit can 
then be written as: 



Tp{k/K^) = y 1 + C2 (^-^) + 62 [-f) , C2 = 0.5238, 62 = 0.3537. (6.28) 

Equation (16.281) permits the accurate evaluation of the spectral energy density of relic gravitons, for 
instance, in the minimal version of the ACDM paradigm. 

Yet another relevant physical situation for the present considerations is the one where the back- 
ground geometry, after inflation, transits from a stiff epoch to the ordinary radiation-dominated epoch. 
In the primeval plasma, stiff phases can arise for various reasons. Zeldovich [122] (see also |123] ) sug- 
gested this possibility in connection with the entropy problem. In [THl [TBI EH EH] it has been suggested 
that the stiff phase could take place after the inflationary phase with the main purpose of identify- 
ing a potential source of high-frequency gravitons. This possibility was also prompted by a possible 
post-inflationary dominance of a quintessence field. 

The simplest consideration leading to the possibility of a post-inflationary phase stiffer than radia- 
tion is connected with our extreme ignorance of the thermal history of the Universe after inflation. In 
a model-independent approach, it is plausible to think that the onset of the radiation-dominance could 
be delayed. This may happen, in concrete models, for various reasons. One possibility could be that 
the inflaton field does not decay but rather changes its dynamical nature by acting as quintessence 
field [78] (see also |126j ). In this kind of situations we are the geometry passes from a stiff phase where 

w,{t) = ^>i, (6.29) 
Pt <J 

2/N K w[ l dln(wt + l) ^ 1 

= ;^ = "^- 37^(^, + l) ="^-3 dlna ^3' ('•'^^ 

to a radiation-dominated phase where c^t = l/"\/3. Note that, according to Eqs. f l6.29p and f l6.30p . 



= Wt iff the (total) barotropic index is constant in time. In the limiting case Wt = 1 = c^^ and 
the speed of sound coincides with the speed of light. As argued in [124] . barotropic indices > 1 
would not be compatible with causality (see, however, [125] ). The presence of a suitable stiff phase 
has been also discussed recently as an effective way of suppressing entropic fluctuations [127] which 
are observationally constrained by the WMAP 5-yr data. 

As in the case of the matter-radiation transition the transfer function only depends upon k which 
is defined, this time, as k = k/k^, where k^. = and parametrizes the transition time. A simple 
analytical form of the transition regime is given by 



a{y) = asJy^ + 2y, y = —, r, = -^J — , (6.31) 



s 




where, by definition, psi = Psi'Ti) and pRi = PR(ri). Equation (16.311) is a solution of Eqs. (I2.2p - (I2.4I) 
when the radiation is present together with a stiff component which has, in the case of Eq. (I6.3ip a 
sound speed which equals the speed of light. In the limit y ^ the scale factor expands as a{y) = \fly 
while, in the opposite limit, a{y) ~ y. In Fig. [7| (plot at the right) Ap(/t,a;) is illustrated for different 
values of k. We shall not dwell here (again) about the possible different forms of the energy momentum 
pseudo-tensor: provided the energy density is evaluated deep inside the Hubble radius the different 
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approaches to the energy density of the rehc gravitons give the same result. The transfer function 
for the spectral energy density is numerically illustrated always in Fig. [7] (plot at the left). The 
semi-analytical form of the transfer function becomes, this time, 

Tl{k/K) = 1.0 + 0.204 (^)'^^ - 0.980 (^)'^' + 3-389(^^) - 0.067 (^^) In^ {k/h), (6.32) 

where = t^^. The value of k^ can be computed in an explicit models but it can also be left as a free 
parameter. Taking into account that the energy density of the inflaton will be exactly psi — H^Mp, 
the value of kg (as well as the duration of the stiff phase) will be determined, grossly speaking, by 
Hi/Mp. In Fig. ([7]) (plot at the left) the full line superimposed to the numerical points (illustrated 




log K log X 

Figure 7: The transition between the stiff phase and the radiation phase is illustrated. The energy 
transfer function increases with the frequency while the opposite is true for the radiation-matter 
transition (see Fig. [5] where the analog results have been presented for the matter-radiation transition). 

by boxes) is the fit of Eq. (16.321) . 



6.4 Analytic results for the mixing coefficients 

The analytic results for the mixing coefficients are rather useful to obtain the final expression of 
the various transfer functions. Indeed, defining as k^, the typical wavenumber of the transition (e.g. 
fc* = keq in the case of the radiation-mattter transition), the slope of the transfer function of the 
spectral energy density can analytically obtained in the limit 1. This observation helps when we 
have, for instance, to fit the numerical data points with an analytical expression which will however 
reproduce the data not only for k > 1 (as Figs. [5] and [7] clearly show). 

For illustration of the method it is practical to consider the transition from a generic accelerated 
phase to a decelerated stage of expansion. In this situation, by naming the transition point — ri, the 
continuous and differentiable form of the scale factors can be written as: 



a; r 



aJT] 



T 



T < -Ti, 









a 




(-- 




+ 1 




r > 













-n, 



(6.33) 
(6.34) 
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In the context of quintessential inflation [7S] (see also [751 EZ]) PRi — Hf [128J 
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where the scale factors are continuous and differentiable at the transition point which has been gener- 
ically indicated as Ti. The pump fields of the tensor mode functions turn out to be: 



aia — r 



r + 



n2 ■ 



(6.35) 



The mode functions can then be written as: 



Mr) 
Mr) 



'2k 



X] 



T < -Ti, 



X = kr, 



[Mc+{k)H^^\y)+M*c.{k)Hi:>{y)] 



where 



X = kr, y = kr + kri ^1 + — ^ 



r > -Ti, 



(6.36) 



^ i(u+l/2)7r/2 

2 



M 



^g-i(A+l/2)7r/2_ 



The continuity of the tensor mode functions at the transition point 

M-ri) = /s(-n), 9i{-ri) = 9s{-ri), 
implies that the mixing coefficients are given by: 



(6.37) 



(6.38) 



c+ik) 



c-{k) 



m2\ 



a 



{2u+l)]Hi'\x,)H^^\y, 



- 2ax,[H^\y,)H%{x,) + H^^\xr)H^lM)]} 



r(l) 



(1) 



{[/3(2A + 1) + a(2v + \)\E^\x^)H^^\y^ 



2ax,[H^^\y,)Hl'l{x,) + hI'\x,)H^^IM]} 



r(l) 



(2) 



(6.39) 



where, according to the notations previously estabhshed, yi = y{—Ti) = {a/P)xi. The case a = [3 = 1 
corresponds to a transition from the infiationary phase to a radiation-dominated phase. In this case 
we do know which are the mixing coefficients. The previous expressions give us: 



c4fc) 



■,2ixi 



2x1 



2 ' 



c+{k) 



1- — + — 



(6.40) 



which clearly agree with previous results [HI HH [50] . In the case of Eq. (16.401) |c+(A;)p — |c_(A;)p = 1 
and fc^|c_(fc)p is exactly scale-invariant. Another interesting situation is the one of the transition from 
inflation to stiff, i.e. /? = 1, a = 1/2, yi = xi/2 which leads to a logarithmic enhancement at small 
wavenumbers [75| [76] . In this situation the mixing coefficients can be written as: 



c.{k) 



c+ik) 



Vf ^^^^"'{^^""'^'[^1 + - 12]//i'^(xi/2) 



2 4xT 

+ (2 + Xi)[iXiH[^\xi/2) - 3iH^^\xi/2)] 



(^ + 1) 
2 4^ 



e''''\xiH!,'\xi/2) + i{i + xi)H["{xi/2) 



(6.41) 
(6.42) 
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The above result can be expanded in for xi ^ 1 and the result is: 

c+{k) = -0-398(1 - ^ ^ ^(Q^^3^ ^ 0.338 - 0.149 (1 - i) Inxi] + 0{xf^), (6.43) 
(7.031 - 1.723 i) - 16.68 (1 + i) Inxi 



c_(/c) -- 



+ y/xl [(-0.621 + 0.265 i) + 0.282 (1 + i) In xi] + C(xi^/^). (6.44) 

The logarithms arising in Eqs. (16.431) and fl6.44p explain why, in Eq. fl6.32p . the transfer function of 
the spectral energy density contains logarithms. In spite of the fact that semi-analytical estimates 
can pin down the slope of the transfer functions in different intervals, they are insufficient for a 
faithful account of more realistic situations where the slow-roll corrections are relevant and when other 
dissipative effects (such as neutrino fee streaming) are taken into account. 

6.5 Exponential damping of the mixing coefficients 

Consider the case of the ACDM paradigm where the inflationary epoch is almost suddenly followed 
by the radiation-dominated phase. By denoting the transition time as T[, it is plausible to think that 
all the modes of the field such that k > UiHi ~ r^^ are exponentially suppressed |129l 1130] . For the 
modes kri > 1, the pumping action of the gravitational field is practically absent. There will be a 
given k, be it fcmax; for which k ^ Tj^^. The latter wavenumber is, in practice, the maximal k to be 
amplified and it can be estimated as : 

/ /)2q X 1/4 

= 7.5959 xl02^(7re^^)V^(^^|°^^) Mpc-\ (6.45) 

/ 7i2o X 1/4 

= 1.1745 x 10" (^^-4^)^^' ( 4 15^10-0 ^^'^^^ 



Equations (16.451) and f l6.46p are derived by assuming that, right after infiation, the radiation-dominated 
phase takes over. Furthermore, recalling the slow-roll dynamics, SH^Mp ~ V and V oc MpeVn. In 
Eqs. (I6.45P and (I6.46P Ati denotes, as already established, the amplitude of the curvature power 
spectrum evaluated at the pivot scale. 

By taking as typical values of the curvature perturbations at the pivot scale the one endorsed, for 
instance, by the WMAP 5-yr data alone we will have we will have that Eqs. (I6.45P and (I6.46P can be 
written as: 



fcmax = 3.5661 X lO^M ^ ' J Mpc-i, (6.47) 

VO.Ol/ V2.41 X 10-9/ V4.15xl0-V ^ ' ^ ' 

where the typical values of the slow-roll parameter have been derived by taking into account that, in 
the absence of running of the tensor spectral index, rx = 16e; since, according to the WMAP 5-yr data 
alone, < 0.43, e < 0.01. For phenomenological purposes it can be also interesting to know what 
kind of exponential suppression we can expect. From the analysis of various transitions it emerges 
that the mixing coefficients for k > fcmax (or u > z/max) will satisfy 

|c+(A;)|2- |c_|2 = 1, |c+(A;)|2 + |c_(A;)|2 = e-'^^ + l. (6.48) 

From Eq. (I6.48P we can easily argue that, for k > /cmax, |c+(/i;)| — > 1 and |c„(A;)| 2^^/^ exp [—Pk/kmax]- 
The point is then to estimate the value of f3 which depends on the nature of the transition regime. 
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Figure 8: The time evolution of the mixing coefficients is reported at the left (on the horizontal axis 
the scale is linear). The exponential decay of the mixing coefficients is illustrated in the plot at the 
right in double logarithmic axes. 

Typically, however, (3 > 2 for sufficiently smooth transitions. To justify this statement it is interesting 
to consider the following toy model where the scale factor interpolates between a quasi-de Sitter phase 
and a radiation-dominated phase: 

a{r) = a-\T^ ^T^ + T^]. (6.49) 

For r — > —oo (i.e. r ^ — Ti) , a(r) ^ —a-Jr and the quasi de-Sitter dynamics is recovered. In the 
opposite limit (i. e. r ^ +Ti), a(r) aiT and the radiation dominance is recovered. In Fig. [6] (plot 
at the left) the exponential damping of the mixing coefficients is numerically illustrated. The curve 
at the top (full line) illustrates the The 2 and k = 3 are barely distinguishable 

at the bottom of the plot. Notice, always in the right plot, the rather narrow range of times which 
are reported in a linear scale. In the plot at the right the asymptotic values of the mixing coefficients 
are reported for different values of k = A;//cmax- By ffiting the numerical data with with an equation 
of the form given in Eq. fl6.48p . the value of (3 = 6.33. Different examples can be presented on the 
same line of the one discussed in Fig. [61 While it is pretty clear that the decay is indeed exponential, 
the value of /3 may well vary. This can be summarized, for instance, in a rescaling of fcmax, i-e. by 
positing, for instance that fcmax ^max/ P- Thus, the dynamics of the transition can slightly shift 
the numerical value of the upper cut-off by a numerical factor which depends upon the width of the 
transition regime. 



6.6 Nearly scale- invariant spectra 

By using the transfer function for the tensor amplitude, the spectral energy density for frequencies 
ly ^ Ucq can be simply given by: 



hlflGw{i^,ro) = A/'hrTl — I e~^^'^ , (6.50) 



. 7.992xlO-^(ggg)-(JgB^)( ^^^^^^^^;^,^,^J -', ,.5i; 
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where dA^z^) is the (comoving) angular diameter distance to decouphng. In Eqs. fl6.50p - fl6.5ip (as 
well as in the program used for the numerical calculations) there are two complementary options. The 
first one is to use the angular diameter distance to decoupling which is directly inferred from the CMB 
data. For instance, in the case of the 5-yr WMAP data alone, dA{z^) = 14115 Mpcl}^^ This is the 
strategy also adopted in other studies |132] (in connection, obviously, with earlier releases of WMAP 
data). At the same time it is also possible to take the best fit value of the total matter fraction (i.e. 
^Mo = 0.258 for the case of the WMAP 5-yr data alone) and compute the comoving angular diameter 
distance according to the well know expression for spatially fiat Universes: 

1 f^* dz 3 375 

dA{z,) = — / , = — — = 14072 Mpc, z, = 1090. (6.52) 

The latter strategy has been used, for instance, in |131j . The two strategies are compatible and, 
moreover, this explains why, in Eq. fl6.5ip the dependence upon ^mo does not cancel. In Eq. (16.501) 
riT denotes, as usual, the tensor spectral index which can be also written as 



nT = -2e + — ln(A;/fcp), «t = — 
Z o 



[ns - 1) + ? 



(6.53) 



If ax = 0, the standard case is recovered. It should be finally mentioned that, in the limit u ^ u^q, 
the oscillating terms have to be appropriately averaged: this is done by setting the terms going as 
cos^ (27rz/ro) to 1/2. The latter procedure has been employed, for instance, in the analyses of jl32l [38]. 
In Fig. [9]the transfer function for the spectral energy density (introduced in Eqs. (16.201) and fl6.28p ) has 
been consistently employed to estimate the spectral energy density itself and the spectral amplitude. 
According to Eq. (15.140 the spectral energy density can be written as 

fiGw(^^,r) = ^z/35,(z/,r), (6.54) 

where we used that k = 2tti' and that Ha = Ti. By making explicit the numerical factors in Eq. 
fl6.54p . Sh{i', To) can be expressed in terms of the spectral energy (in critical units) 

/lOO Hz\ 

5,,(z/,ro) = 7.981 X 10"^=^ (^^j /i^(^gw(^^, tq) Hz-\ (6.55) 

where, we recall, Hq = 3.24078 x 10^^^ ho Hz. In Fig. [9] the spectral energy of the relic gravitons as 
well as Sh{i^, To) are reported for different values of r^. The spectra of Fig. [9]have been obtained from 
the direct integration of the mode functions but can be parametrized, according to Eq. fl6.28p as 



hlnGw{i^,ro) = MpTl{v/v,q)rJ-\'e~^^^ (6.56) 
AC = 4.165xl0-"f ,,5""° V (6.57) 



^-15 ( ^O^RO 

x4.15 X 10-5, 

By comparing Eqs. fl6.50p - fl6.5ip to Eqs. fl6.56p - fl6.57p . the amplitude for z/ 3> t'eq differs, roughly, by 
a factor 2. This coincidence is not surprising since Eqs. f l6.50p -f l6.5ip have been obtained by averaging 
over the oscillations (i.e. by replacing cosine squared with 1/2) and by imposing that \gk\ = k\fk\. 
These manipulations are certainly less accurate than the procedure used to derive the transfer function 
for the spectral energy density. 

So far the evolution of the tensor modes has been treated in the absence of anisotropic stress. This 
approximation is, strictly speaking, unrealistic. Indeed we do know that there are sources of anisotropic 
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Figure 9: The spectral energy density of the rehc gravitons (plot at the left) and the related ^^(z/, tq) 
(plot at the right) for different values of rx and for the same set of fiducial parameters illustrated in 
Fig. 



stress. After neutrino decoupling, the neutrinos free stream and the effective energy-momentum tensor 
acquires, to first-order in the amplitude of the plasma fiuctuations, an anisotropic stress, i.e. 

6Ti = -6p6i+Ui, d,U)=Ul = 0. (6.58) 

The presence of the anisotropic stress clearly affects the evolution of the tensor modes. To obtain the 
wanted equation we perturb the Einstein equations to first-order and we get: 

/if + 2nh{' - V^hi = -levrGa^n^ (6.59) 

Equation fl6.59p reduces to an integro-differential equation which has been analyzed in (see also 
[371 [381 [39] ) . The overall effect of coUisionless particles is a reduction of the spectral energy density of 
the relic gravitons. Assuming that the only coUisionless species in the thermal history of the Universe 
are the neutrinos, the amount of suppression can be parametrized by the function 

J^{R^) = 1 - 0.539i?^ + O.UARl (6.60) 

where Ri, is the fraction of neutrinos in the radiation plasma. In the case R^ = (i.e. in the absence 
of coUisionless patrticles) there is no suppression. If, on the contrary, R^ ^ the suppression can even 
reach one order of magnitude. In the case Ni^ = 3, R^ = 0.405 and the suppression of the spectral 
energy density is proportional to J-"^ (0.405) = 0.645. This suppression will be effective for relatively 
small frequencies which are larger than z/gq and smaller than the frequency corresponding to the Hubble 
radius at the time of big-bang nucleosynthesis, i.e. 

z^bbn = 2.252 X 10-'' ( ( ^\ ( -J^h^\ Hz. (6.61) 



VlO.75/ V MeVy V4.15 x 10- 

The effect of neutrino free-streaming (as well as all the other late time effects addressed in this subsec- 
tion) has been taken into account in Fig. [2] (see section [1]) but it is absent in Fig. [91 The second effect 
which has been taken into account in Fig. [2] is the damping associated with the (present) dominance 
of the dark energy component. Indeed the redshift of A-dominance is simply defined by 

,,,,^(2l)^l^f, (6.62) 
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Consider now the mode which will be denoted as /ea, i.e. the mode reentering the Hubble radius at ta. 
By definition /ca = H\a\ must hold. But for r > ta is constant, i.e. H\ = Hq where Hq is the present 
value of the Hubble rate. Using now Eq. (I6.62p . it can be easily shown that /ca = (f^Mo/^A)^^'^^H 
where /ch = cloHq. The frequency interval between and i^a is rather tiny. Indeed, it turns out that 

= 1.686 X 10- (Jl^)(I}m.y"(J}±.y" Mpc-. (6.63) 



0.719/ V0.258/ V0.742 
= 2.607 X 10-- ( ^) (^) Hz. (6.64) 



0.719/ V0.258/ V0.742 

For the same choice of parameters of Eq. fl6.64p . z/h = 3.708 x 10~^^ Hz which is not so different 
than i/A = 2.607 x 10~^^ Hz. It would therefore seem that this branch of the spectrum could be easily 
neglected. However, it turns out that the adiabatic damping of the mode function across t\ reduces the 
amplitude of the spectral energy density by a factor (^mo/^a)^- For the typical choice of parameters 
of Eqs. (16.631) and (16.641) we have that the suppression is of the order of 0.12. Again this is a tiny 
number which is, anyway, comparable with the suppression due, for instance, to the neutrino free 
streaming. This class of effects has been repeatedly in a number of recent papers |137j (see also |145j ). 
The essence of the effect is captured by the following observation. Consider a mode k which reenters 
before t\. The present value of the amplitude F].{t) = fk{T)/a{T) will be adiabatically suppressed 
since, as repeatedly stressed, in this regime /fe(r) will simply be plane waves. Consequently, defining 
as Fk^ the amplitude at fc* = iJ^,a* when the given mode crosses the Hubble radius, we will also have 
that 

F.(r.)-m r-^]F...ar(^]h., (6.65) 



CJA / mat \ / A VA;h/ V ^A 

where the subscripts (in the first equality) denote the time range over which the corresponding redshift 
is computed, i.e. either matter-dominated or A-dominated stages. The second equality follows from 
the first one by appreciating that a(fc*) ~ ~ k'"^ and by using Eq. (16.621) . Equation (16.651) implies, 
immediately, that the spectral energy density of relic gravitons is corrected in two different fashions. 
For z/ < z/h the frequency dependence will be different and will be proportional to flGy^{h',TQ) oc 
('^/'^h)"'^~^(^mo/^a)^- Vice versa, in the range z/ > z/h the frequency dependence will be exactly 
the one already computed but, overall, the amplitude will be smaller by a factor {Qmo/^aY ■ The 
modification of the frequency dependence is only effective between 0.36 aHz and 0.26 aHz: this effect 
is therefore unimportant and customarily ignored (see, for instance, |132[ 1137] ) for phenomenological 
purposes. On the other hand, the overall suppression going as (Omo/^a)^ must be taken properly into 
account on the same footing of other sources of suppression of the spectral energy density. There is, in 
principle, a third effect which may arise and it has to do with the variation of the effective number of 
relativistic species. Recall, indeed, that the total energy and the total entropy densities of the plasma 
can be written as 

Pt = 9piT)-T\ s, = g^{T)—T\ (6.66) 

For temperatures much larger than the top quark mass, all the known species of the minimal standard 
model of particle interactions are in local thermal equilibrium, then Qp = Qs = 106.75. Below, T ~ 175 
GeV the various species start decoupling, the notion of thermal equilibrium is replaced by the notion 
of kinetic equilibrium. The time evolution of the number of relativistic degrees of freedom effectively 
changes the evolution of the Hubble rate. In principle if a given mode k reenters the Hubble radius at 
a temperature the spectral energy density of the relic gravitons is (kinematically) suppressed by a 
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factor which can be estimated as (see, for instance, |137j ) 



( 



g,{T,)\/g,{T,)y^/^ 



(6.67) 



S'pO /V fi'sO / 



At the present time Qpo = 3.36 and gso = 3.90. In general terms the effect parametrized by Eq. (16.671) 
will cause a frequency-dependent suppression, i.e. a further modulation of the spectral energy density 
^Gw{^, To)- The maximal suppression one can expect can be obtained by inserting into Eq. (I6.67P the 
largest Qs and gp. So, in the case of the minimal standard model this would imply that the suppression 
(on Qq^j^IujTq)) will be of the order of 0.38. In popular supersymmetric extensions of the minimal 
standard models gp and gg can be as high as, approximately, 230. This will bring down the figure given 
above to 0.29. 

All the three effects estimated in the last part of the present section (i.e. free streaming, dark 
energy, evolution of relativistic degrees of freedom) have common features. Both in the case of the 
neutrinos and in the case of the evolution of the relativistic degrees of freedom the potential impact 
of the effect could be more pronounced. For instance, suppose that, in the early Universe, the particle 
model has many more degrees of freedom and many more particles which can free stream, at some 
epoch. At the same time we can say that all the aforementioned effects decrease rather than increasing 
the spectral energy density. Taken singularly, each of the effects will decrease ^gw by less than one 
order of magnitude. The net result of the combined effects will then be, roughly, a suppression of 
^Gw(i^; To) which is of the order of 3 x 10"^ (for 10"^^ Hz < z/ < 10"^^ Hz) and of the order of 4 x 10~^ 
for u > lO^^^Hz. These figures are comparable with the possible inaccuracies stemming from the 
calculation of the transfer function and, therefore, this is a further motivation, to use the transfer 
function of the spectral energy density. Finally the late time effects reduce a quantity which is already 
pretty small, i.e., as computed, hlflQT^{h',TQ) ~ 10~^^ for z/ ^ u^q- 
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7 B-modes induced by long wavelength gravitons 



In the minimal realization of the ACDM scenario the scalar fluctuations of the geometry induce an 
15-mode polarization which has been observed and which is now subjected to closer scrutiny [HI HJ 
El El |7j The tensor modes of the geometry not only induce an E-mode polarization but also a B- 
mode polarization. The detected angular power spectra due to the presence of a putative (adiabatic) 
curvature perturbation are the temperature autocorrelation (TT angular power spectrum) the E-mode 
autocorrelation (EE angular power spectrum) and their cross correlation (i.e. the TE angular power 
spectrum). The various angular power spectra of the temperature and polarization observables have 
been already defined in section [2] (see, in particular, Eqs. fl2.48p - fl2.49p and discussions therein). Long 
wavelength gravitons contribute not only to the TT, EE and TE angular power spectra but also to the 
B-mode autocorrelations, i.e. the BB angular power spectra. The effect of long wavelength gravitons 
on the temperature and polarization observables can be studied by deriving the evolution equations 
of the brightness perturbations which are related, in loose terms, to the fluctuations of the Stokes 
parameters. The tensor nature of the fluctuation defined in Eq. (12. 8p plays, in this respect, a decisive 
role. In particular the following two points should be borne in mind: 

• in the case of the scalar modes of the geometry the heat transfer equations have an azimuthal 
symmetry; 

• in the case of the tensor modes the fluctuations of the brightness do depend, both, upon /i = cos 
as well as upon (p; this is ultimately, the rationale for the existence of a B-mode polarization. 

The heat transfer equations can be schematically written as 



lOTT J 



.dr, 

where e' = XeUeO'TO' is the differential optical depth In the expression for the differential optical 
depth ax = (8/3)7rr2 where tq = e'^/rrie is the classical radius of the electron. The (differential) cross 
section for Compton scattering of polarized photons can be written in terms of aj- and it takes the 
usual Klein-Nishina form: 



(7 2) 

dn 32n\EiJ [Ei ' E{ " ' ^ ' r ^ ' ' 

where e and e' are, respectively, the outgoing and ingoing polarizations of the photons; Ef and E^ 
are, respectively, the energies of the outgoing and of the ingoing photons. In the limit E{ ^ E^ the 
differential cross section becomes, as anticipated 

^ = ^(e ■ ir, (7.3) 
ail 87r 

The right left side of Eq. (17. ip constitutes the collisionless term while the right hand side is the col- 
lisional contribution. At the right-hand side of Eq. (17.30 is, in general, a matrix whose di- 

mensionality depends upon the specific problem. As it will be shown Q') can be easily computed 

from Eq. fl7.3p . In similar terms f{^) should be understood as a column matrix whose components 
are the various Stokes parameters. 



^^In the differential optical depth enters not only the cross section but also the electron concentration and the 
ionization fraction Xc- The notation for the differential optical depth varies: some authors prefer k' some other f. Given 
the notations used for the conformal time coordinate we will stick to the choice made in Eq. (|7.ip . 
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7.1 CoUisionless Boltzmann equation for the tensor modes 



The coUisionless part of the Boltzmann equation can be written as: 

dl^dl dx^dl dq^dl 

dr dr dr dx^ dr dq^' ^ ' ' 

where g* is the comoving three-momentum of the photon. If the tensor fluctuation of the geometry is 
parametrized as in Eq. fl2.8p 

Si^^Qij = -a'^{T)hij{x,T), (7.5) 
then, the right hand side of Eq. (17.41) can be written as 

where the prime denotes a derivation with respect to r and where the following identities have been 
used 

dq' dq' dr \q\' dr 2^^^' ^'''^ 

In Eq. (17.71) denotes the direction of the photon momentum. The third identity appearing in Eq. 
(17.71) stems directly from the definition of comoving three-momentum, i.e. 

- 6,,pY = g,,P'P\ q' = ap' (7.8) 

where is a generic spatial component of the canonical momentum obeying 

-J ^G. 

g^pP^pP = 0, = — . (7.9) 

Using Eq. (12. Sp the full metric Qij appearing in Eq. (17. 8p becomes 

gij = -a^{T)6ij + 6[^^gij = -a^{T){6ij + hij). (7-10) 

Consequently the comoving three-momentum and its conformal time derivative can be expressed 
as: 

9' = a=(*- + ift5)pJ, (7.11) 



f ^ [V«(.- + ift-)p'4A;> + a=(.; + ift-);?^]. (7.12) 

To obtain an explicit expression for dq^/dr the conformal time derivative of the canonical momentum 
should be made explicit by using the geodesic equation, i.e. 

To first order in the tensor fluctuations of the geometry the Christoffel connection will then lead to 
the following expression: 

Using Eq. (17.140 into Eqs. (I7.12p and (17.130 . the third identity reported in Eq. (17.70 is quickly 
recovered. Depending upon the problem at hand the coUisionless part of the Boltzmann equation can 
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be written in different (but equivalent) forms. Equation (12.91) allows to write h'^ in terms of the two 
polarizations: 

A 

Using Eqs. (17. 7p and (I7.15p . Eq. (17. 6p can be written, in Fourier space, as 

^ + ik^f - ■ " ■ + 2(^ ■ a){n ■ h)h'^ = Ccou, (7.16) 

where Ccou denotes, for completeness, the collisional part of the Boltzmann equation; /i = /c ■ n is 
the projection of the photon momentum on the direction of the Fourier mode. If the direction of 
propagation of the gravitational wave coincides with the third (Cartesian) component, i.e. k = z. It 
is then easy to see that the unit vectors a and h must be directed, respectively, along the x and y 
Cartesian directions. This particular choice of the coordinate system implies then 

{n ■ af - {n ■ hf = nl - nl = {1 - /i^) cos2(^, (7.17) 
2(n ■ d){n ■ b) = 2n^ny = (1 - /i^) sin2(/?. (7.18) 

Equations (17.171) and (17. 181) are not symmetric for ip — * —ip: while Eq. (I7.17P is left unchanged, Eq. 
(I7.18P acquires a minus sign. Conversely, the evolution equations of the scalar modes of the geometry 
are symmetric for ip —ip. 



7.2 Azimuthal structure of the collisional contribution 

The differential cross section of Eq. (17. 3p depends upon e and e' which are, respectively, the polariza- 
tions of the outgoing and of the ingoing photons. If k coincides with the radial direction e and e' can 
be written, respectively, as 

k = {sin'd cos V9, sin '(9 simp, cost!}), 

ix = (cos?? cosip, cos'd simp, — sin-;?), iy = {siiup, —cosip, 0), (7-19) 
k' = (sint?' cosip', siia-d' simp', cosi)'), 

= (cos-i?' cosip', cos{}' simp', —sin{}'), e'y = {simp', —cosip', 0). (7.20) 

In Eqs. (I7.19P and (17.201) the components of e and e' have been identified with the x and y directions. 
This is just to guide the intuition since the components of e and e' should only be orthogonal to 
each others and orthogonal to the direction of propagation. Furthermore notice that, in the case 
1^' = ^' = 0, 

k' = (0, 0, 1), e', = (1, 0, 0), e'y = (0, -1, 0). (7.21) 

Under the transformation ip — ^ —(p and ip' —f', Eqs. (I7.19P and (17.201) lead to another valid choice 
of orthogonal frame. The intensity and the polarization of the (outgoing) radiation field will be given, 
respectively, as 

J(r, if) = Ix{t, ip) + Iy{T, i3, if), Q{t, i3, ip) = Ix{t, ip) - lyir, i3, ip), (7.22) 

Since Q and U are not invariant for a rotation in the polarization plane and as soon as Q is generated 
also U follows (see, e.g., Eqs. (I2.37P and discussion therein). For pedagogical purposes, it is useful to 
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consider, as a warm-up the simplification provided by Eq. (17.221) . Consider the case where the ingoing 
polarization is fixed (e. g. {}' = ip' = 0); suppose also, for sake of simplicity, that the dependence of 
the Stokes parameters upon the conformal time coordinate is trivial. Using Eq. (17. 3p the outgoing 
Stokes parameters are given by 

3e' 



6TT 

" 8^ 



a; '^x) -^X ~^ i^X S^y) 



I ^x 6^ 



{eye',)X + ie,e:fi:, 



(7.23) 
(7.24) 



In the generic situation the ingoing Stokes parameters do depend upon ip' and and, therefore, the 
outgoing Stokes parameters are 



J,(^9,y.) = ^ JdQ' 



(7.25) 
(7.26) 



From the definitions of e and e', simple trigonometric identities lead to the following expressions for 
the four products appearing in Eqs. (17.251) and (I7.26P : 



(e^ e^.) = /i' cos ((^' - (^) + 1 - /i - 
i^y eL) = -/i' sin {if' - if), {iy e'y) = cos (</?' - if). 

As already mentioned, in Eqs. (I7.27P and (I7.28P the notations 
fi = cos'd, /i' = cos-i?'. 



ex e'y) = n sin(v?' - f), 



sin 'd 



sin 



1 - f^", 



(7.27) 
(7.28) 

(7.29) 



have been used. The results of Eqs. ( I7.25p -( 17.26p and of Eqs. ( I7.27p - (l7.28p allow for an explicit 
evaluation of the the coUisional part of the Boltzmann equation which can be written, component by 
component, in the simplified case where /(r, i?, ip) has only two components which can be identified 
with the intensities of the radiation field along the x and y axes 



/(r, ^, Lp) = (X^(r, ^, Lp), Xy{T, ^,ip)). 



(7.30) 



The result is: 



dXr 



-e'Tx 



dr 

dfi' V[-Mll(/i,/x^¥^,¥^02:;(r,/i^</^') +-Mi2(/i,/i',¥',¥'')2:;(r,/x',^')] (7-31) 



IGvr J- 



dXo, 



dr 

3e' /•! r27r 

+ 



167r 

where, by definition. 



"ci<^'[^12l(/i,/x^<^,<^o^^(r,/i^<^o +^^22(/^,/x^v',<^02■;(r,/i^<^0], (7.32) 



2(e^ 



/ \2 



+ 2(1 - ^^)(1 - ^'^) +/iV^cos2((/?' - y?) 



+ ^'yl ~ fi'^yl - fi'"^ cos {if' - ip), 

A^i2(/x,/i',<^,<^') = 2{e^e'y)^ = fi^[l-cos2{^' -^)], 

M2i{fi,fi',(p,(p') = 2{eye^)^ = fi'^ - fi'^ cos 2{ip' -If), 

A^22(/U,/u', (/),(/)') = 2(ej^e' )^ = 1 + cos2(v9' - (/}). 



(7.33) 
(7.34) 
(7.35) 
(7.36) 
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The second equality appearing in Eqs. fl7.33l) - (l7.36l) follows immediately from Eqs. (I7.27P and fl7.28p . 
If T'^ and do not depend on Eqs. f l7.3ip and (17.320 can be simplified by integrating explicitly 
upon if': 

^ = -eX + 1^ £ dn'{[^?n'' + 2(1 - - /i'2)K(r, /i') + ^,X{r, /z')}, (7.37) 
(7T 3(=' /■! 

^ = + |^d/i'{[/i'2j;(r,/i') +X;(r,/i')}. (7.38) 

Equations (I7.37P and (17.381) is just a useful warm-up in view of the realistic situation where: 

• the components of /(^9, ip) are not 2 but 3, i.e. Xj.('(9, y?) and Xj^('(9, ip) are supplemented by ?7('i9, v^); 

• all the 3 components of f{{},ip) do depend upon ip; in the analog of Eqs. (I7.37P and (17.381) on 
top of the integration over n' an integration over ip' will appear. 

Since f{i), ip) has now three components 

/(r, ^, ip) = iUr, ^9, ip), lyir, ^9, ip), U{t, ^, ^)). (7.39) 

M.{VL,VL') will be a 3 X 3 matrix whose entries can be computed in full analogy with what has been 
already done in Eqs. (I7.33p - (l7.36p . In general terms, the ingoing electric field can be expressed as 

E\r, ip') = E',{t, ^') + ^(r, {}' , ip') e'^, (7.40) 

implying that the outgoing Stokes parameters can be written as 

I,{T,'d,ip) = ^ fdn'ie,-E')^, (7.41) 

Iy{T,^,ip) = ^ [dn'iey -E'f, (7.42) 

oTT J 

UiT,^,ip) = dn'[2{e, ■ E'){ey ■ E')]. (7.43) 

where, for the moment, the collisionless part of the Boltzmann equation has been neglected. Inserting 
Eq. (17.401) into Eqs. (I7.4ip -( 17.43p the explicit form of Jv[{VL,Q!) since the Stokes parameters of the 
ingoing radiation field are nothing but 

X', = E'^, T'y = E'l U' = 2E',E'y. (7.44) 

The terms M-u, J^u, M.21, M.22 will be exactly the ones already evaluated in Eqs. (I7.33p - (l7.36p . 
The remaining entries are found to be 

Miz{^, n', ip, if') = 2 (e^ e^) (e^ e'y) = /iV' sin 2((^' - ip) 



+ 2fi^l-fi^^l-fi'^sm{y^'-y^), (7.45) 
M23{f^,fi',V,¥^') = 2{e',ey){eye'y) = -fi'sm2{^' -^), (7.46) 
A^3i(/i,/x', (/?,(/?') = 4:{e'^e^){eye'^) = -2fi'^fism2{ip' - ip) 



- 4/iVl-/iVl-/i'2sin(<^'-(^), (7.47) 
A^32(/i,/i',¥P,(/^') = 4(e,e;)(ej,e;) = 2/isin2(^'-y;), (7.48) 
M33{f^,f^',<f,<f') = 2{e^ e'^){ey e'y) + 2{e^ e'y){ey e'^) = 2fifi' cos2{ip' -ip) 



+ 2^1- fi^^l- fi'^cos{ip' -^). (7.49) 
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As in Eqs. f l7.33p - (l7.36p . the second equality in each of Eqs. fl7.45l) - fl7.49l) follows immediately from 
Eqs. ( 17.271) and (17.281) . A final remark on the symmetry properties of the various entries of A4{Q, Q') 
is in order: 



M.12, ■M.211 M.22 and AI33 are all symmetric under the simultaneous transformation Lp 
—if and if' —V?'; 

for the same transformation, the remaining entries flip their respective sign. 



7.3 Different parametrizations of the full equation 

The phase space distribution obeying the coUisionless part of the Boltzmann equation (see, e. g., Eq. 
(I7.16P ) does depend upon k = \k\ (i.e. the Fourier mode). In the previous subsection the dependence 
upon k has been suppressed (for sake of simplicity) since all the aforementioned considerations could 
be separately repeated for each Fourier mode. Similarly, from Eq. (I7.16p . it is clear that the phase- 
space distribution does depend not only upon n (i.e. the direction of the photon) but also upon its 
comoving three-momentum. The phase space distribution consists of 3 independent functions whose 
dependence can be written, in Fourier space, as: 

f{k,T,fi,ip,q) = [I^{k,T,^M,ip,q),Iy{k,T,fi,ip,q),U{k,T,fi,ip,q)]. (7.50) 

As already discussed it is more handy to use directly ^ = cos^ as independent variable. The full 
Boltzmann equation can then be solved by iteration around the equilibrium configuration provided by 
the Bose-Einstein distribution fo{q); this means that the 3 components of /(r, /i, ip) can be written as: 

I^{k,T,fi,ip,q) = /o(g)[l +J^(/c,r,/i, (/?)], (7.51) 
Iy{k,T,^i,ip,q) = fo{q)[l+Ty{k,T,fi,<p)], (7.52) 
U (fc, r, yu, if, q) = foiq)U{k, r, /x, (p). (7.53) 

The solution of the problem will then require the determination of Tx{k,T, fi,ip), Xy{k,T, fj,,ip) and 
U{k, r, /i, ip). The form of the final solution for X^i^k, r, /i, ip), Xy{k, r, /i, ip) and U{k, r, /i, ip) will depend 
upon the polarization of the relic graviton. This aspect cal be appreciated by looking at the coUisionless 
part of the full Boltzmann equation, i.e. Eq. ( I7.16p . Consider, for sake of concreteness, the case of 
the © polarization. Inserting Eq. (17.500 into Eq. (I7.16P and using the perturbative scheme of Eqs. 
(I7.5ip - (l7.53p . the three independent components of the Boltzmann equation become: 

dX — 1 dh d In /q 

— ^ + {ikfi + e')X^ - -— ^cos2(y9(l - /i^)^^ — - = C^{k,T, fx, ip), (7.54) 
or 2 OT olnq 



+ {ikfi + e')Xy - --^cos2v9(l - fi = Cy{k,T, fi,ip), (7.55) 

+ {ikfi + e')U = Cuik,T,fx,ip), (7.56) 



where, for convenience, the collision terms Cx{k,T, fi,ip), Cy{k,T, fi,ip) and Cu{k,T, fi,ip) have been 
expressed as: 

3e' f — — 

Cx{k,r,iJ,,ip) = -— / dij,'dip'[Mn{fi,iJ,',ip,ip')Xx{k,T,iJ,',ip')+Mi2{fi,fi\(p,(p')Xy{k,T,ij,',ip') 
ion J 
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+Mn{fi, /x', v^, ip')U{k, r, /i', ip% (7.57) 

3e' /■ — — 
ibvr J 

+X23(/^, /i', V', V'Of^l^, ^, /^', V'')], (7.58) 

3e' f — — 

Cuik,T,iJ,,(p) = — / dfi'dip'[M3i{ix,fi',ip,(f')I.^{k,T,ix',ip') + M32{l^, fJ^' , ^, ^')Iyik,T, n', ip') 



IGtt 

+M33{fi, fi', ^, ^')U{k, T, fi', ^')]. (7.59) 

In the case of the ® polarization Eqs. f l7.54l) - fl7.56l) lead to a similar result where, however, the relevant 
part of the coUisionless contribution is different and it is given by the replacement 

-^cos2v.(l-n^-2^sm2c,(l-/.)^, (7.60) 

as it can be argued directly from Eq. (17.161) . By inspecting Eqs. (17.541) . (17.551) and (17.561) it is possible 
to argue that the azimuthal structure of the equations, appearing in connection with the polarization 
of the gravitational wave, can be decoupled from the radial structure by appropriately writing the 
various Stokes parameters. In the case of the © polarization, symmetry considerations demand that a 
sound ansatz for the full solution can be written in terms of two independent functions, i.e. I3{k, r, /i, Lp) 
and C(A;,r,/i,v9): 

J.(fc, r, ^) = ^I^^ (1 - 2^ + ^i^^(i + cos 2^, (7.61) 

J,(fc, r, /i, ^) = ^(^1^(1 - cos 2^ - ^^tlllA^i + (7_g2) 
U{k, T, /i, (/)) = 2fiP{k, T, /i) sin2(y9. (7.63) 

In the case of the polarization the azimuthal structure of the coUisionless Boltzmann equation 
changes because of the replacement of Eq. (I7.60p . The analog of Eqs. (I7.61I) - (I7.63I) become, in the 
case of the polarization, 

Uk, r, /i, ^) = ^i^^(i - ;,2) sin 2^ + ^i^^(i + gin 2^, (7.64) 

X,(A:, r, ^, V.) = ^i^I^(l - ^^) sin 2^ - ^i^^(i + 2yp, (7.65) 

U{k,T, ^,ip) = — 2/i/?(A;, r, /i) cos2(y9. (7.66) 

Using, alternatively, either Eqs. (I7.6ip - (l7.63p or Eqs. (I7.64l) - (l7.66p the exphcit evolution equations 
obeyed by P{k,T, fj.,ip) and ({k,T, ^,ip) can be derived. For sake of concreteness consider again the 
case of the © polarization and insert Eqs. dZSH), (17:^^ and fl?:^ into Eqs. fl737p . fl73g]) and (173^]) . 
In the collision terms the integrations over (f' can be performed by recalling that 

/ cos 2(^9' — </)) cos 2^9' c/i/)' = / sin 2 ((/)' — (/)) sin 2<^' (if/)' = 7r COS 2<^, (7.67) 
JO Jo 

I-2TT 

/ sin 2 ((^9' — (y9) COS 2(y9' = —TT sin 2(y9, (7.68) 

JO 

p2it 

/ COS 2(^9' — </)) sin 2(y9'(iy9' = vr sin 2(y9. (7.69) 

JO 

Consequently, Eqs. (17.570 . (I7.58p and (17.590 become 

C^{k,fx,^) = -e' fx^cos2ip J^^lCil-x^ - P{l + xy -Apx^]dx, (7.70) 
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Cu{k, /i, = -2fie' sin 2<^ J ^[({1 - - + x^ - 4Px^]dx. 
Now the essential steps of the derivation are the following: 



(7.71) 
(7.72) 



Eqs. f l7.61l) - fl7.63l) must be inserted, respectively, at left hand side of Eqs. f l7.54l) - fl7.56l) : 
Eqs. f l7.70l) - (l7.72l) must be inserted, respectively, at the right hand side of Eqs. fl7.54p - fl7.56p . 



-2ey^ (7.73) 
2e^ (7.74) 
(7.75) 



The final result of the previous pair of manipulations can be explicitly written as cj 



(1-/X 

d(3 



, , '91n/o 
^ + (.A:^ + e)C-^/., 



+ (l + /i' 



^ + {zkfi + t')(3 
or 



Or 



+ {ik^i + e')(3 = -e'^ 



where "^{k, r) is given by 

^{k,T) 



3 



32 j-i 



dx[-il + x'YP + (1 - x'YC - 4x^/3]. 



(7.76) 



As expected, one of the three equations appearing in Eqs. (17.730 - 07. 741) is redundant. Inserting Eq. 
f l7.75p into Eq. f ]7.74p the two independent equations turn out to be 



d(3 



+ {ikjj, + e')l3 = -e'^. 



OT omq 
Concerning the result obtained in Eqs. fl7.77p and fl7.78p few comments are in order: 



(7.77) 
(7.78) 



in Eqs. fl7.77l) - fl7.78l) /i^ denotes indifferently either /i'^ or h'^: the derivation reported in the case 
of /i® can be repeated in the case of bearing in mind the differences in the angular structure 
(i.e. Eqs. (EMD-dHSD should be used instead of Eqs. ([7SI])-([7l3])); 



by changing ip —Lp and if' — >• —(/?', Eqs. fl7.61l) - fl7.63l) and Eqs. fl7.64p - fl7.66p will be different 
because if the various sines appearing the various expressions; therefore the angular structure of 
the Stokes parameters will change but Eqs. f)7.77p and fl7.78p will keep their form. 



The rationale for the latter statement is that while (p —if and (p' —if' changes the angular 
structure of the Stokes parameters, also some of the entries of M.{VL^VL') will change. The net result, 
as already mentioned, will be that Eqs. fl7.77p and fl7.78p will still be vahd. The result expressed by 
Eqs. f l7.77p and fl7.78p has been firstly obtained by Polnarev |133j and then exploited for different 

^^As in the previous sections the h'^ and h'^ denote the time derivatives of the polarizations with respect to the 
conformal time coordinate t. This notation has been avoided in the previous equations of the present section since it 
could have been confused with the angular variables describing the polarizations of the outgoing photons. From now on 
this possible clash of notations does not arise. 
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semi-analytical discussions of the problem (see, e. g. |134[ 11351 1136] ). The polarization decomposition 
leading to Eqs. fl7.77p and fl7.78p can be related to slightly different treatments (see, for instance, 
|91j ) using the brightness perturbations rather than (3 and C,. The connection between the different 
formalisms will be explored later in this section. By summing up Eqs. (17.771) and (I7.78P the collision 
terms cancel: 

(7.79) 



+ [ikii + e )(C + p) — = 0. 



Defining ^ = P + ( the collision term of Eq. fl7.76p can also be written as: 

3 /•! 



"^{k, r) = - J ^ dx[{l - x^mk, r) - 2(1 + x^f (3{k, r)]. 



(7.80) 



The collision term can be further simplified by expanding in multipoles the unknowns (3{k,^,T) and 

C{k,^^,T) 

f3{k, /i, r) = E(-0'(2^ + r)Pe{f^), ({k, /x, r) = E(-^)'(2^ + ^)Uk, r)Pi{fi), (7.81) 



where Pf(/i) denote, as usual, the Legendre polynomials. By recalling that the first three Legendre 
polynomials of even order are 



P2(x) 



3x' 



PAx] 



35x^ - 30x2 ^ 3 



(7.82) 



the orthogonality of the Legendre polynomials imply that the integral over x in Eq. fl7.80p leads to 

(7.83) 



So far the discussion has been conducted in terms of one single polarization at the time. For both 
polarizations the evolution equations of ( and (3 turn out to be, of course, the same. However the 
azimuthal structure of the relevant Stokes parameters will be different for different poolarizations. In 
view of specific applications it is useful to introduce a unified notations allowing for the simultaneous 
treatment of the two different cases: 



I{k,T,fi,ip) = I^{k,T,fi,ip)+Iy{k,T,fi,(p) 



e2'%(A;)+e-2*%(A;) 



2iip 



Q{k, r, /i, V3) = I^{k, r, /i, ^) - Iy{k, r, /i, (p) 



e2^^ei(A;) + e-2*^e2(fc) 



2iip 



C{k,r,fx) 
V2 ■ 



Pik,T,fx) 



U{k,T,fi,^) 



e''%(A;) -e-"'%(fc) 



V2 

l3{,k,T,ii) 

V2 ' 



(7.84) 

(7.85) 
(7.86) 



where, following Eqs. f l5.1ip -( l5.13p (see also [91j) the stochastic variables 



ei{k) 



e@(fc) - ie^(fc) 



62 (fc) 



(i®{k) + ie^jk) 
V2 



(7.87) 



have been introduced. Note that, Eqs. fl7.85p and fl7.86p reproduce exactly Eqs. (17.610 - 07.630 (when 
e(^{k) = 0). Similarly, when e^{k) = 0, Eqs. (17. 640 - 07. 66p are correctly recovered. The variables of 
Eq. 07.87P assume (see also Eqs. 05.120 - 05.130 ) that the time dependence can be factorized by in 
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terms of an appropriate transfer function for the amplitude. By linearly combining Eqs. (17.851) and 
fl7.86p it is also easy to obtain 



Q{k,T,fi,(p) +iU{k,T,fi,<^) 
Q{k, T,^,ip) -iU{k, T,^,ip) 



;i + i^fe^'^'eiik) + (1 - fife-^'^e2{k) 
;i - ixfe''''e,ik) + (1 + ixfe-^'^e2{k) 



l3{,k,n,T) 
V2 



(7.88) 
(7.89) 



The dependence of Eqs. (I7.78P and (17.791) on the derivative of /o(g') with respect to the comoving 
three-momentum q can be further simplified by rephrasing the Boltzmann equations in terms of the 
appropriate brightness perturbations. Recalling that, by definition of brightness perturbations, 



/o 



1 + Ai 



1 - 

(9 In g 



(7.90) 



the intensity I introduced in Eq. (I7.84p can be written as X = — Ai(91n /o/c?lng). Consequently, Eqs. 
(EMI) and dTSHD-dZSnD can also be written as: 



Ai{k,T,fi,ip) -- 
A+{k, T,fi,ip) 
A_(/c,r,/i,¥?) 



(1 - ix')e'''^e,{k) + (1 - /i2)e-2^^e2(fc) 
(1 + fxfe^'^eiik) + (1 - ^)2e-2^^e2(fc) 
(1 - fxfe^'^eiik) + (1 + ^)2e-2^^e2(fc) 



At(A;, r,/i), 
Ap(A;,r,/x), 
Ap(/c, r,/i). 



where the relation of At (A;, r, /i) and Ap(A;, r, /i) to P{k, r, /i) and ({k, r, /i) is 

C{k, r, /i) = -V2^^^AT{k, r, /i), /3(A;, r, /i) = -V2^i^Ap(A;, r, /i). 



(7.91) 
(7.92) 
(7.93) 

(7.94) 



Using Eq. (17.941) inside Eqs. (I7.78I) - (I7.79I) . the evolution equations obeyed by Ap(A;, r, /x) and 
At (A;, r, /x) can be obtained and it is: 



dA^ 

S{k,T) 



+ {ikfi + e')Ap = —e'S, 
+ {tkfi + e')AT = -/i' + e'S, 



Ato At2 3 ^ \ /3 ^ 6 ^ 3 ^ 

— - H ^ At4 - -Apo Ap2 H Ap4 

10 7 70 7 V5 " 7 ^ 70 



(7.95) 
(7.96) 
(7.97) 



In Eq. (I7.96P h is the canonical amplitude of the graviton. Recalling the results of section [3] it has been 
defined that = \/2i-ph and /i® = \f2t-ph. The ^J2 factor simplifies once the brightness perturbations 
of Eqs. (I7.94P are used. Using the common practice Eqs. (I7.95p - (l7.97p have been written in units 
£p = 1. 



After Eqs. (I7.19P and (I7.20p it has been stressed that by flipping the sign of and y?' the outgoing 
and ingoing polarizations change but the choice of orthnormal frame is still valid. The discussion 
conducted so far can be repeated by assuming the polarizations of Eqs. (I7.19P and (I7.20p but with 
(y9 — i> —Lf and Lp' — i> —(/}'. The same sign flip can be followed throughout the derivation and it can be 
appreciated that the results of Eqs. (I7.95p . (I7.96P and (I7.97P are left unchanged (see also discussion after 
Eqs. (I7.48p - (l7.49p ). The flip in the sign of <y9 and affects, however, the brightness perturbations. 
The intensity of the radiation field, as expected, is left invariant when >^ —ip (see Eq. (I7.9ip ). 
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On the contrary, by repeating all the steps of the derivation, Eqs. (I7.92p and (17.931) are modified. 
Consequently, by fiipping the sign of y?, Eqs. f l7.9ip . (17.921) and (17.931) become 



Ai(fc,r,/i,v?) = 



1 - fi')e^'^e,{k) + (1 - /i2)e-2^^e2(^) 
(1 - /i)2e2^^ei(fc) + (1 + fife-^'^e^ik) 
(1 + fxfe^'^eiik) + (1 - /i)2e-2*^e2(fc) 



AT(/c,r,/i), 
Ap(fc,r,/i), 
Ap(A;, r,/x). 



(7.98) 
(7.99) 
(7.100) 



In [91] the brightness perturbations for the poolarization have been assigned as in Eqs. (I7.99p - (I7.100p 
and not as Eqs. (I7.92p and (17.930 . The ladder operators of Eqs. (I2.52p and (I2.53P are not invariant 
under the transformation if — > —if. More specifically the spin weight of a given function changes the 
sign whenever ip — »• —cp. Equations (I2.52p and (12.530 are consistent with the brightness perturbations 
written as in Eqs. (17.990 and (I7.100p . To be consistent with the customary notations we will therefore 
adhere to the conventions stipulated in \91\ . 



7.4 B-modes from relic gravitons 

The analytic and numerical solutions of the evolution equations for the brightness perturbations allow 
for an explicit evaluation of the temperature and polarization observables. The results obtained so 
far imply that long wavelength gravitons will produce both temperature and polarization anisotropies. 
More specifically, following the terminology of section [2] (see, in particular, Eqs. (12.480 - 02.500 ) the 
relevant angular power spectra induced by the relic gravitons will be the TT, EE, TE and BB angular 
power spectra. 

The explicit form of a\J^ and a}^ can be derived from the general expressions already encountered, 
respectively, in Eqs. (I2.72p and (12.730 . The integrands of Eqs. (I2.72p and (12.730 can be computed 
from the action of the the differential operators of Eqs. (12.520 - 02.530 . i.e., more specifically. 



'(2), 



dlA, 



+ 2i 



d^A'^ 



1 — fi 



-,d^A+ 



K^-^\n)[K'x^\n)A_{n)] = (1 - /x2)A': - 4/xA'_ - 2A_ 



2i 



d^A'_ 



1 — H 

Equations (I2.72p and (12.730 become then 



(7.101) 



(7.102) 



-N, J dnY;^{n)[{l - ^2)A'^ - 4/xA'u - 2Au - ^ 



dlA^ 



+ 2 



fx 



-N, J dnY;^{n)[{l - f?)A'l^ - 4/iA'Q - 2Aq - ^ 



_Q 
/z2 



/i 



1-/X2 



:d^Aij 



(7.103) 



(7.104) 
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In Eqs. f l7.103l) - (l7.104p the prime denotes a derivation with respect to /i = cos'd. This notation will 
be consistently followed, for sake of simplicity, but only in this subsection. Equations (17.1031) and 
(I7.104P can be easily written in even more explicit terms and, as before, it is easier to focus on a single 
polarization, e. g. ©: 

Aq(A;, ra,r) = V2{1 + /i^) cos 2ip Ap {k, ij,t), Au(fc, ri,r) = -2V2 jj, sm2ipAp{k, ij,t). (7.105) 
Inserting Eq. (17. 1051) into Eqs. (I7.103P and (I7.104p the resulting expressions will be 

- 2V-2N, J dnY;jn) =i„2^(l - ,^)|,AJ + 4Ay, (7.106) 



"t m 
.(E) 



(7.107) 



In Eqs. (I7.106P and (I7.107P the explicit integration over the angular variables can be performed in 
explicit terms. Consider, in particular, Eq. (I7.106p . i.e. 

aS = 2V2N^ dfi d^Y;^{fi, ^) sin 2^ (1 - /i^) [^A'^ + 4A'p]. (7.108) 

Bearing in mind the connection between spherical harmonics and the associated Legendre functions, the 
^£m(^) will be essentially given by the product of P™(/i) and of e*"*"^, up to a well know normalization 
coefficient. In Eq. (I7.108P the integration over can be carried on in explicit terms and since there 
is a sin2<^ in the integrand the whole integral will be proportional to {6m, 2 — 6m -2)', in more precise 
terms the integration over (p will bring Eq. (17. lOSp in the form: 



(£-2)! 2i+l 
(£ + 2)!V 47r 



{Sm,2-6m-2) I pf (/x) [(1 - /x')/xA;i + 4(1 - /i^) A^,] . (7.109) 



Equation (I7.109P can be further simplified by using the following three relations: 

Ap(fc, /i, r) = J2{2j + l)(-z)^ Ap,(A;, r) P,(/i), 
j 

(j - m + l)Pp^,{fi) + (j + m)P-L,{fi) 



(2j + 1) 
P,-i(/i)-P,-,(/i) 

2j + 1 



(7.110) 

(7.111) 
(7.112) 



Equations (17.1 lip and (I7.112p are well known recurrence relations for the for the associated Legendre 
functions (see p. 24 of Ref. [HS] for a derivation). Equation (I7.110p can be written in slightly different 
ways. Some authors do not include, at the right hand side of Eq. (I7.110p . the factor {—iY inside the 
sum (see for instance [93]). The analog of Eq. (I7.109P can also be obtained within the approach of 
Ref. [93] whose conventions (for the © polarization) are 



To 



P^\k,n) = ^Y.i'^j + l)A^^{k)Mf 



(7.113) 



with 



Mg^(ri) = P,(/i) 



/ (1 + /i2)cos2(/? -2 



sm 



: sin 2(p 



cos 2ip 



2^ \ 



(7.114) 
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For X polarization, the angular structure of Eq. (17.1141) will be different and can be obtained from 
the results for the © polarization by replacing cos2y9 sm2ip and sin2(y9 —cos2ip. Eqs. (17.1131) 
and ( I7.114P are nothing but the explicit expression of the polarization tensor on the sphere introduced 
in Eq. (12.811) with the difference that the indices (in Eq. (I2.8ip ) are both covariant while in Eq. 
(17.1131) they are both contravariant. The conventions of [93] imply also a difference of a factor of 1/8 
in comparison with the conventions adopted here. Equations (17.1 lip and (I7.112p can be used inside 
Eq. (I7.109p . As an example the term (1 — /i^)/iAp, after using Eq. (I7.110p will produce a factor 
/i(l-/i')i^", i.e. 

where the first equality follows from the definition of the associated Legendre functions (i.e. Eq. (12.550 ) 
while the second follows from Eq. (17.1 lip with m = ±2. The second term in the integrand of Eq. 
f lTTOD . i.e. (1 - fi^) A'p will produce a factor (1 - fi'^)P^, i.e. 



where, again, the first and second equalities follow, respectively, from Eqs. ( ]2.55p and Eq. ( ]7.11ip 
both in the case m = ±2. Inserting Eqs. (I7.115P and (17.1160 inside Eq. (I7.1U9P the integration over 
/i can be performed by simply using the orthogonality of the associated Legendre functions, i.e. 

Pr(f^)mf^) = ^ 1^ ^n'- (7.117) 
The final result for the B-mode will then be: 

.(B) 



\ 27TI ("'^'^'^"^'2 - V-2)[Ap,_i -{i- l)Ap,+i]. (7.118) 

If the factor {—iY would be absent from the expansion of Eq. (17. HOD there would not be {—i)^ at 
the right-hand side of Eq. (17.1180 : for the same reason, the relative sign of the two terms inside the 
squared bracket would be plus instead of minus. Finally, the factor (1/8) appearing at the right-hand 
side of Eq. (I7.113P will also modify slightly the prefactor of Eq. (I7.118P which will be 



^ ?' / 27r r ~ ~ 1 

a(^)'®(fc)z™ = --{6^2 - [{i + 2)A^,_,{k) + l)A^e+,{k)] , (7.119) 

where Eq. (12.1020 has been used to relate and af2.- In the case of the E-mode, always working 
with the © polarization, Eqs. (I7.104p and (I7.107p will give, after integration over ip, 



where Eqs. (12.540 and (12.550 have been used. To complete the calculation it is necessary to perform 
the integration over fi. This can be done, as previously shown, by using wisely the recursion relations of 
the associate Legendre functions. From Eq. (I7.120p we have that, using the expansion of Eq. (I7.110p . 

J'^ d/i(l - /i2)[(l + fi^)A'^ + 8/iA'p + 12Ap] = 

Y: Ap,(-^)^ (2j + 1) £ Pf H/i)(l - /i^) [(1 + fi')P;' + 2f^P; + 12P,] rf/i. (7.121) 



70 



As in the case of the B-mode, the main idea of the calculation is to express the relevant part of the 
integrand appearing at the right hand side of Eq. (17.1211) in terms of products of associated Legendre 
functions with the same m but different values of j so that the orthogonality relation (17.1170 can be 
exploited. Using the definition of the associated Legendre functions of Eq. (12.551) it is easy to show 
that 



(1 - /i')[(l + /i')P/ + 2/iP; + 12P,] = (1 + /i2)Pf ^ - S/i^/l - f^^Pj^^ + 12(1 - fi^)Pj, (7.122) 

where the minus sign in the second term at the right hand side of Eq. (I7.122p arises since the Condon- 
Shortley phase has been included in the definition of the associated Legendre functions (see Eqs. 
(I2.54p - (l2.55p and discussion therein). The right hand side of Eq. (I7.122p now be simplified by using 
the following steps: 

• Eq. (I7.11ip can be used to simplify the term /x^Pj^'* (/i) ; 

• Eq. (I7.112P can be used to simplify the term 8/xa/1 — n'^Pj^\n); 

• the term (1 — fi'^)Pj{fi) can be simplified by using the equation of the Legendre polynomials and 
the recursion relations (17.1 lip and (I7.112p . 

Using these three steps, a rather lengthy but straightforward algebra leads to the following expression: 

(1 /.)[(l+/x)P,. +2/xP,. + 12P,J - ^2_^._^^^2j + 3)^^- ^^^+(2j + l)(2j + 3)^^+^^^^ 

(j-3)(j-2) (2) 
+ (2j + l)(2,-l)^-^^^^- ^^-^^^^ 

To get to Eq. (I7.123P it is useful to recall, on a side, that 

U 1-^)1^,^) (2j + l)(2j + 3) (2j-l)(2i + 3) + (2j + l)(2i-l)' 

which can be derived by recalling the equation obeyed by the Legendre polynomials Pj{fi) as well as 
the recursion relations reported in Eqs. fl7.llip - fl7.112l) . Using Eq. (ILIMI) inside Eqs. fl7.120p - fl7.121l) 
it is easy to obtain, after simple algebra, that 



6(£ -!)(£ + 2) 

-Ap< 



{£ + !){£ + 2) , 



(2£- l)(2£ + 3) 

(7.125) 



(2£ + l)(2£-l) (2£+l)(2£ + 3) 

As in the case of the B-mode the presence of the factor (— «)•' in Eq. fl7.110p leads to a sign difference 
in the terms Ap^±2 appearing in Eq. fl7.125p . The absolute values of the coefficients of the three terms 
appearing (inside the the square bracket) in Eq. fl7.125p are independent on the conventions related to 
Eq. (17. nop . The coefficient of the term Ap^ appears to be different from the homologue term reported 
in Eq. (4.41) of [HS]- After many checks and cross-checks of the derivations leading to Eq. fl7.125p it 
has been concluded that the difference (i.e. 6£(£+ 1) instead of 6(£ — 1)(£ + 2) in the numerator of the 
coefficient of Ap^) is just the result of a typo and the correct expression is 6(£ — l)(i' + 2) as in Eq. 

f Tmsl) . 
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7.5 Angular power spectra induced by long wavelength gravitons 



The tensor modes of the geometry produce both, temperature and polarization anisotropies. Consider, 
for simphcity, the ACDM parameters determined from the best fit of the WMAP 5-yr data alone. These 
parameters have been reported in Eq. (\2.7\\ . In Fig. [10] the TT angular power spectra stemming from 
the standard adiabatic mode and the temperature autocorrelations induced by the tensor modes are 
illustrated. In both plots of Fig. [10] the full line denotes the TT angular power spectra corresponding 



n =0.742, n =0.214, 
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Figure 10: The temperature autocorrelations of the standard adiabatic scenario (full line) are compared 
with the TT angular power spectra induced by the tensor modes. 

to the best fit of the WMAP 5-yr data alone. In the plot at the left the dashed and dot-dashed lines 
illustrate, respectively, the cases rx = 0.1 and rx = 0.2. The tensor spectral index is fixed according 
to Eqs. (11.21) and (I1.6P (see also Eqs. (16.41) and (16. 5p ). In both plots of Fig. [10] the spectral index 
does not run, i.e. ax = in Eq. (II. 6p . Always in Fig. [10] (plot at the right) the dashed line denotes 
the tensor contribution in the case rx = 1 while the dot-dashed line denotes the total TT angular 
power spectrum. The temperature and polarization observables can be obtained numerically (see, e.g. 
[911 [93]) but useful analytic results do exist in the literature [1331 11341 1135[ [T36| I137j . An essential step 
for both approaches is to derive the angular power spectra in more explicit terms and as a function 
of the solution of the heat transfer equations. The solution of Eqs. ( ]7.95p . ( 17.96^ and (17.97P can be 
formally written by using the integration along the line of sight which is also common to the scalar 
case (see, e.g. |119j . for an introduction): 

/■TO 

Ap{k,fi,To)= c/rSp(A;,r)e^*^^ (7.126) 
Jo 

Ax(A;,/x,ro)= T rfr5x(A;, r)e-*^^ (7.127) 
Jo 

where x = /c(ro — r). Note that S^ik, r) and Sp{k, r) are related to the source terms of, respectively, 
Eqs. (17:951) and ^Mj, i.e. 

Sp{k, t) = -/C(r) S{k, r), S'x(fc, r) = -h'e-<^'^°^ + 1C{t) S{k, r), (7.128) 

where S{k, r) is given by Eq. (I7.97P and where 

/C(r) = e'e '^'^'^''^oj^ ^^^^ ^^-j _ / aneXea^dr, e = aheXeCTT- (7.129) 
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is the visibility function expressed in terms of the differential optical depth e' and in terms of the optical 
depth e(r, tq). The visibility function /C(r) gives that probability that a photon is last scattered between 
r and r + dr. In loose terms the visibility function is peaked around the redshift of recombination. In 
analytic discussions the visibility function is often approximated by means of a Gaussian profile (see, 
e. g., [nni [1371 [nn] and also [1381 [ffl [HQl fUT] ): 



/C(r) =U{a,,,) 



(t — Tree) 
0^ 



/C(r)rfr = 1. 



(7.130) 



As indicated in Eq. fl7.130p . A/'(crrec) is determined by requiring that the integral of /C(r) over r is 
normalized to 1. The WMAP data suggest a thickness (in redshift space) Azrec — 195 ± 2 [8j which 
would imply that (Tree, in units of the (comoving) angular diameter distance to recombination, can be 
estimated as a^cc/To — 1-43 x 10^'^. When tq ^ t^-^c and tq ^ a^^c the normalization appearing in 
Eq. (17. 1301) can be estimated as A/'(cTrec) 



,-1 



(Tree Y2/7r. In |137j it has been suggested that a better 
approximation, for the case of the tensor modes, is to assume that /C(t) has two different widths, 
respectively, for r < Tree and for r > Tree- Recalling now Eq. fl2.44p and (12.491) the af^, can be written 
as 



rfnr;„Ai(n) 



(27r)3/2 



27r 



d'k I dip I c//i(l-/i')[e2*'^ei(A:) + e-'*^e2(A;)]F;„(/i,(/^)AT(A;,/i, To), (7.131) 



2i(p 



where the second equality follows from the Fourier transform of Ai(n) and by using Eq. f l7.98p in 
the obtained expression. The angular integrations can now be performed with the approach already 
exploited in previous subsection for the polarization observables. It is however useful to solve Eq. 
fl7.13ip with a slightly different method (see, e. g. [91]) where explicit use is made of the solutions 



fl7:T26l) and ^JSH). Indeed, inserting Eq. fTTM into Eq. fTrml) a£ can be written as 
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Figure 11: The polarization autocorrelations induced by the tensor modes are compared among them- 
selves and with the EE angular power spectra produced by the adiabatic mode. 



''em 



^ ''d^k r dTST{k,r)[lP{i,m,x)ei{k)+I^J\i,m,x)e2{k)], {7.132) 



(27r)3/2 



2tt 



(7.133) 
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In X^'' the integration over is trivial since the (ordinary) spherical harmonics depend upon as 
^irrup j^ggg^ g_ -gqg_ f|2.54l) and (12.551) ). Thus, as expected on the basis of the results of the previous 
subsection, Zj_ oc 2n5m,±2- The latter results allows for a simplification so that, for instance. 



jf^(£,m,x) = 2n5, 



'2£ + 1 



m, 2 



-A7T(-iY6:, 



An \\ {£ + 2)1 J-i ^ ^ dfi^ 

+ dinx%ix)], (7.134) 



I2i- 



m, 2 



An \{i + 2)\ 



where the second equality follows by first integrating by parts and by then expanding e~*^^ in series 
of Legendre polynomials. Indeed, in Eq. (I7.134p . je{x) are the spherical Bessel functions. Note, finally 
that, inside the integral of Eq. (17.1341) expressions like (1 — /i^)^e~*^'^ can be traded for (1 + 9^)^ e~^^^. 
Recalling the properties of the associated Legendre (see second relation of Eq. (12.551) ) the other 
integral, i.e. X_ (£, m, x) can be performed exactly in the same way. By making explicit the ensemble 
averages and by using Eq. (12.491) . the angular power spectrum illustrated in Fig. [10] can be obtained 
and it is given by: 

= 4vr/^PT(A:)|AT.(A:,ro)p, 

AT,(fc,ro) = Ne r S^ik^r)^-^, (7.135) 
Jo x^ 



where, as in Eq. (12.461) . A^^ = ^J{£ — 2)\/{£ + 2)\. Having discussed the temperature autocorrelation 
induced by the long wavelength gravitons it is now the moment of discussing the polarization observ- 
ables. In Fig. [TT] (plot at the left) the EE angular power spectra (full line) are compared with the 
B-mode autocorrelations (dashed line) induced by the tensor modes in the case rx = 0.1 while the 
other cosmological parameters are fixed to the best-fit values of the WMAp 5-yr data alone. Always 
in Fig. [11] (plot at the right) the full line illustrates the B-mode autocorrelation in the case rx = 0.1 
while the dashed line illustrated the EE angular power spectrum stemming from the standard adia- 
batic mode. In Fig. [12] (plot at the left) the BB angular power spectra are illustrated for different 
values of tt compatible with current bounds on (see Tabs. [I] and [2]). 

The E-mode and the B-mode angular power spectra can be obtained from Eqs. (I7.99P and (17.1001) 
by carefully following all the steps leading to the temperature autocorrelation of Eq. (I7.135p . The 
crucial difference will be, of course, that and arise, respectively, in the expansion of A-^ifi) 
and Ab(5^) as reported in Eq. (12.460 . To compute and in terms of Sp{k, r) the steps are, in 



short, the following: 

• and should be first expressed in terms of A-^in) and Asin) as in Eq. (I2.46P : 

• then Eqs. (17:991) and fITJOOD should be inserted into Eqs. (12351) and ^M)] 

• the evolution for Ap(A;, /i,r) can then be solved in Fourier space as in Eq. (17.1260 . 

The first and the second step has been already (partially) performed in the previous subsection and 
the results have been given in Eqs. ( I7.103p -( l7.104p with the difference that, now, the two polarizations 
can be treated simultaneously so that, in Fourier space, 

AQ{k,T,fi,^) = (l + /i2)[e2^'^ei(^) + e-2^^e2(^)]Ap(fc,r,/i), (7.136) 
Au(A;,r,/i,(/?) = 2fii[e^''^ei{k) - e-^''^e2{k)]Ap{k,T, fi). (7.137) 
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Note that Eqs. (\Tim and (171371) follow directly from Eqs. and (\TJm by definition of 

A±{k,T, fj,,(p). It should be appreciated that one derivation with respect to ip changes the azimuthal 
structure of Aq(A;, r, /i, and A\j{k,T, fi,ip), i.e. 



5^Au = -Me^'^e,{k) + e-^'^e2{k)]Ap{k,T,fi), 
9^Aq = 2i{l + fi')[e^'"^ei{k)-e~^'''e2{k)]Ap{k,T,fi). 



(7.138) 
(7.139) 



Using now Eqs. (I7:T36D - (17137]) and (171381) - (171391) inside Eqs. (17l03|) and (17104]) we do get, after 
Fourier transform, 



.(B) 



-2i- 



d^k 



2tt 



^"^ ^''(27r)3/2j JO 

X il-fi^)[fidl + Ad,]Apik,fi,To) 



dip I ^ dfi[e''^ei{k) - e 



-2i(p 



(7.140) 



,(E) 



''tm 



(27r)3/2 



d'^k 



2tv 



dip I ^ dfi[e^'^e,ik) + e-^'^e2ik)]Y;^ifi,ip) 
X {l-fi')[{l + fi')dl + 8fxd^ + l2]Apik,fx,To). (7.141) 
According to Eq. (I7.126P Ap{k, fi,TQ) can be related to Sp{k,T) and the resulting coefficients are: 

Np 



(27r)3/2 



2tt 



-1 



dip / di2Y;^{ij,ip) / d^kle^'^'eiik) - e-2'^e2(A;)] 



,(E) 



X (1 - fi^)OB{x) / Sp{k, r) e-^'^^dr 



(7.142) 



dnY;^{fi,^) / (i3A;[e2^'^ei(fc) + e-2^'^e2(^)] 



(27r)3/2 

X (l-/i2)OE(x) / 5p(A;,r)e-*'^^'cir, 



(7.143) 

where Ob{x) and Oe(3^) are two differential operators, i.e. 

Ob{x) = 2x'^dl + 8x, OE{x) = x\l-dl)- 8x0^-12. (7.144) 
Now Eqs. (I7.142P and (I7.143P can be used into Eqs. (12.480 and the angular power spectra become: 

dk„ " (7.145) 

(7.146) 

(7.147) 
(7.148) 



where 



AE^(A;,ro) 
ABe{k, To) 



Cr^ = 4nJ^VT{k)\ABe{k,ro)\^' 
Cr^ = 4vr/^PT(A:)|Arf,ro)p 



4 2 2 

-dx]t{x) + -i,ii{x) + dje,{x) -je{x) 

4 



drSp^k, r). 



2d^ji{x) + -irjiix) 







drSpi^k, t) 



The tensor modes also induce a TE angular power spectrum (i.e. a cross-correlation between the 
temperature and the E-mode polarization). In Fig. 112] (plot at the right) the absolute value of the TE 
(tensor) angular power spectrum is compared with the corresponding angular power spectrum induced 
by the standard adiabatic mode. The cross-correlation between temperature and polarization can be 
obtained inserting Eqs. (I7.133P and (17.1410 into Eq. (12.490 and the result is 

Cl™) = 47r| ^PT(A:)AEKfc,ro)AT,(fc,ro), (7.149) 
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Figure 12: The B-mode autocorrelation for different values of tt (plot at the left). The TE cross 
correlation induced by the tensor modes is compared with the corresponding angular power spectrum 
induced by the standard adiabatic mode. 

where, ATi^kjTo) and AEeik,To) are given, respectively, by Eqs. fl7.135p and (I7.147p . It should finally 
be mentioned that the derivatives with respect to x appearing in Eqs. (I7.147P and (I7.148P can be 
transformed into derivatives with respect to r by making use of integration by parts [HI] . This step is 
carried on in full analogy with what happens with the scalar modes of the geometry |119j . 

Various experiments provided, so far, direct limits on the B-mode polarization. The limits on 



Experiments 


Years 


Pivot frequencies 


Upper limits on BB 


Dasi [211 Ea EH] 
Boomerang |142j 
Maxipol [1431 1144] 
Quad [281 [291 [30] 

Cbi [261 [27] 
Capmap [23 [25] 
Wmap P [1 [6] 


2000- 2003 
2003 
2003 
2007 

2002-2005 
2004-2005 

2001- 2006 


26-36 GHz 
145 GHz 
140 GHz 
100/ 150 GHz 
26-36 GHz 
35-46/84-100 GHz 
23-94 GHz 


£(£ + l)Cf^^(27r) < 50 (/iK)2 
£(£+l)Cf''V(2vr) <8.6(/iK)2 
£(£ + l)Cf^V(27r) < 112.3 (/iK)2 
£(£ + l)Cf''V(27r) < 10(/iK)2 
iii + l)Cf''V(2vr) < 3.76 K)^ 
^(£+l)Cf''V(27r) <4.8(/iK)2 
+ l)Cf''V(2vr) < 0.15 (/iK)2 



Table 3: The main polarization experiments, their typical observational frequencies, and the upper 
limits on the B-mode polarization are illustrated. 



reported in Tabs. [T] and [2] follow from a combined analysis of the TT, EE and TE angular power 
spectra which allows for a tensor contribution. As Tab. [3] indicates the upper limits on the B-mode 
polarization are still rather loose and, often, derived on the basis of a limited range of harmonics. 
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8 High-frequency spikes in the relic graviton background 



In the ACDM paradigm long wavelength gravitons can affect the CMB polarization. As the frequency 
increases towards the region accessible to wide band interferometers, the ACDM signal can only de- 
crease for number of independent reasons: 

• in the ACDM paradigm the spectral energy density is only approximately scale-invariant (see, 
e.g. Fig. [2]) but the scaling violations always tend to make the spectral energy density smaller 
at high frequencies; 

• there are various secondary effects associated, for instance, with the variation of the effective 
number of relativistic degrees of freedom, with the neutrino anisotropic stress and with the 
transition to a phase dominated by the dark-energy contribution: all these features reduce the 
spectral energy density in different frequency regions; 

• the addition of supplementary scalar fields driving the inflationary phase does not change the 
two previous statements. 

The previous conclusions are all based (either directly or indirectly) on the assumption that the 
thermal history of the Universe is minimal in the sense that, after inflation, the Universe soon becomes 
dominated by relativistic particles so that the sound speed of the plasma soon reaches values Cst = 
The post-inflationary thermal history might not be minimal. For instance it could happen that 
the transition to the radiation-dominated regime is not instantaneous [75] • More specifically it can 
happen that, after inflation, the sound speed of the plasma is such that Cst > l/"\/3. When the sound 
speed is larger than 1/a/3, the fluid is said to be stiff. In the system of units used in the present paper 
the speed of light is such that c = 1. A natural upper limit for the sound speed is exactly 1 which 
is the maximally stiff fluid compatible with causality |124j (see, however, also |125j ). If the thermal 
history of the Universe contemplates a post-inflationary phase stiffer than radiation, a spike in the relic 
graviton spectrum is expected at high frequencies [75] (see also [761 EZ] as well as Eqs. fl6.29l) - fl6.30l) 
and discussion therein). The possibility of having a post-inflationary phase stiffer than radiation has 
been also investigated in different contexts such as in |146l 11471 1148] . 

In the early Universe, the dominant energy condition might be violated and this observation will 
also produce scaling violations in the spectral energy density [149] . If we assume the validity of the 
ACDM paradigm, a violation of the dominant energy condition implies that, during an early stage of 
the life of the Universe, the effective enthalpy density of the sources driving the geometry was negative 
and this may happen in the presence of bulk viscous stresses [149j (see also [150[ 1151] for interesting 
reprises of this idea). In what follows the focus will be on the more mundane possibility that the 
thermal history of the plasma includes a phase where the speed of sound was close to the speed of 
hght. 

Absent any indirect tests on the thermal history of the Universe prior to the formation of light 
nuclear elements, it is legitimate to investigate situations where, before nucleosyntheis, the sound speed 
of the plasma was larger than 1/v^, at most equalling the speed of light. In this plausible extension 
of the current cosmological paradigm, hereby dubbed Tensor-ACDM (i.e. TACDM) scenario, high- 
frequency gravitons are copiously produced [101 SI]. Without conflicting with the bounds on the 
tensor to scalar ratio stemming from the combined analysis of the three standard cosmological data 
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sets (i.e. cosmic microwave background anisotropics, large-scale structure data and observations of 
type la supenovae), the spectral energy density of the relic gravitons in the TACDM scenario can 
be potentially observable by wide-band interferometers (in their advanced version) operating in a 
frequency window which ranges between few Hz and few kHz. 

The presence of a stiff phase increases the spectral energy density for frequencies larger than a 
pivotal frequency which is related to the total duration of the stiff phase. If the stiff phase takes 
place before BBN, then z/g > 10~^ nHz. If the stiff phase takes place for equivalent temperatures larger 
than 100 GeV, then > /iHz. If the stiff phase takes place for T > 100 TeV, then > mHz. The 
frequency Us marks the beginning of a branch of the spectrum where the tilt of the spectral energy 
density is blue (i.e. increasing in slope) rather than nearly scale invariant or slightly red (as it is the 
case in the conventional scenario). 



8.1 Scaling violations 

In the ACDM paradigm, for frequencies larger than z/gq the spectral energy density of the relic gravitons 
is, approximately, scale invariant (see Fig. [2]). In the context of the TACDM scenario the approximate 
scale-invariance of the flat plateau is violated. This situation is illustrated, for instance, in Fig. HI 
In the present subsection a more detailed account of the typical frequencies of the problem will be 
presented. If there is some delay between the end of inflation and the onset of radiation the maximal 
wavenumber of the spectrum will be given by: 



Op 



where a = 2/[3{wt + 1)] (with Wt > 1/3) is related to the specific kind of stiff dynamics. Equation 
(18.11) can also be written as 



where 

In the case E = 0{1) (as it happens in the case a = 1/3 if the initial radiation is in the form of 
quantum fluctuations) i^max = ^max/(27r) ^ 100 GHz, more precisely: 

= 1.177 xlO"E-'(jJ^)'" Hz. (8.4) 

On top of the standard parameters of the ACDM scenario (see, e. g., Eq. (12. 7p ) the minimal TACDM 
scenario demands two supplementary parameters 

• the frequency defining the region of the spectrum at which the scaling violations take place; 

• the slope of the spectrum arising during the stiff phase. 

The frequency Ug can be dynamically related to the frequency of the maximum and, consequently, the 
first parameter can be trated for S. As it will be shown in a moment, typical values for S range between 
0.01 and 0.5. The slope of the spectrum during the stiff phase depends upon the total barotropic index 
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and can therefore be traded for Wt- The curvature scale Hj. determines kg (or i/g), i.e. the frequency at 
which the spectral energy density starts increasing. Supposing that from the end of inflation there is 
a single stiff phase (as it is natural to assume in a minimalistic persepective) the value of kg is 



Using the relation of to T,, Eq. (18. 5p can also be written as 
and as 

z/s = 1.173 X io^^j:m-2a) ( ^o^Ro Hz. (8.7) 

^ ^ V4.15 X lO^V ^ ^ 

The quantity S which parametrizes the location, in the spectrum, where the scaling violations 
appear must be smaller than 1 or, at most, of order 1. This is what happens within specific models. 
For instance, if the radiation present at the end of inflation comes from amplified quantum fluctuations 
(i.e. Gibbons-Hawking radiation), quite generically, at the end of inflation pr ~ H'^. More specifically 

In Eq. (18. 8p A^ofr is the number of species contributing to the quantum fluctuations during the quasi-de 
Sitter stage of expansion. In |128j (see also [ISIEZIES]) it has been argued that this quantity could be 
evaluated using a perturbative expansion valid in the limit of quasi-conformal coupling. It should be 
clear that A'^efr is conceptually different from the number of relativistic degrees of freedom Qp. Given 
H and Nes the length of the stiff phase is fixed, in this case, by [78] 



a-. 



4 



where we used the fact that a = 2/[3{w + 1)] and where we defined A = A'^efj /(4807r^) . Equation (18.91) 
implies that 



Recalling that H,. = H{a^/ a^Y'"- ^ we also have that Eq. (I8.10p implies 



\W^)( \ 8.11 



Mp/ VMi 



'p 



Using Eq. (18.111) and the definition of S (see Eq. (18. 3p ) it turns out that S = A^/"^, which is always 
smaller than 1 and, at most, 0{1). 

Instead of endorsing an explicit model by pretending to know the whole thermal history of the 
Universe in reasonable detail, it is more productive to keep S as a free parameter and to require that 
the scaling violations in the spectral energy density will take place before BBN. Consequently the 
variation of S, w and rx can be simultaneously bounded [101 E]. If the stiff dynamics takes place 
before big-bang nucleosynthesis, then > ^hhn (see also Eq. (I6.6ip ). This requirement guarantees 
that the stiff dynamics will be over by the time light nuclei start being formed. In a complementary 
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Figure 13: The bounds on z/g in the (S,Wt) plane; the spectral energy density of the relic gravitons in 
the TACDM scenario. 

approach one might also require that z/g > z/ew where z^ew corresponds to the value of the Hubble rate 
at the electroweak epoch, i.e. 

( 9p y/y 

V 106.75/ UOO GeV. 



2.5764 X 10' 



0"R0 



3.998 X 10" 



4.15 X 10-5 



1/4 



Mpc-^ 



(8.121 



1/4 



Hz. 



V106.75/ V 100 GeV/ V4.15 X 10-5. 
Finally, yet a different requirement could be to impose that u > z/tcv where z/xev is defined as 



Z/TeV = 4.819 X 10" 



9p 



Hz. 



(8.13) 



(8.14) 



.228.757 VlOO TeV J V4.15 x 10-5, 
The latter requirement would imply that the stiff age is already finished by the time the Universe has 
a temperature of the order of 100 TeV when, presumably, the number of relativistic degrees of freedom 
was much larger than in the minimal standard model (in Eq. (18.141) the typical value of gp is the one 
arising in the minimal supersymmetric extension of the standard model). 

In Fig. [12] (plot at the left) the constraints on S, Wt are illustrated. The value of S controls the 
position of the frequency at which the nearly scale-invariant slope of the spectrum will be violated. 
The barotropic index Wt is taken to be always larger than 1/3 and with a maximal value of 1. The 
shaded region corresponds to the region excluded in the most constraining case, i.e. the one demanding 
> ^hhn- Also the values of tt are bounded by the same kind of considerations. Indeed, z/g depends 
upon e which is related to tt (see, e.g., Eq. (18. 7p ). Therefore, a lower bound on z/g also implies a 
bound in the (rx, Wt) plane. 



8.2 Spectral energy density in the minimal TACDM scenario 

In Fig. [13] (plot at the right) two examples of the scaling violations on the spectral energy density are 
illustrated. Both examples are compatible with the bounds illustrated in the plot at the left. Similar 
examples have been already illustrated in Fig. H] These examples will now be discussed in greater 
detail. The spectral energy density has been computed by using the numerical approach presented in 
section [6] 
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In the two examples of Fig. [13] (plot at the right) the ACDM parameters are fixed to the values 
reported in Eq. (12. 7p (see [21 HI O [6l [7]). The spectral index has been allowed to run, i.e. ax 7^ 
(see Eqs. (II. 6p and (16.530 ). The two supplementary parameters should be identified with the sound 
speed during the stiff phase (i.e. Cst) and the threshold frequency (i.e. z/g). Besides Cgt and z/g, there 
will also be rx which controls, at once, the normalization and the slope of the low-frequency branch 
of the spectral energy density. At the moment wide band interferometers have sensitivities which are 
insufficient for cutting through the phenomenologically interesting region [631 EH [65]. In the near 
future, however, there is the hope of a dramatic improvement of the sensitivity: even 5 or 6 orders of 
magnitude at least heeding the original design (see e.g. [55]) together with the recent proposals for an 
advanced Ligo program. 

As specifically discussed in Eqs. f l8.5p -( l8l6l) the frequency of the elbow, i.e. fg, is fully determined 
by S (see Eq. (18.31) and discussion therein). The two supplementary parameters z/g and Cgt can be 
traded for S and Wt as already done in Fig. [13] (plot at the left). In doing so there is also a potential 
advantage since, according to Eq. (18.40 . S shifts the maximal frequency of the spectrum. 

As soon as the frequency increases from the aHz up to the nHz (and even larger) the spectral 
energy density increases sharply in comparison with the nearly scale-invariant case where the spectral 
energy density was, for u > nHz, at most O{10~^^). In the case of Fig. [13] the spectral energy density 
is clearly much larger. The accuracy in the determination of the infra-red branch of the spectrum is 
a condition for the correctness of the estimate of the spectral energy density of the high-frequency 
branch. The plots of Fig. [13] (see also Fig. ^ demonstrate that the low-frequency bounds on tt do 
not forbid a larger signal at higher frequencies. 

A decrease of tt implies a suppression of the nearly scale-invariant plateau in the region z/gq < 
u < Ug. At the same time the amplitude of the spectral energy density still increases for frequencies 
larger than the frequency of the elbow (i.e. z/g). The latter trend can be simply understood since, 
at high frequency, the transfer function for the spectral energy density grows faster than the power 
spectrum of infiationary origin. For instance, in the case Wt = 1 and neglecting logarithmic corrections, 
^Gwi^y'i'o) oc z/"T+i for u ^ z/g. Now, recall that is given by Eq. (11.41) . If tt — > 0, the combination 
{ut + 1) will be much closer to 1 than in the case when, say, tt — 0.3. This aspect can be observed 
in Fig. [H where different values of rx have been reported. By decreasing the Wt from 1 to, say, 0.6 the 
extension of the nearly flat plateau gets narrower. This is also a general effect which is particularly 
evident by comparing the two curves of Fig. [13] (plot at the right). The slope of the high-frequency 
branch of the graviton energy spectrum can be easily deduced with analytic methods and it turns out 
to to be 

dlnQcw Qwt-2 
ami' 6Wt + 1 

up to logarithmic corrections. The result of Eq. (I8.15P stems from the simultaneous integration of the 
background evolution equations and of the tensor mode functions according to the techniques described 
in section [3] The semi-analytic estimate of the slope (see [75]) agrees with the results obtained by 
means of the transfer function of the spectral energy density. 

To conclude the discussion it is appropriate to elaborate on the interplay between the stiff spectra 
and the phenomenological bounds mentioned in subsection II. 6[ Let us start with millisecond pulsar 
bound of Eq. ([LTi]) . 

Assuming the maximal growth of the spectral energy density and the minimal value of z/g, i.e. z^bbn 
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Figure 14: The bounds stemming from the amount of extra-relativistic species at the epoch of the 
synthesis of hght nuclei are apphed to the rehc graviton spectra from the stiff epoch (plot at the left). 
The detectability constraint (full line in the plot at the right) stemming from the putative sensitivities 
of wide-band interferometers in their advanced version. The points corresponding to the spectral 
energy density should lie above the full lines to be potentially interesting for those instruments. 
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we will have 

hlncwil^, To) OCU, Z/ > Z/s ~ Z/bbn- (8.16) 

Since z^puisar — lO^i^bbn, Eq- (I8.16P implies that /io^Gw(^'puisar, To) — 10^^^ or even 10^^^ depending 
upon tt- But this value is always much smaller than the constraint stemming from pulsar timing 
measurements (see Eq. (I1.14p ). If either Us ^ ^hhn or Cst < 1 the value of /io^Gw(z^puisar, tq) will be 
even smalleiEll. Consequently, even in the extreme cases when the frequency of the elbow is close to 
z^bbm the spectral energy density is always much smaller than the requirement of Eq. fll.l4p . 

As noticed in the past |121l 175] . the most significant constraint on the stiff spectra stems from 
BBN. The models illustrated in Fig. [13] (plot at the right) are on the verge of saturating the bounds 
of Eqs. (Il.18p - fll.19p . This conclusion stems directly from the form of spectral energy density: the 
broad spike dominates the (total) energy density of relic gravitons which are inside the Hubble radius 
at the time of big bang nucleosynthesis. For consistency with the low-frequency determinations of the 
tensor power spectrum must be bounded from above according to the values reported, for instance, 
in Tabs. [1] and [21 Once the the value of has been selected, the constraints of Eqs. fll.18p - fll.19p can 
be imposed. From Fig. [13] a large signal is expected for uu/ — 100 Hz for 0.35 < Wt < 0.6. This range 
turns out to be compatible with the bounds of Eqs. fll.18p - fll.19p . In the opposite limit (e. g. Wt — 1) 
the spike becomes narrower, the elbow frequency augments and the signal at the interferometer scale 
diminishes. 

In Fig. [13] the energy density of the relic gravitons inside the Hubble radius at the nucleosynthesis 
epoch is reported in the case = 0.1 and for different values of S. In the plot at the left ris = 0.963 as 
implied by the WMAP 5-yr data alone. The acceptable region of the parameter space must stay below 
the horizontal lines which illustrate different values of ANi, (see Eqs. fll.18p - fll.19p ). As the scalar 
spectral index diminishes, the constraints are better satisfied since rig controls and, consequently, 

^^This conclusion follows immediately from the hierarchy between I'puisar and i-'bbn- If either I's 3> lA^hn, ^o^gw can 
only grow very little and certainly much less than required to violate the bound of Eq. ()1.14p . 
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the frequency dependence of the tensor spectral index nx (see Eq. (11.41) ) in the case ax 7^ 0. 




Figure 15: The spectral energy density is illustrated for small values of Wt and different values of S 
(plot at the left). In the plot at the right E = 0.2 and Wt = 0.6. 



8.3 Detectability prospects 

By lowering Wt, /io^gw('^, ^o) increases for u = z/lv — 0.1 kHz. This trend can be inferred from Fig. [H] 
(plot at the right) where the spectral energy density is evaluated exactly for u = z/lv- To be detectable 
by wide band interferometers the parameters of the TACDM must lie above the full lines. The region 
of low barotropic indices emerging neatly from Fig. [TH leads to spectral energy densities which are 
progressively flattening as Wt diminishes towards 1/3. Low values of Wt bring the frequency of the 
elbow, i.e. below 10~^° Hz which is unacceptable since it would mean that, during nucleosynthesis, 
the Universe was dominated by the stiff fluid. In Fig. [13] (plot at the left) the region above the 
full hne corresponds to a range of parameters for which z/g > Uhhn- in such a range a decrease of Wt 
demands an increase of E. The latter occurrence is illustrated in Fig. [15] where, at the left, Wt = 0.5. 
The full, dashed and dot-dashed curves illustrated in Fig. [15] (plot at the left) are incompatible with 
the phenomenological constraints since the frequency of the elbow is systematically smaller than z/bbn- 
Once more, this choice of parameters would contradict the bounds of Fig. [12] and would imply that the 
stiff phase is not yet finished at the BBN time. In the left plot of Fig. [15] the diamonds denote a model 
which is compatible with BBN considerations but whose signal at the frequency of interferometers 
is rather small (always three orders of magnitude larger than in the case of conventional inflationary 
models). 

The compatibility with the phenomenological constraints demands that the parameters of the 
TACDM paradigm must lie above the full lines of Fig. [TH] (plot at the left). The requirements of 
Fig. [T3] suggest, therefore, that E should be raised a bit. In this case the frequency of the elbow gets 
shifted to the right but, at the same time, the overall amplitude of the spike diminishes. The putative 
amplitude remains still much larger than the conventional inflationary signal reported, for instance, 
in Fig. [31 In Fig. [TH] the spectral energy density of the relic gravitons is illustrated as a function 
of tt for a choice of parameters which is compatible with all the bounds applicable to the stochastic 
backgrounds of the relic gravitons. The three curves refer to three different frequencies, i.e. 0.1 kHz, 
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Figure 16: The graviton energy spectrum is illustrated, in the TACDM scenario, for u = vym and as a 
function of rx (plot at the left) in the case ^ 0. In the plot at the right = 0. 



1 kHz and 10 kHz. Indeed, if the spectrum is nearly scale- invariant (as in the case o Fig. [3]) we can 
compare the potential signal with the central frequency of the window. If the signal increases with 
frequency it is interesting to plot the same curve for some significant frequencies inside the window of 
wide-band interferometers. Even if the frequency window extends from few Hz to 10 kHz the maximal 
sensitivity is in the central region and depends upon various important factors which will now be 
briefly discussed. 

To illustrate more quantitatively this point we remind the expression of the signal-to-noise ratio 
(SNR) in the context of optimal processing required for the detection of stochastic backgrounds: 



SNR' 




3Hi 



(8.17) 



{F depends upon the geometry of the two detectors and in the case of the correlation between two 
interferometers F = 2/5; T is the observation time). In Eq. fl8.17p . S^''\f) is the (one-sided) noise 
power spectrum (NFS) of the k-th {k = 1, 2) detector. The NFS contains the important informations 
concerning the noise sources (in broad terms seismic, thermal and shot noises) while 7(z/) is the 
overlap reduction function which is determined by the relative locations and orientations of the two 
detectors. In [68] Eq. fl8.17p has been used to assess the detectability prospects of gravitons coming 
from a specific model of stiff evolution with Wt = 1. At that time the various suppressions of the low- 
frequency amplitude as well as the free-streaming effects were not taken into account. Furthermore, 
the evaluation of the energy transfer function was obtained, in [70], not numerically but by matching 
of the relevant solutions. We do know, by direct comparison, that such a procedure is justified but 
intrinsically less accurate than the one proposed here. It would be interesting to apply Eq. (18.171) 
for the (more accurate) assessment of the sensitivities of different instruments to a potential signal 
stemming from the stiff age @. 



For intermediate frequencies the integral of Eq. (|8.17p is sensitive to the form of the overlap reduction function 
which depends upon the mutual position and relative orientations of the interferometers. The function ^{v) effectively 
cuts-off the integral which defines the signal to noise ratio for a typical frequency v ~ l/(2(i) where d is the separation 
between the two detectors. Since XIgw increases with frequency (at least in the case of relic gravitons from stiff ages) 
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Equation (18.171) assumes that the intrinsic noises of the detectors are stationary, Gaussian, uncor- 
related, much larger in amphtude than the gravitational strain, and statistically independent on the 
strain itself [71], [72l [73], [74] . The integral appearing in Eq. (18.170 extends over all the frequencies. 
However, the noise power spectra of the detectors are defined in a frequency interval ranging from 
few Hz to 10 kHz. In the latter window, for very small frequencies the seismic disturbances are the 
dominant source of noise. For intermediate and high frequencies the dominant sources of noise are, 
respectively, thermal and electronic (i.e. shot) noises. The wideness of the band is very important 
when cross-correlating two detectors: typically the minimal detectable h^Qcw will become smaller 
(i.e. the sensitivity will increase) by a factor 1/V AuT where Ai/ is the bandwidth and T, as already 
mentioned, is the observation time. Naively, if the minimal detectable signal (by one detector ) is 
hlQcw — 10~^, then the cross-correlation of two identical detector with overlap reduction 7(z/) = 1 
will detect hl^GW ^ 10"^° provided Au ~ 100 Hz and T ~ O(lyr) (recall that lyr = 3.15 x lO^Hz"^). 
The achievable sensitivity of a pair of wide band interferometers crucially depends upon the spectral 
slope of the theoretical energy spectrum in the operating window of the detectors. So, a fiat spectrum 
will lead to an experimental sensitivity which might not be similar to the sensitivity achievable in 
the case of a blue or violet spectra. Previous calculations [68l [69] [70] showed that, however, to get 
a reasonable idea of the potential signal it is sufficient to compare the signal with the sensitivity to 
fiat spectrum which has been reported in Eq. fll.l3p . Of course any experimental improvement in 
comparison with the values of Eq. fll.l3p will widen the detectability region by making the prospects 
of the whole discussion more rosy. 

In the TACDM paradigm the maximal signal occurs in a frequency region between the MHz and the 
GHz. This intriguing aspect led to the suggestion (68] [69] that microwave cavities [152] can be used as 
GW detectors precisely in the mentioned frequency range. Prototypes of these detectors |153j have been 
described and the possibility of further improvements in their sensitivity received recently attention 
[154] 11551 11561 1157] I158[ I159j . Different groups are now concerned with high-frequency gravitons. In 
[155] the ideas put forward in |152[I153[I154] have been developed by using electromagnetic cavities (i.e. 
static electromagnetic fields). In |156l 11571 1158] dynamical electromagnetic fields (i.e. wave guides) 
have been studied always for the purpose of detecting relic gravitons. Yet a different approach to the 
problem has been described in [159j . In [158] an interesting prototype detector was described with 
frequency of operation of the order of 100 MHz (see also ^160j ). It is not clear if, in the near future, the 
improvements in the terrestrial technologies will allow the detection of relic gravitons for frequencies, 
say, larger than the MHz. It is, however, a rather intriguing possibility. 



at most as u and since there is a in the denominator, the main contribution to the integral should occur for v < 0.1 
kHz. This argument can be explicit verified in the case of the calculations carried on in [68] and it would be interesting 
to check it also in our improved framework. 
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9 Final remarks 



The incoming score year might witness direct experimental evidences of rehc gravitons either from small 
frequency experiments or from high-frequency experiments. By low-frequency experiments we mean, 
as in the bulk of this review, the CMB experiments possibly analyzed together with the two remaining 
cosmological data sets (i.e. large-scale structure determinations of the matter power spectrum and 
type la supernova observations). By high-frequency experiments we mean the appropriately advanced 
versions of wide band interferometers such as Ligo, Virgo, Geo and Tama. 

The main observables related to relic gravitons have been reviewed in a self-contained manner and 
in the framework of the ACDM paradigm with specific attention to the complementarity between the 
low-frequency and the high-frequency branches of the relic graviton spectrum. Any model claiming a 
signal coming from high-frequency gravitons should be compatible with the ACDM paradigm in the 
low-frequency branch of the relic graviton spectrum. 

It is instructive to go back to the comparison drawn, in Fig. [1], between the electromagnetic 
spectrum and the relic graviton spectrum. The gap between the graviton frequencies explored by 
CMB experiments and the graviton frequencies probed by wide-band interferometers is of the order of 
the frequency gap between radio waves and 7-rays. 

A detection of long wavelength gravitons in CMB experiments can be direct (i.e. from the B-mode 
polarization) or indirect (i.e. from some global fit of CMB observables including the tensor contribu- 
tion). In the context of the ACDM paradigm the CMB detection of long wavelength gravitons will fix 
the overall normalization of the spectral energy density. Even in the absence of such a direct detection, 
the current upper limits on the contribution of long wavelength gravitons to CMB observables implies 
a minute signal at higher frequencies. The (hoped) sensitivities achievable by the advanced wide-band 
interferometers are still insufficient for a direct detection of high-frequency gravitons. 

The latter statement summarizes the standard lore of the problem which may well be realized. 
At the same time it might be unwise to assume (or presume) that the current success of the ACDM 
paradigm also fixes the whole thermal history of the Universe for temperatures larger than the MeV. 
The general ideas conveyed in the present review suggest that the high-frequency branch of the relic 
graviton spectrum is rather sensitive to the whole post-infiationary thermal history of the Universe. 
If the post-infiationary evolution is dominated by stiff sources, for instance, it is not impossible, as 
explicitly shown, to have positive detection of relic gravitons both at small and high frequencies even 
enforcing the current bounds on the tensor contribution to CMB observables. 

If we assume (or strongly believe) the standard lore (i.e. that relic gravitons will probably not be 
seen by wide-band detectors) it is useless to demand theoretical accuracy. For instance it is useless to 
ask what would be the effect of changing Q\ on a signal which is anyway 6 or even 7 orders of magnitude 
smaller than the most optimistic sensitivities. A positive detection of relic graviton backgrounds at 
high-frequencies would demand, however, more accurate estimates of the theoretical signal in different 
models. Absent a direct detection of relic gravitons by wide-band inteferometers, accurate theoretical 
calculations can be used to set bounds on possible deviations of the post-infiationary thermal history. 
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